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We study the dual Higgs theory for the confinement mechanism based on Quantum 

rH I Chromodynamics (QCD) in the 't Hooft abehan gauge. In the abehan gauge, QCD is 

^ ■ reduced into an abehan gauge theory including color-magnetic monopoles, which appear 

corresponding to the nontrivial homotopy group n2(SU(A^c)/U(l) ''~"^) = Z ''~^. With 

the two conjectures of abelian dominance and monopole condensation, QCD in the abelian 

gauge becomes the dual Higgs theory, and then color confinement can be understood as 

^ . one-dimensional squeezing of the color-electric flux due to the dual Meissner effect. In 

(<— ^ I the basis of the dual superconductor picture, confinement phenomena are systematically 

^^ ' studied using the lattice QCD Monte-Carlo simulation, the monopole-current dynamics 

^^ . and the dual Ginzburg-Landau (DGL) theory, an infrared effective theory of QCD. 

Q>^ , First, we study the origin of abelian dominance for the confinement force in the maxi- 

0^ I mally abelian (MA) gauge in terms of the gluon-field properties using the lattice QCD. In 

the MA gauge, the gluon-field fluctuation is maximally concentrated in the abelian sector. 

As the remarkable feature in the MA gauge, the amplitude of the off-diagonal gluon is 

Q^l strongly suppressed, and therefore the phase variable of the off-diagonal (charged) gluon 

tends to be random, according to the weakness of the constraint from the QCD action. 

Using the random-variable approximation for the charged-gluon phase variable, we find 

the perimeter law of the charged-gluon contribution to the Wilson loop and show abelian 

p\ ' dominance for the string tension in the semi-analytical manner. These theoretical results 

j^ . are also numerically confirmed using the lattice QCD simulation. 

Second, we study the QCD-monopole appearing in the abelian sector in the abelian 
gauge. The appearance of monopoles is transparently formulated in terms of the gauge 
connection, and is originated from the singular nonabelian gauge transformation to realize 
the abelian gauge. We investigate the gluon field around the monopole in the lattice QCD. 
The QCD-monopole carries a large fluctuation of the gluon field and provides a large 
abelian action of QCD. Nevertheless, QCD-monopoles can appear in QCD without large 
cost of the QCD action, due the large cancellation between the abelian and off-diagonal 
parts of the QCD action density around the monopole. We derive a simple relation between 
the confinement force and the monopole density by idealizing the monopole contribution 
to the Wilson loop. 

Third, we study the monopole-current dynamics using the infrared monopole-current 
action defined on a lattice. We adopt the local current action, considering the infrared 
screening of the inter-monopole interaction due to the dual Higgs mechanism. When the 
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monopole self-energy a is smaller then Oc = ln(2L' — 1), monopole condensation can be 
analytically shown, and we find this system being in the confinement phase from the 
Wilson loop analysis. By comparing the lattice QCD with the monopole-current system, 
the QCD vacuum is found to correspond to the monopole-condensed phase in the infrared 
scale. We consider the derivation of the DGL theory from the monopole ensemble, which 
would be essence of the QCD vacuum in the MA gauge because of abelian dominance and 
monopole dominance. 

Finally, we study the QCD phase transition at finite temperatures in the DGL the- 
ory. We formulate the effective potential at various temperatures in the imaginary-time 
formalism. Thermal effects reduce the QCD-monopole condensate and bring a first-order 
deconfinement phase transition. We find a large reduction of the self-interaction among 
QCD-monopoles and the glueball mass near the critical temperature by considering the 
temperature dependence of the self-interaction. The string tension is also calculated at 
finite temperatures. We apply also the DGL theory for the bubble formation process in 
early Universe and quark-gluon-plasma (QGP) formation process in the ultra-relativistic 
heavy-ion collision. 
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Chapter 1 
Introduction 



Quantum Chromodynamics (QCD) is the fundamental theory of the strong interac- 
tion, and describes the properties and the underlying structure of hadrons in terms 
of quarks and gluons. The QCD Lagrangian has the local SU(A'"c) symmetry and is 
written by the quark q and the gluon field A^ as 

^QCD = -^tr {G^^Cn + q i^l.D^ - m) q, (1.1) 

where (j^,y is the SU(A^c) field strength Gfj^u = j^[Dfj_,Diy] and D^ is covariant- 
derivative D^ = <9^ + ieA^^ |]1], ^, |^. In the chiral limit m -^ with Nf flavor, this 
Lagrangian has also global chiral symmetry \J{Nf)L x U[Nf)ji, although U(1)a is 
explicitly broken by the U(l)^ anomaly at the quantum level |T|. 

Due to the asymptotic freedom, the gauge-coupling constant of QCD becomes 
small in the high-energy region and the perturbative QCD provides a direct and 
systematic description of the QCD system in terms of quarks and gluons. On the 
other hand, in the low-energy region, the strong gauge-coupling nature of QCD 
leads to nonperturbative features like color confinement, dynamical chiral-symmetry 
breaking [^ |§, §[ and nontrivial topological effect by instantons [0, §, 0|, and it is 
hard to understand them directly from quarks and gluons in a perturbative manner. 
Instead of quarks and gluons, some collective or composite modes may be relevant 
degrees of freedom for the nonperturbative description in the infrared region of 
QCD. As for chiral dynamics, the pion and the sigma meson play the important 
role for the low-energy QCD, and they are included in the effective theory like the 
(non-) linear sigma model |l|, |T^, the chiral bag model |jTl], |12|] and the Nambu- 



Jona-Lasinio model |]T3|, H, where these mesons are described as composite modes 



of quarks. Here, the pion is considered to be the Nambu-Goldstone boson relating 
to spontaneous chiral-symmetry breaking and obeys the low-energy theorem and 
the current algebra [|l|. On the other hand, confinement is essentially described by 
dynamics of gluons rather than quarks. Hence, it is quite desired to extract the 
relevant collective mode from gluon for confinement phenomena. 
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In 1970's, Nambu, 't Hooft and Mandelstam proposed an interesting idea that 
quark confinement can be interpreted using the dual version of the superconductivity 



fTSl , |T6|, ^Tj (see Fig. |l.l|) . In the ordinary superconductor, Cooper-pair condensa- 
tion leads to the Meissner effect, and the magnetic flux is excluded or squeezed like 
a quasi-one-dimensional tube as the Abrikosov vortex, where the magnetic flux is 
quantized topologically. On the other hand, from the Regge trajectory of hadrons 
and the lattice QCD, the confinement force between the color-electric charge is char- 
acterized by the universal physical quantity of the string tension, and is brought by 



one-dimensional squeezing of the color-electric flux [18| in the QCD vacuum. Hence, 
the QCD vacuum can be regarded as the dual version of the superconductor based 
on above similarities on the low-dimensionalization of the quantized flux between 
charges. In this dual-superconductor picture for the QCD vacuum, as the result of 
condensation of color-magnetic monopoles, which is the dual version of the electric 
charge, the squeezing of the color-electric flux between quarks is realized by the dual 
Meissner effect. However, there are two following large gaps between QCD and the 
dual superconductor picture. 

1. This picture is based on the abelian gauge theory subject to the Maxwell- 
type equations, where electro- magnetic duality is manifest, while QCD is a 
nonabelian gauge theory. 

2. The dual-superconductor scenario requires condensation of (color-) magnetic 
monopoles as key concept, while QCD does not have such a monopole as the 
elementary degrees of freedom. 

As the connection between QCD and the dual superconductor scenario, 't Hooft 
proposed concept of the abelian gauge fixing [|l^, the partial gauge fixing which only 
remains abelian gauge degrees of freedom in QCD. By definition, the abelian gauge 
fixing reduces QCD into an abelian gauge theory, where the off-diagonal element 
of the gluon field behaves as a charged matter field and provides a color-electric 
current in terms of the residual abelian gauge symmetry. As a remarkable fact in the 
abelian gauge, color-magnetic monopole appears as topological object corresponding 
to nontrivial homotopy group H2(SU(A^c)/U(l)^''~^) = Z^'^"^. Thus, by the abelian 
gauge fixing, QCD is reduced into an abelian gauge theory including both the electric 
current j^ and the magnetic current k^, which is expected to provide a theoretical 
basis of the monopole-condensation scheme for the confinement mechanism. 

For irrelevance of the off-diagonal gluons, Ezawa and Iwazaki assumed abelian 
dominance that the only abelian gauge fields with monopoles would be essential for 
the description of nonperturbative phenomena in the low-energy region of QCD, and 
showed a possibility of monopole condensation in an infrared scale by investigating 
"energy-entropy balance" on the monopole current pO| in a similar way to the 



Kosterlitz-Thouless transition in the 1+2 dimensional superconductivity 0]. Ezawa 



and Iwazaki formulated the dual London theory as an infrared effective theory of 



Figure 1.1: Correspondence of the QCD vacuum to the superconductor. In the 
QCD vacuum, the color confinement is reahzed by monopole condensation. In the 
superconductor, the Meissner effect occurs by Cooper-pair condensation. 



QCD, and later Maedan and Suzuki reformulated it as the dual Ginzburg-Landau 
(DGL) theory in 1988 [||]. 

Furthermore, such abelian dominance and monopole condensation have been 
investigated using the lattice QCD simulation in the maximally abelian (MA) gauge 
23, 24, 25, 26 1 . The MA gauge is the abelian gauge where the diagonal component 
of gluon is maximized by the gauge transformation. In the MA gauge, physical 
information of the gauge configuration is concentrated into the diagonal components 
as much as possible. The lattice QCD studies indicate abelian dominance that the 
string tension 12^, E^, ^ and chiral condensate 12^, |2^ are almost described only by 



abelian variables in the MA gauge. In the lattice QCD, monopole dominance is also 
observed such that only the monopole part in the abelian variable contributes to the 
nonperturbative QCD in the MA gauge ||28|, ^. Thus, the lattice QCD simulations 



show strong evidence on the dual Higgs theory for the nonperturbative QCD in the 



MA gauge |31, 32 



As the result, DGL theory is expected as a reliable infrared effective theory of 
QCD. Recently, RCNP group studied the DGL theory |2^, Q, and derived the sim- 
ple formula for the string tension and also pointed out the relevant role of monopole 
condensation to chiral symmetry breaking by solving the Schwinger-Dyson equation 
with the nonperturbative gluon propagator |]33| , p4| , |35[| . This abelian dominance for 
chiral-symmetry breaking is confirmed by Miyamura and Woloshyn in more rigid 
framework of the lattice QCD p8|, ^. Considering the fact that instanton needs 
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the nonabelian component for existence, RCNP group pointed out the correlation 
between instantons and monopoles in the abehan dominant system in terms of the 
remaining large off-diagonal component around the monopole, which is the topologi- 
cal defect. The evidence on the strong correlation between instantons and monopoles 



have been observed also in the lattice QCD calculation [^ and in the analytical 



demonstration using the Polyakov-like gauge |^ ^ and the MA gauge 



Thus, the monopole seems to play the essential role for the nonperturbative QCD 
like confinement and chiral symmetry breaking and topology [^, |3^, 5T | . 



A question arises if the color is confined in the QCD vacuum by monopole con- 
densation all the time? As the superconducting state at low temperature is changed 
into the normal phase at high temperature, the QCD vacuum would also change 
from the confinement phase to the quark-gluon-plasma (QGP) phase, where the 
quark and gluon can move freely. This phase transition is called QCD phase transi- 
tion, and becomes one of the most important subject related to various fields such 
as the early Universe and relativistic heavy-ion collision. According to the big bang 
scenario, which is a successful model for cosmology, the quark gluon plasma phase 
at high temperature is changed into the hadron phase 10~^ second after the big 
bang. At the Brookhaven National Laboratory in the United States, some physi- 
cists are trying to form the quark gluon plasma as a new matter by colliding high 
energy heavy- ions using RHIC. RCNP group also has studied the QCD vacuum 
at finite temperature in terms of confinement-deconfinement phase transition and 
chiral phase transition using the DGL theory. 

In this way, the study of color confinement phenomena based on the dual Higgs 
picture is divided into two categories in terms of the method of approach. 

1. Study with the lattice QCD, which is an useful method for the direct calcula- 
tion of the QCD partition functional Zqcb [^, |3 . 



2. Study with the infrared effective theory, the DGL theory ^, ^, which con- 
sists of the essential degrees of the freedom for infrared phenomena. 

In the lattice formalism, space-time coordinate is discretized with the lattice 
spacing a, and the theory is described by the link variable U^{s) = e^"■^'^^^■i^) ^ which 
corresponds to the line integral Pexp{i J^'^'^ dx^eA^{x)} along the link. The lattice 
QCD partition functional Z in the Euclidean metric is given as 

rf[/„e-^^[^-] (1.2) 



'M 



where /? = J^ is the control parameter related to the lattice spacing a |^ (Ap- 
pendix A.l). Here, e(a) is the QCD running coupling constant. The standard lattice 
action is given as 

S= E [1 - ^tr{n,,(s) + nt^(s)}], (1.3) 



where □^;^(s) is plaquette defined as 

0,uis) ^ U,{s)U,{s + fi)Ul{s + 0)Ul{s). (1.4) 

In the hmit a — >■ 0, i.e. (3 -^ oo, and the lattice action (3S becomes the QCD action 
S'qcd in the continuum theory. The gauge configuration of QCD is generated on 
the lattice using the Monte Carlo method (Appendix A. 2). In the lattice QCD, the 
abelian gauge fixing can be performed after the generation of gauge configurations, 
and the abelian link variable is extracted by neglecting the off-diagonal part, which 
is called as the abelian projection. In the lattice formalism, the abelian link variable 
can be separated numerically into the photon part and the monopole part corre- 
sponding to the separation of the electric current and the monopole current as will 
be shown in Chapter ^ The dual-superconductor scenario for confinement has been 
examined in the lattice QCD by measurements of the abelian variable, monopole 
and so on in the abelian gauge. 

On the other hand, the DGL theory is derived from the gluon sector of QCD, 
and is composed of the dual gauge field i3^ and the monopole field x ^^ the pure 
gauge case. The Lagrangian is expressed as 

^DGL = tr£ 

£ = -^id^B, - d,B,f + [b,, x]n^M> X] - KX^X - v'f, (1.5) 

where Z)^ = 5^ + igB^ is the dual covariant derivative. Imposing the abelian gauge 
fixing on QCD, the monopole appears as the line-like object in 4-dimensional space. 
The monopole field is derived by summing all the paths of the monopole trajectories 
and monopole-field interaction is introduced taking the lattice result "monopole 
condensation" into the consideration. Here, the off-diagonal component is neglected 
due to the lattice QCD result "abelian dominance", and the dual gauge field B^ 
is used instead of the abelian gauge field A^ adopting the Zwanziger formalism in 
order to describe gluon dynamics without the singularity in the gauge field. 

The DGL theory provides an useful method of studying the various confinement 
phenomena such as inter-quark potential, hadron fiux-tube system and the phase 
transition, since it gives not only just the numerical results on the various quantities 
but also their reasons. This is largely different from the lattice simulation. The 
quark-antiquark potential arises from the strong correlation between two quarks in 
the infrared region, which is revealed by the DGL gluon propagator with double 
pole. The hadron flux is constructed by a massive dual gauge fleld, whose mass 
is obtained by the dual Higgs mechanism of monopole condensation. However, in 
the process of the construction of DGL theory, abelian dominance and monopole 
condensation is assumed, and its origin is not clear. Namely, we can not answer 
the question what feature of the monopole degrees of freedom is important for such 
conflnement phenomena or where the effect of nonabelian feature appears. 



Figure 1.2: Comparison among QCD, abelian projected QCD (AP-QCD) and QED 
in terms of the gauge symmetry and fundamental degrees of freedom. (A^^ = 2) 

In this paper, we try to understand the confinement mechanism based on the 
dual Higgs mechanism using both methods, the lattice QCD simulation and the DGL 
theory. In the first half of this paper, we study the origin of abelian dominance and 
monopole condensation in terms of the gluon configuration using the lattice QCD, 
and in the second half we apply the DGL theory to the confinement phenomena 
such as QCD phase transition and multi-hadron flux tube system. 

The point at issue in the first half is the following three subjects. 

1. What is the origin of abelian dominance for infrared quantities like the string 
tension in the MA gauge? 

2. Why does monopole appear in abelian projected QCD(AP-QCD), although 
AP-QCD is an abelian gauge theory like QED? 

3. What is the role of monopoles for the confinement phenomena? What is the 
role of the off-diagonal component of gluon in QCD in the MA gauge? 

In the MA gauge, AP-QCD neglecting the off-diagonal gluon component almost 
reproduces essence of the nonperturbative QCD, although AP-QCD is an abelian 
gauge theory like QED. One may speculate that the strong-coupling nature leads 



to the similarity between AP-QCD and QCD, because the gauge couphng e in AP- 
QCD [Q is the same as that in QCD in the lattice simulation. However, the 
strong-coupling nature would not be enough to explain the nonperturbative feature, 
because, if monopoles are eliminated from AP-QCD, nonperturbative features are 
lost in the remaining system called as photon part, although the gauge couphng e 
is same as that in QCD. For further understanding, let us compare the theoretical 
structure of QCD, AP-QCD and QED in terms of the gauge symmetry and the fun- 
damental degrees of freedom as shown in Fig.1.2. As for the interaction, the linear 
confinement potential arises both in QCD and in AP-QCD, while only the Coulomb 
potential appears in QED. On the symmetry, QCD has a nonabelian gauge sym- 
metry, while both AP-QCD and QED have abelian gauge symmetry. The obvious 
difference between QCD and QED is existence of off-diagonal gluons in QCD. On the 
other hand, the difference between AP-QCD and QED is existence of the monopole, 
since the magnetic monopole does not exist in QED because of the Bianchi iden- 
tity. This indicates the close relation between monopoles and off-diagonal gluons. 
In particular, off-diagonal gluon components play a crucial role for existence of the 
monopole [^, and the monopole itself is expected to play an alternative role of 
off-diagonal gluons for confinement. 

In Chapter 0, we review the abelian gauge fixing in line with 't Hooft to discuss 
the confinement phenomena in terms of the monopoles based on the dual super- 
conductor picture. We show the gauge invariance condition in the abehan gauge. 
As the abelian gauge fixing, we introduce the maximally abelian gauge, which is 
considered to be the best abelian gauge for the infrared physics according to recent 
lattice QCD studies. In addition, we generalize the maximally abelian gauge fixing. 

In Chapter ^, we investigate the origin of abelian dominance in the MA gauge in 
the SU(2) gauge theory. We introduce the abelian projection rate, which is defined 
as the overlapping factor between SU(2) link variable and abelian link variable, 
and investigate abelian dominance for the abelian link variable. We study abelian 
dominance for the Wilson loops in terms of by approximating the off-diagonal angle 
variable as a random variable. Using the U(l)3 Landau gauge, we study the abelian 
gauge field and abelian field strength directly in the MA gauge. Then, we compare 
the abelian gauge field in the SU(2) Landau gauge, where the gauge field is fixed 
most continuously and study the feature of the MA gauge in terms of the gluon field 
fluctuation. 

In Chapter §, we investigate the mechanism of the appearance of monopole in 
AP-QCD both in the continuum theory and in the lattice QCD formalism. We show 
the appearance of monopole using the covariant derivative in the abelian sector of 
QCD, clarifying the role of the off-diagonal components of QCD. 

In Chapter ^, we investigate the relations between monopoles and gauge-field 
fluctuations in the MA gauge, measuring the probability distribution of gluons 
around the monopole. We show the distribution of the action density on the SU(2), 
U(l)3 and off-diagonal part around the monopoles and consider the appearance of 
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the monopole in terms of the action density. 

In Chapter ^, we study extraction of the monopole degrees of freedom from 
11(1)3 gauge theory. The abelian gauge field is decomposed into the monopole and 
photon parts. We investigate the properties on the magnetic and electric currents, 
action density, field variable itself in both parts. We show the scaling properties on 
the monopole current and the Dirac sheet and obtain the good scaling on variables 
related to the Dirac sheet. Furthermore, we obtain the simple relation between the 
string tension and monopole current density. 

In Chapter |^, we study monopole condensation in the QCD vacuum by com- 
paring to the monopole- current system. We first generate the monopole- current 
system on the lattice using a simple monopole current action, and study the role of 
monopole current to confinement. 

In the second part of this paper, making the best use of the DGL theory, we 
investigate the QCD phase transition using effective potential formalism, and multi- 
flux-tube system, which cannot be studied by the lattice QCD simulation. 

In Chapter p, we first review the derivation of the DGL theory starting from the 
QCD Lagrangian. We derive the monopole field by summing all of the monopole 
trajectories, and the dual gauge field is introduced using the Zwanziger formalism 
in order to describe the gluon dynamics without the singularities originated from 
the monopole. Then, using the DGL theory, we show the dual Meissner effect 
by monopole condensation and the structure of the flux tube as the dual version 
of the Abrikosov vortex in the superconductor. We also demonstrate the quark 
confinement potential using the DGL gluon propagator including nonperturbative 
effect. Finally, we discuss the asymptotic behavior in terms of the DGL theory. 

In Chapter ^, we consider the QCD vacuum at finite temperature in the DGL 
theory. We formulate the effective potential at various temperatures by introducing 
the quadratic source term, which is a new useful method to obtain the effective 
potential in the negative- curvature region. We find the thermal effects reduce the 
monopole condensate and bring a first-order phase transition. 



In Chapter |10|, we apply the DGL theory to the multi-fiux tube system. We 
formulate this system by regarding it as the system of two flux tubes penetrating 
through a two-dimensional sphere surface. We find the multi-fiux-tube configuration 
becomes uniform above some critical fiux-tube density. 



Chapter 2 

Abelian Gauge Fixing 



In infrared QCD, there appear interesting nonperturbative phenomena such as color 
confinement and chiral symmetry breaking due to the strong couphng nature. At 
the same time, however, the large gauge- coupling constant leads to breaking down 
of the perturbative technique. As the result, it becomes difficult to treat the in- 
frared phenomena analytically. We have to use an effective model, which includes 
essence. Otherwise, we can perform the partition functional of QCD directly using 
the huge supercomputer. Historically, the dual superconductor picture was proposed 
by Nambu and 't Hooft more than 20 years ago to understand the confinement 
mechanism. The QCD vacuum is regarded as a dual version of the superconductor 
and color confinement is understood as the exclusion of the color-electric field by 
monopole condensation. Later, 't Hooft and Ezawa-Iwazaki showed that QCD is 
reduced into the U(l) gauge theory with monopoles by the abelian gauge fixing. 
If "abelian dominance" and "monopole condensation" occurs in the QCD vacuum, 
color confinement would be realized through the dual Higgs mechanism. Recently, 
these assumptions are supposed by the Monte Carlo simulation of the lattice QCD. 
In this chapter, we study the abelian gauge fixing in QCD in terms of the residual 
gauge symmetry. In the abelian gauge, the SU(A^c) gauge theory is reduced into the 
U(l)^'=~^ gauge theory |]T^ ^, and confinement phenomena can be studied by the 



dual abelian Higgs theory. In this gauge, the diagonal gluon component remains 
to be a U(l)^'=~^ gauge field, and the off-diagonal gluon component behaves as 
a charged matter field and provides the color electric current j^ in terms of the 
residual U(l)^'=~^ gauge symmetry. In the abelian gauge, color-magnetic monopoles 
also appear as topological defects provided by the singular gauge transformation. 
Thus, QCD in the abelian gauge is reduced into an abelian gauge theory including 
both the electric and monopole currents, which is described by the extended Maxwell 
equation with the magnetic current. For simplicity, we concentrate ourselves on the 
Nr = 2 case hereafter. 
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2.1 Residual Symmetry and Gauge Invariance Con- 
dition in the Abelian Gauge 

The abelian gauge fixing, the partial gauge fixing which remains the abelian gauge 
symmetry, is realized by the diagonalization of a suitable SU(A^c) -gauge dependent 
variable as <l>[y4^(a;)] e su{Nc) by the SU(A^c) gauge transformation. In the abelian 
gauge, $[A^(a:)] plays the similar role of the Higgs field on the determination of the 
gauge fixing. 

For an hermite operator $[A^(x)] which obeys the adjoint transformation, $(x) 
is transformed as 

= diag(Ai(a;),---,A^^(a;)), (2.1) 

using a suitable gauge function Q{x) = exp{z^'^(x)T"} G SU(A^c)- Here, each di- 
agonal component A* (^=l,■ ■ ^A^c) is to be real for the hermite operator $[y4^j(x)]. 
In the abehan gauge, the SU(A''c) gauge symmetry is reduced into the 11(1)^'=^^ 
gauge symmetry corresponding to the gauge-fixing ambiguity. The operator $(x) 
is diagonalized to H ■ ^diagi^) also by the gauge function Q'^{x) = uj{x)Q{x) with 
uj{x) = exp{-iH- 0{x)) G U(l)^=-\ 

$(X) ^ n'^{xMx)n^^{x) = UJ{X)H ■ $d^ag{xWix) = H ■ $d^ag{x), (2.2) 

and therefore 11(1)^'=^^ abelian gauge symmetry remains in the abelian gauge. 

In the abelian gauge, there also remains the global Weyl symmetry as a "relic" 
of the nonabelian theory [^^ ^. Here, the Weyl symmetry corresponds to the 



subgroup of SU(iVc) relating to the permutation of the basis in the fundamental 
representation. Then, the Weyl group is expressed as the permutation group Pjy^ 
including NJ elements. For simplicity, let us consider the Nc = 2 case. For SU(2) 
QCD, the Weyl symmetry corresponds to the interchange of the SU(2)-quark color, 
1+) = (q) and |— ) = (5), in the fundamental representation. The global Weyl 
transformation is expressed by the global gauge function, 

^ ^ e^^^^°^"+^^'"^">" = z{riCos« + r2sina} 



= ^[^^a ) GP2CSU(2) (2.3) 

with an arbitrary constant a G R. By the global Weyl transformation W, the 
SU(2)-quark color is interchanged as Vr|+) = ie*°|— ) and W^|— ) = ie~*"|+) except 
for the global phase factor. This global Weyl symmetry remains in the abelian 
gauge, because the operator $(x) is also diagonalized by using Q^{x) = WQ{x), 

$(X) ^ n^(x)$(x)n^t(3.) = W^d^ag{x)^W^ = -^d^ag{x)^. (2.4) 
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Here, the sign of ^diag{x), or the order of the diagonal component A*(a;), is globally 
changed by the Weyl transformation. It is noted that the sign of the U(l)3 gauge 
field Afj, = A^-^ is globally changed under the Weyl transformation, 

A,^A^ = WAl^W^ = -4| = -A,. (2.5) 

Therefore, all the sign of the abelian field strength, electric and magnetic charges 
are also globally changed: 



JfJ. — '^ ''afi ^ Ji_i J/1) 

k^ = d^^T^^-^k";^ = -k^. (2.6) 

In the abelian gauge, the absolute signs of the electric and the magnetic charges are 
settled, only when the Weyl symmetry is fixed by the additional condition. When 
$[^^(a;)] obeys the adjoint-type gauge transformation like the nonabelian Higgs 
field, the global Weyl symmetry can be easily fixed by imposing the additional 
gauge-fixing condition as ^diagi.^) > for SU(2), or the ordering condition of the 
diagonal components A* in H -^diag as A^ > .. > A^'^ for the SU(A^c) case. As for the 
appearance of monopoles in the abelian gauge, the global Weyl symmetry P^v^ is 
not relevant, because the nontriviality of the homotopy group is not affected by the 
global Weyl symmetry. However, the definition of the magnetic monopole charge. 



which is expressed by the nontrivial dual root of SU(A^c)duai [19, is globally changed 
by the Weyl transformation. 

Now, we consider the abelian gauge fixing in terms of the coset space of the fixed 
gauge symmetry. The abelian gauge fixing is a sort of the partial gauge fixing which 
reduces the gauge group G =SU(A^c) local of the system into its subgroup H =U(1) 
io^~i^(xP|r° ^) including the maximally torus subgroup of G. In other words, the 
abelian gauge fixing freezes the gauge symmetry relating to the coset space G/H, and 
hence the representative gauge function Q which brings the abelian gauge belongs 
to the coset space G/H: fl G G/H. In fact, Q G G/H is uniquely determined 
without the ambiguity on the residual symmetry H, if the additional condition on 
H is imposed for Q. 

However, such a partial gauge fixing makes the total gauge invariance unclear. 
Here, let us consider the SU(A^c) gauge-invariance condition on the operator defined 
in the abelian gauge |^. To begin with, we investigate the gauge-transformation 
property of the gauge function Q G G/H which brings the abelian gauge (see 
Fig. |2.1| ). For simplicity, the operator $ to be diagonalized is assumed to obey 
the adjoint gauge transformation as $ -^ $^ = g^g^ with ^g G G. After the 
gauge transformation by ^g E G, Q^ G G/H is defined so as to diagonalize $^ as 
Q^^^{Q^y = ^diagi and hence the gauge function Vt^ G G/H which realizes the 
abelian gauge is transformed as 

Vl^Q.3 = h[g]^g^ (2.7) 
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under arbitrary SU(A^c) gauge transformation hj g E G. Here, h[g] E H is chosen 
so as to make Q^ belong to G/H, i.e., Q^ G G/H. (If the additional condition on H 
is imposed to specify Q G G/H, Qg"^ does not satisfy it in general.) This is similar 
to the argument on the hidden local symmetry |T0[ in the nonlinear representation. 



In general, the gauge function f2 G G/H is transformed nonlinearly by the gauge 
function g due to h[g] G H. Thus, the gauge-transformation property on the gauge 
function Q G G/H becomes nontrivial in the partial gauge fixing. 

Owing to the nontrivial transformation ( [2. 71 ) of Q eG/H, any operator O^ de- 
fined in the abelian gauge is found to be transformed as O^ — >■ (O^)'^'^' by the 
SU(iVc) gauge transformation of ^g G G. We demonstrate this for the gluon field 
A^ = Q{A^ + ^9^)r2"'' in the abelian gauge. By the gauge transformation of ^g G G, 
A^ is transformed as 



Here, we have used 

{A^;y^ = g,{A^; + -d,)gl = {g2g,){A, + -d,){g2giy = {A,r^^ (2.9) 



-d^,)g2 = {g2gi){Af, + -i 

e e 



for the successive gauge transformation by gi, g2 E SU (Nc). Similarly, the operator 
O^ defined in the abelian gauge is transformed by ^5^ G G as 

qQ _^ ^Qg^w ^ n909Q9l = h[g]ng^ ■ gOg^ ■ gQ^h^lg] 

= h[g]nOn^h^[g] = h[g]0%^[g] = {Oy^^\ (2.10) 



as shown in Fig. |2.1| . Here, O is assumed to obey the adjoint transformation as 
O^ = gOg'^ for simplicity. 

Thus, arbitrary SU(A^c) gauge transformation by (? G G is mapped into the partial 
gauge transformation by h[g\ G H for the operator O^ defined in the abelian gauge, 
and O^ transforms nonlinearly as O^ -^ i^Q^'^h[g] |^y ^^iq SU(A'"c) gauge transforma- 
tion g. If the operator O^ is if-invariant, one gets (O^)'^'^^ = O^ for any h[g] G H, 
and hence O^ is also G- invariant or total SU (iVc) gauge invariant, because O^ is 
transformed into {^O^Y^^^^ = O^ by ^g G G. Thus, we find a useful criterion on the 
SU(A'"c) gauge invariance of the operator defined in the abelian gauge ^^: If the 



operator O^ defined in the abelian gauge is if-invariant, O^ is also invariant under 
the whole gauge transformation of G. 

Here, let us consider the application of this criterion to the effective theory of 
QCD in the abelian gauge, the dual Ginzburg-Landau (DGL) theory p2|, |33 



see 



Chapter |^). In the DGL theory, the local U(l)^'=~^ and the global Weyl symmetries 
remain, and the dual gauge field B^ and the monopole field Xa ['^=1) " ■ ■> ^Nc^Nc—l)] 
are the relevant modes for infrared physics. Although Bf^ is invariant under the local 
transformation of U(l) ''" C SU(A^c), B^ = B^- H is variant under the global Weyl 
transformation, and therefore B^j^ is SU(iVc) -gauge dependent object and does not 
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$ 3 ^ G ^ ^g 



n e G/H 



*^diag 



VLS = h[g\VLg-^ e G/H 




n^ ^ h[g]ng-^ 



Q 



Figure 2.1: The gauge transformation property of $ and gauge function ^2 G G/H. 
(a) After the gauge transformation by "^g G G, the operator $^ is diagonalized by 
the gauge function fl^ = hlglQg'^ G G/H. (b) The gauge transformation property 
of the operator O^ defined in the abelian gauge. If O^ is iZ-invariant, O^ is found 
to be invariant under the whole gauge transformation of G. 

appear in the real world alone. As for the monopole field, there exists one Weyl- 
invariant combination of the monopole field fluctuation, x — So Xa [0 ? which 
is also U(l)^'=~^-invariant. Therefore, the monopole fluctuation x is completely 
residual-gauge invariant in the abehan gauge, so that x is SU(A^c)-gauge invariant 
and is expected to appear as a scalar glueball with J^'-' = 0++, like the Higgs particle 
in the standard model. 



2.2 Maximally Abelian Gauge 

The abelian gauge has some arbitrariness corresponding to the choice of the op- 
erator $ to be diagonalized. As the typical abelian gauge, the maximally abelian 
(MA) gauge, the Polyakov gauge and the F12 gauge have been tested on the dual 
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superconductor scenario for the nonperturbative QCD. Recent lattice QCD studies 
show that infrared phenomena such as confinement properties and chiral symmetry 
breaking are almost reproduced in the MA gauge p3|, |2^, ^, ^ ^ |2^, |30|. In the 



SU(2) lattice formalism, the MA gauge is defined so as to maximize 



S,fl 



S,fJ, 

= 2Y^[l-2{U;,{sr + Ul{sf}] (2.11) 

S,fJ. 

by the SU(2) gauge transformation (Appendix B). Here, the link variable Ufj_{s) = 
f/°(s) +iT^U^{s) G SU(2) with f/°(s), U^Ks) G R relates to the (continuum) gluon 
field A^ = v4^T" G su{2) as U^{s) = e^^'^^f^'*), where e denotes the QCD gauge 
coupling and a the lattice spacing. In the MA gauge, the absolute value of off- 
diagonal components, U}^{s) and [/"^{s), are forced to be small. In the continuum 
limit a ^ 0, the link variable reads f/^(s) = e*"''^^^'*) = 1 + iaeA^{s) + 0{a'^), and 
hence the MA gauge is found to minimize the functional, 

Rck[A,] ^ le' I d^xiAlixf + Alix)'} = e' j d'xAl{x)A-{x), (2.12) 

with A'^{x) = -j^{Aj^{x) ±iAf^{x)}. Thus, in the MA gauge, the off-diagonal gluon 
component is globally forced to be small by the gauge transformation, and hence the 
QCD system is expected to be describable only by its diagonal part approximately. 
The MA gauge is a sort of the abelian gauge which diagonalizes the hermite 
operator 

m,is)]^J2U±,is)TsUUs). (2.13) 

Here, we use the convenient notation U^^^^s) = Uj^{s — jl) in this paper. In the 
continuum limit, the condition of the MA gauge becomes '^{id^ ± eA^)^^ = 0. 

This condition can be regarded as the maximal decoupling condition between the 
abelian gauge sector and the charged gluon sector. 

In the MA gauge, $(s) is diagonalized as $ma('S) = ^diag{s)T^ with ^diag{.s) G 
R, and there remain the local U(l)3 symmetry and the global Weyl symmetry |]46| . 
As a remarkable fact, $(s) does not obey the adjoint transformation in the MA 
gauge, and the sign of ^diag{,s) is not changed by the Weyl transformation by W in 
Eq.(i3), 

*Ma(s) = ^diag{s)Tz 

^ <a(s) = E WU^,{s)W^T,WUi^{s)W^ 

= - E WU^^{s)r^Ui^{s)W^ = -W^d^ag{s)T:,W^ = ^d^ag{sW (2.14) 
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Thus, the Weyl symmetry is not fixed in the MA gauge by the simple ordering con- 
dition as ^diag > 0, unlike the adjoint case. We find the gauge invariance condition 
on the operator O^ defined in the MA gauge: if O^ is invariant both under the 
local U(l)^'=~^ gauge transformation and the global Weyl transformation, O^ is also 
invariant under the SU(iVc) gauge transformation. 

In the continuum SU(A^c) QCD, it is more fundamental and convenient to define 
the MA gauge fixing by way of the SU(A^c)-covariant derivative operator D^ = 
d^ + ieAf^, where d^ is the derivative operator satisfying [d^,f{x)] = d^f{x). The 
MA gauge is defined so as to make SU(A^c) -gauge connection D^ = d^ + ieA'^T'' 
close to U(l)^'=~^-gauge connection D^^^^ = <9^ + ieAf^ ■ H by minimizing 

= e' d'xj:crC^,a-(3tT{EiE,) = - d'x ^ jC^^, (2.15) 



which expresses the total amount of the off-diagonal gluon component. Here, we 

Na{Nc-l) 

have used the Cartan decomposition, A^ = A^^T" = A^ ■ H + ^ C')^-E"; 

Q = l 

H = (T3, Tg, ■ ■ ■ , T/v2_i) is the Cartan subalgebra, and E'^{a = 1, 2, ■ ■ ■ , N^ — Nc) 
denotes the raising or lowering operator. In this definition with £)^, the gauge 
transformation property of Rch becomes quite clear, because the SU(A^c) covariant 
derivative D^ obeys the simple adjoint gauge transformation, D^ -^ QD^Q\ with 
the SU(iVc) gauge function ^2 G SU(A^c)- By the SU(A^c) gauge transformation, i?ch 
is transformed as 



Rch -^ Rch = j d^x tr {pD^n\ H]^[nD^n\ H]'^ 

d^x tr ([b^, nWn\^[b^, nWn]) , (2.16) 



and hence the residual symmetry corresponding to the invariance of -Rch is found 
to be U(l)i^=-i^ X Pgtbai cSU(iV,)iocai, where Pg^^^ai denotes the global Weyl group 
relating to the permutation of the Nc basis in the fundamental representation. In 
fact, one finds lu^Hlu = H ioi uj = e''^^''>^ G U(l)(^^;:i\ and the global Weyl 
transformation by W^ G -Pgiobai only exchanges the permutation of the nontrivial 
root dj and never changes -Rch- In the MA gauge, arbitrary gauge transformation 
hj^Q E G is to increase -Rch as -R^ > -Rch- Considering arbitrary infinitesimal gauge 
transformation f2 = e*^ ~ 1 + ie with ^e Gsu{Nc), one finds Q.'^Hi^ ~ -/? + i[H,€] 
and 



< ~ Rch + 2ij d^xti [[b^, [H, e]] [b^, H]] 

= Rch + 2tJ d^xti {e[H, [b^, [b^, H]]]) . (2.17) 
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For the Nc=2 case, the MA gauge extremum-condition of R^^ on "^e Esu{2) provides 

[t3,[D^,[D^,T3]]] = 0, (2.18) 

which leads to I]^(*<9^ i ^^fi)^^ = 0- Thus, the operator $ to be diagonahzed in 
the MA gauge is found to be 

^A^] = [D^,[D^,Ts]] (2.19) 

in the continuum theory. Here, $[^^] is hermite as ^^[A^] = $[A/x] because of 
Z)jj = —D^, and hence the diagonal elements of ^[A^] should be real. 

In the commutator form, the diagonal part of the variable 0[A^(x)] is expressed 

as iH 

0^ = d-[H,[H,d]]. (2.20) 



For the covariant derivative operator, one finds 

Df = D,- [H, [H, D,]] =d^ + leA^ix) ■ H (2.21) 

with An{x) = Af^{x) ■ H + C'^{x)E°'. Then, the abelian projection, D^ -^ Dj^ , is 

expressed by the simple replacement as A^{x) G su{Nc) — > A^{x) = Af^{x) ■ H e 
m(1)^^-i. 

2.3 Generalization of the Maximally Abelian Gauge 

In the MA gauge, i?ch[^/j] in Eq. (|2.15|) is forced to be reduced by the MA gauge 
transformation fiuAix) G G/H [^, and therefore the gluon field A^{x) is maximally 
arranged in the diagonal direction H in the internal SU(A''c) color space. In the 
definition of the MA gauge, H is the specific color-direction, since H explicitly 
appears in the MA gauge-fixing condition with i?ch[^/^]- On this point of view, 
the MA gauge can be called as the "maximally diagonal gauge". However, for the 
extraction of the abelian gauge theory from the nonabelian theory, we need not take 
the specific direction as H in the internal color-space, although the system becomes 
transparent when the specific color-direction as H is introduced on the maximal 
arrangement of the gluon field A^{x). 

In this section, we consider the generalization of the framework of the MA gauge 
and the abelian projection, without explicit use of the specific direction H in the 
internal color-space on the gauge fixing [^]. Instead of the special color-direction H, 



we introduce the "Cartan frame field" (f){x) = {(f)i{x),(f)2{x),- ■ ■ ,(f)Nc~i{x)), where 
0j(x) = 0"(x)T" (0"(x) G R) commutes each other as [0j(x), 0j(x)] = 0, and satisfy 
the orthonormality condition 2tr(0j(a;)0j(x)) = Y.a=i^ '^ii^)'P'ji^) — ^ij- ^^ ^^^^ 
point x^, (j){x) forms the Cartan sub-algebra, and can be expressed as 

0(a;) = nl.{x)Hnc{x) (2.22) 
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using flc{x) £ G/H. For the fixed Cartan frame field (/)(x), we define the generalized 
maximally abelian (GMA) gauge so as to minimize the functional 

i?0ch[^^] = j d^xti[b^J{x)]'^[b^J{x)] (2.23) 

by the SU(A^c) gauge transformation. Here, the Cartan frame field 0(x) is defined at 
each x^ independent of the gluon field like H, and never changes under the SU(A^c) 
gauge transformation. For the special case of 0(x) = H, the GMA gauge returns 
to the usual MA gauge. In the GMA gauge, the SU(iVc) covariant derivative D^ is 
maximally arranged to be "parallel" to the (^(a;)-direction in the internal color-space 
using the SU(A^c) gauge transformation. 

In the GMA gauge, the gauge symmetry is reduced from SU(A^c) into U(l)^''~^, 
and the generalized AP-QCD leads to the monopole in the similar manner to the MA 
gauge. In the GMA gauge, the remaining U(l)^''~^ gauge symmetry corresponds to 
the invariance of R^c\i[-^i^ under the 11(1)^''"^ gauge transformation by 

uj^{x) = e*^^(">^» e U{l)f-\ x{x) e R^^-\ (2.24) 

In fact, using lu'^{x)(I){x)lUij,{x) = 0(x), U(l)^''^^ invariance of -R^ch[^/i] is easily 
confirmed as 



iR^ch[A^]r = I d^xtT[uj{x)^b^ujl{x)J{x)]^[uj{x)^b^ujl{x)J{x)] (2.25) 

d^xtT[b^,ujl{x)(j){x)uj^{x)]^[bf„ujl{x)(j){x)uj^{x)] = R^ch[A^,]. 



There also remains the global Weyl symmetry Pjv^ similarly in the usual MA gauge, 
although the gauge function takes a complicated from. 

Here, we consider the generalized abelian projection to 0(x)-direction. Similar 
to the "diagonal part" in Eg. ( p. 201) , we define the "0(a;)-projection" of the operator 



0{x) as 

0<^(x) = d{x) - [0(x), [0(x), d{x)]], (2.26) 

using the commutation relation. For the SU(A^c) covariant derivative operator D^ = 
d^ + leA^, its 0(x)-projection is defined as 

bi^b,- [0(x), [${x),b,]] = d, + teAf,{x) + [^{x), d,${x)] (2.27) 

with^^(x) = A'f^{x)-(f) = 2tT{(f){x)A^{x))-(j){x). Here, the nontrivial term [0(s), (9^0(x) 
appears in Dj^ owing to the x-dependence of the Cartan-frame field (j){x). The 
U(l)^''~^ gauge field is defined as the difference between £)^ and d^, 

Af,{x) ^ -{Di - d,) = Af^ix) + ^[$(x),dj{x)] e suiN,). (2.28) 

6c- to 
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Here, Af^i^x) includes both the (/)(x)-component Af^i^x) = 2tT{A^{x)(j){x)) ■ (j){x) and 
the non-(/)(x)-component r^[0(a;), 9^0(x)], because [0(x), (9^0(a;)] does not include 
0(x)-component as tr f0j(x)[0(x), 9^0(x)] j = 0. Here, A'^^{x) is the image of A'^^{x) 
mapped into the U(l). ''"^-manifold. The generalized abelian projection for the vari- 
able 0[A^(x)] is defined via the two successive mapping, 0[A^(x)] — > 0[yl^(x)] — > 
Oap = 2tr{(j){x)0[Af,{x)]), after the GMA gauge fixing. 

Under the U(l)^''~^ abelian gauge transformation by ^^(^(x) = e^'^^^'^''^^^'^ E U(l)^''~\ 
Afi^ix) or A'l^{x) behaves as the U(l)^'=~^ abelian gauge field, 

Atix) -^ Mj(x))- = Aiix) + -^d,U^) ■ 0(x). (2.29) 

The abelian field-strength matrix is defined as 

to 

= d,At{x) - d^Aiix) + ze[Ai{x),At{x)], (2.30) 

which generally includes the non-0-component as well as Af^lx). The 0-component 
of ^fyix) is the image of T'jlyix) projected into the U(l)^'=~^ gauge manifold, and is 
observed as the "real abelian field-strength" in the abelian-projected gauge theory. 
The explicit form of F^^(a;) is derived as 



F^ix) ^ 2tT(^Ux)${x 



= d,At{x)-d,A^^{x) + -ti{${x)[d,U^),dM^)]) (2.31) 

to 

= d,At{x) - d^A^x) + -Jabcrd.cl^ldA'i, (2.32) 

where the last term breaks the abelian Bianchi identity and provides the monopole 
current. The magnetic monopole current is derived as 



1 



X 



fi \""/ ^ ajj, 



d'^^FUx) = —e,^p,fatcd"r{x)d^<P-{x)d'^<P'i{^), (2.33) 



which is the topological current induced by 0. Hence, the monopole appears from 
the center of the hedgehog configuration of as shown in Fig. [4.1| in the SU(2) case. 
Next, we investigate the properties of the GMA gauge function ficuAix), which 
brings the GMA gauge. Here, ficMAix) is a complicated function of A^{x) and is 
expressed by an element of the coset space G/H = SU(iVc)/{U(l)^''~^ x Weyl} as 
the representative element because of the residual gauge symmetry. For instance, 
we impose here 

tr(fiGMA(x)0(a;)) = (2.34) 
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for the selection of (^gma ^ G/ H. Similarly to the MA gauge function [^, fiGMA[^A«] 
obeys the nonlinear transformation as 

nGMA{x)eGlH^{nGMA{x)f' = d^{x)nGMA{x)V\x) G G/H (2.35) 

by the SU(A^c) gauge transformation with V{x) E G. Here, d^{x) e H = U(l)^"~^ 
X Weyl appears to keep (I^gma)^ belonging to G/H. Therefore, the gluon field 
^GMA ^ f2gj^^(^^ + i^^A')^GMA ^ 5' i^ the GMA gauge is transformed as 



'gmaI^J 



to 

= (f{x){Af'^ + -d,)(f\x) = {Af'^Y' (2.36) 



le 

GMA\dy 



by the SU(A^c) gauge transformation. As a remarkable feature, the SU(A^c) gauge 
transformation by V{x) G G is mapped as the abelian sub-gauge transformation 
by d^{x) e H in the GMA gauge: (A^^^)^ = (A™)'^''. In particular, for the 

residual gauge transformation by cu{x) = e'^'^^^^'^^^^ E H, we find d'^{x) = uj{x) to 
keep the representative-element condition tr(f2Qj^^(x)0(x)) = imposed above, and 
then A^^^ obeys the ordinary iJ-gauge transformation 

^r^(^) - (^r^(^))^ = ^(^)(^r^ + -o.wix). (2.37) 

For the arbitrary variable ©[A^'^^j = ©[A^^maj defined in the GMA gauge, we find 
d[A^^^]^ = dlA^^^Y"" with d^ eH from Eq.(|;3|), and hence we get a similar cri- 
terion on the SU(A^c) gauge invariance: if 0[A^] is if -invariant as 0[A|^'^'^]'^ = 0[A^] 
for ^u E H, d[A^^^] is also G-invariant, because of 0[/l™]^ = ©[Aj^^]^"" = 

The correspondence between ^gma and ^ma is straightforward. Using flc{x) E 
SU(A^c) satisfying 0(x) = ^(j{x)HQ,c{x), ^gma is expressed as 

fi^^^(x) = fi^(x)fi^^(x). (2.38) 

Then, for regular 0(x), ^c{x) becomes regular, and the singularity of Qua is directly 
mapped to that of f^GMA- However, if singular 0(x) is used, the singularity of f^MA 
can be mapped in </)(x) or Qc{x) instead of I^gma- In this case, the gluon field A^^^ 
is kept to be regular, and the Cartan frame field (f){x) includes the multi-valuedness 
or the singularity, which leads to the monopole. 



Chapter 3 

Origin of Abelian Dominance in 
the MA gauge 



In the abelian gauge, the diagonal and the off-diagonal gluons play different roles 
in terms of the residual abelian gauge symmetry: the diagonal gluon behaves as the 



abelian gauge field, while off-diagonal gluons behave as charged matter fields p3 
Under the U(l)3 gauge transformation by a; = exp{—iip^) G U(l)3, one finds 

1 



A^ ^ iA-r = A'^ + -d,cp (3.1) 

4 - {A-f = A^e^^^ (3.2) 

with A^^ = 75 (^i ± ^^u)- The abelian projection is simply defined as the replace- 
ment of the gluon field A^ = AJ^^ G 5^(2) by the diagonal part Afj, = ^j^-y G 
m(1)3 C sm(2). 

We call "abelian dominance for an operator 0[A^]", when the expectation value 
(O) is almost unchanged by the abelian projection A^ -^ A,^i as (0[yl^]) ^ (0[^^])a.g.) 
when ( )a.g. denotes the expectation value in the abelian gauge. Ordinary abelian 
dominance is observed for the long-distance physics in the MA gauge, and this would 
be physically interpreted as the effective-mass generation of the off-diagonal gluon 
induced by the MA gauge fixing |^, ^ . 

In the lattice formalism, the SU(2) link-variable Uf^{s) can be factorized as 

M,{s) = exp{z{n9l{s) + T29l{s)}) eG/H, 

u^{s) = exp{iTX{s)) eH (3.3) 

with respect to the Cartan decomposition of G = G/H x H into G/H =SU(2)/U(1)3 
and H =U(1)3. Here, the abelian link variable, 

u,{s) = e-^^2 W = ( '^'l^^^ ^_,23 (.) ) G U(l)3 C SU(2), (3.4) 

20 
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plays the similar role as the SU(2)-link variable U^{s) G SU(2) in terms of the resid- 
ual 11(1)3 gauge symmetry in the abelian gauge, and 0^{s) G (— vr, n] corresponds to 
the diagonal component of the gluon in the continuum limit. On the other hand, 
the off-diagonal factor M^(s) gSU(2)/U(1)3 is expressed as 

(3.5) 




with e^{s) = mod|y'(ei)2 + (02)2 ^ ^g, f ] and Xm(s) ^ {-n,n]. Near the contin- 
uum limit, the off-diagonal elements of Mfj^{s) correspond to the off-diagonal gluon 
components. Under the residual U(l)3 gauge transformation by uj(s) = e"*'^^*-' 2 g 
11(1)3, u^{s) and M^(s) are transformed as 

M^) -^ K(^)=^{sW{sW{s + ft) eH (3.6) 

M^{s) ^ M;:{s) = uj{s)M^{s)J{s) eG/H (3.7) 

so as to keep M^(s) belonging to G/H. Accordingly, 0^{s) and c^(s) G C are 
transformed as 

OU^) - el-is) = mod2A0lis) + Wis + ix)-cp{s)}/2] (3.8) 

c,is) -. c-(s)=c,(.)e^^(^). (3.9) 

Thus, on the residual U(l)3 gauge symmetry, u^{s) behaves as the U(l)3 lattice 
gauge field, and 9^{s) behaves as the U(l)3 gauge field in the continuum limit. On 
the other hand, M^(s) and Cfj_{s) behave as the charged matter field in terms of the 
residual U(l)3 gauge symmetry, which is similar to the charged weak boson Wj^ in 
the standard model. 

In this parameterization ( p.3|) , there are two U(l)-structures embedded in SU(2) 
corresponding to e*^*^ and e**''. To clarify this structure, we reparametrize the SU(2) 
link variable as 

^''^'^ ^ (sin^^e^^M cos^^l-^jj' ^^'^^^ 

or equivalently 

[/; = cos 0, cos ^;, Ul = sme,smx„ 

Ul = cos^^sin^j^, f/2 = sin^^cosx^, (3.11) 

with Xfj. = Xfi + (^fi- The range of the angle variable can be redefined as < 6*^ < | 
and —71 < 9f^, Xfi < tt. Here, ([/°, Uj^, U"^, Uf^) forms an element of the 3-dimensional 



22 3 Origin of Abelian Dominance in the MA gauge 



hyper-sphere S^ ^SU(2), because of (f/°)2 + (f/^)^ + (f/2)2 + (f/3)2 _ ^_ p^^ ^ g^^^ 
6^^, both ([/°, U^) and (f/^, f/^) form the two 5*^ ~U(1) subgroups embedded in S^ 
in a symmetric manner. From the parametrization in Eq. (|3.10| ), the SU(2) measure 
can be easily found as 

3 



fdu, ^ jdulululul5{Y.{u;f-i) 

■' ■' a=0 

= ^JJde,sme^,cose,iyx,J\d9l. (3.12) 

In the lattice formalism, the abelian projection is defined by replacing the SU(2) 
link variable U^{s) G SU(2) by the abelian link variable u^{s) G U(l)3. 

3.1 Microscopic Abelian Dominance in the MA 
gauge 

In the MA gauge, the off-diagonal gluon component is strongly suppressed, and 
the SU(2) link variable is expected to be U(l)3-like as U^{s) ~ M/i(s) in the rel- 
evant gauge configuration. In the quantitative argument, this can be expected as 
(C/^(s)m|^(s))ma — 1, where {)ma denotes the expectation value in the MA gauge. In 
order to estimate the difference between U^{s) and u^{s), we introduce the "abelian 



projection rate" -RAbei [0, 0' Hill' which is defined as the overlapping factor as 



i?Abei(s,/i) = -Re tr{U^is)ul{s)} = -Re trM^(s) = cos0^(s). (3.13) 

This definition of -RAbei is inspired by the ordinary "distance" between two matrices 
A,B e GL(iV,C) defined as d^{A,B) = |tr{(A - B)^{A - B)} 0, which leads 
to d^{A,B) = 2 - Re tr(ASt) for A, B gSU(2). The similarity between U^{s) and 
Ufj_{s) can be quantitatively measured in terms of the "distance" between them. For 
instance, if cos 6*^(3) = 1, the SU(2) link variable becomes completely abelian as 



e'^M 




^m(*) = I n .-ie?, 



while, if cos 6*^(5) = 0, it becomes completely off-diagonal as 
-e"**" 



^m(^) ^ ,^x. 



We show in Fig. 3.2 and Fig.p.3| the spatial distribution of the abelian projection 



rate -RAbei = cos 6* as an arrow (sin 6*, cos^). In the MA gauge, most of all SU(2) link 



variables become U(l)3-like. We also show in Fig.|3.4|(a) the probability distribution 



Figure 3.1: Geometrical explanation of the abelian projection rate -RAbei in terms 
oi U = Uq + lUaTa- Tlic deviation from the vertical axis indicates the magnitude of 



the off-diagonal component. This description will be used in Figs 3.2 and 3.3 



P(-RAbei) of the abelian projection rate Rxhc\{.s, fi) in the MA gauge. Here, (-RAbei)/3=o 
in the strong coupling limit (/? = 0) [Q ^ is analytically calculable as 



(^Abel)/3=0 = (cos6'^(s))^=o 



J dU^{s) COS 9 ^{s) 



JdU.is) 

gde^js) sin e^is)cos^e^is) _2 

/o^t/^^(s)sin^^(s)cos^^(s) 3' 



(3.14) 



using Eq. (|3.12|) . In the MA gauge, -RAbei approaches to unity as shown in Fig.|3^(a). 
The off-diagonal component of the SU(2) link variable is forced to be reduced. As 
a typical example, one obtains (-RAbei)MA — 0.926 on 16^ lattice with /3 = 2.4. In 
Fig.p^(b), we show the abelian projection rate (-RAbei)MA as the function oi p. For 



larger (3, (cos6'^(s))ma becomes slightly larger. Without gauge fixing, the average 
(-RAbei) is found to be about | without dependence on /3. In the continuum limit 
in the MA gauge, Uj^{s) and U^{s) become at most 0{a), and therefore (-RAbei)MA 
approaches to unity as (-RAbei)MA = 1 + 0(a^) due to the trivial dominance of U^{s), 
which differs from abelian dominance in the physical sense. The remarkable feature 
of the MA gauge is the high abelian projection rate as (-RAbei) ma — 1 in the whole 
region of p. In fact, we find (-RAbei)MA — 0.88 even for the strong coupling limit 
/? = 0, where the original link variable t/^ is completely random. Thus, abelian 
dominance for the link variable f/^ is observed at any scale in the MA gauge. 

To understand the origin of the high abelian projection rate as (-RAbei)MA — 1, 
we estimate the lower bound of (-RAbei)MA in the MA gauge using the statistical 



□ 



Figure 3.3: Local abelian rate -RAbei = cos^^(s) at /3 =2.4 on 16"^ lattice in the MA 
gauge. The meaning of the arrow is shown in Fig. |3.1| . 
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Figure 3.4: (a) The probability distribution P(i?Abei) of the abehan projection rate 
i?Abei at /5 = 2.4 on the IG"' lattice from 40 gauge configurations. The solid curve 
denotes P(-RAbci) in the MA gauge, and the dashed line denotes P(-RAbci) without 
gauge fixing, (b) The average of the abelian projection rate (i?Abci) in the MA gauge 
as the function of j3. For comparison, we plot also (-RAbci) without gauge fixing. 
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consideration. The MA gauge maximizes 

RMA[U^]=T.^T{U,is)T,Ulis)Ts}=tT{nY.Ms)), (3.15) 

where (f>i_i{s) = f/^(s)r3f/^(s) is an su{2) element satisfying 0^ = 1. Denoting 0^(s) = 
0^(s)r", we parameterize the 3-dimensional unit vectors 0^ = (0^, 0^, 0^) E S^ {fj, = 
1,2,3,4) as 0^ = (sin26'^cosxAt; sin26'^sinx^, cos26'^) using Eqs. (|3.3| ) and ( p.5|) . 



The MA gauge maximizes the third component 0^ using the gauge transformation. 

Under the local gauge transformation by V{s) = 1 + {V{so) — l}Ssso ^ SU(2), 0^(so) 
is transformed as the unitary transformation, 

M^o) ^ 0;(so) = Viso)Mso)V-\so), (3.16) 

which leads to a simple rotation of the unit vectors 0^. In the MA gauge, ^0/^ is 

S,fl 

"polarized" along the positive third direction. On the 4-dimension lattice with A^ 
sites, 4A^ unit vectors 0^(s) are maximally polarized by A^ gauge functions V{s) in 
the MA gauge. Then, (-RAbei)MA is expressed as the maximal "polarization rate" 
of 4A^ unit vectors 0^ by suitable A^ gauge functions V{s). On the average, this 
estimation of (-RAbei)MA is approximately given by the estimation of the maximal 
polarization rate of 4 unit vectors 0^ by a suitable rotation with V G SU(2). The 
lower bound of (-RAbei)MA is obtained from the strong-coupling system with /? = 0, 
where link variables Ufj_{s) are completely random. Accordingly, 0^ can be regarded 
as random unit vectors on S"^. The maximal "polarization" of 4 unit vectors 0^ 

is realized by the rotation which moves =^0^/|^0^| to the unit vector (p^ = 

(0,0, 1) in third direction. Here, cos 2^^ after the rotation is identical to the inner 
product between 0^ and 0, because of ■ 0^ = 0^^ ■ 0^ = (0^)^ = cos 2^^. Then, we 
estimate (i?Abei)MA = (cos 6'^) ma at /? = as 

(cos^Xt ^in [dU,} f^Ecos^] 

^ 1 rf p /l ^ 1 _ - - \ 

= {n~/ de^cosO^sme^, rfx^J -Ecos{- cos ^(0-0^)} . (3.17) 



Using this estimation ( |3.17| ), we obtain (-RAbei)MA — 0.844, which is close to the 



lattice result (-RAbei) — 0.88 in the strong coupling limit (/3 = 0). Such a high 
abelian rate (-RAbei)MA in the MA gauge would provide a microscopic basis of abelian 
dominance for the infrared physics. 
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3.2 Abelian Dominance for Confinement Force in 
the MA Gauge 

In this section, we study the origin of abehan dominance on the string tension as 
the confinement force in a semi-analytical manner, considering the relation with 
microscopic abelian dominance on the link variable |^, ^. 

In the MA gauge, the diagonal element cos^^(s) in M^{s) is maximized by the 
gauge transformation as large as possible. For instance, the abelian projection rate 
is almost unity as -RAbei = (cos6'^(s))ma — 0.93 at /3 = 2.4. Then, the off-diagonal 
element e*^''*^*-' sin^^(s) is forced to take a small value in the MA gauge due to 
the factor sin^^(s), and therefore the approximate treatment on the off-diagonal 
element would be allowed in the MA gauge. Moreover, the angle variable x^(s) 
is not constrained by the MA gauge-fixing condition at all, and tends to take a 
random value besides the residual 11(1)3 gauge degrees of freedom. Hence, Xfj,{s) 
can be regarded as a random angle variable on the treatment of M^(s) in the MA 
gauge in a good approximation. 

Let us consider the Wilson loop {Wc[U^{s)]) = (tr]Jf/^(s)) = {tiY[{M^{s)u^{s)}) 

c c 

in the MA gauge. In calculating {Wc[U^{s)]), the expectation value of e*^^*^*-* in 

M^{s) vanishes as 

(e^XM(^)) ^ / dx,{s) exp{ix,{s)} = 0, (3.18) 

<J u 

when x^l{s) behaves as a random angle variable. Then, within the random- variable 
approximation for Xf^is), the off-diagonal factor M^(s) appearing in {Wc[Uf^{s)]) is 
simply reduced as a c-number factor, M^{s) — > cos6'^(s) 1, and therefore the SU(2) 
link variable U^{s) in the Wilson loop {Wc[U^{s)]) is simplified as a diagonal matrix, 

Uf,{s) = M^{s)u^{s) -^ cos^^(s)u^(s). (3.19) 

Then, for the I x J rectangular C, the Wilson loop ^^^[[/^(s)] in the MA gauge 
is approximated as 

i=l i=l j=l 

L 



(exp{^ln(cos6'^.^(si))})MA (W^cK('5)])ma 



i=\ 



~ exp{L(ln(cose^(s)))MA} (W^cK(s)])ma, (3.20) 

L 

where L = 2(7 + J) denotes the perimeter length and W^c[^a«('^)] — ^^11 '"/^^(■^i) 

L 

the abelian Wilson loop. Here, we have replaced ^ ln{cos(6'^^(sj)} by its average 



i=l 
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L(ln{cos6'^(s)})MA in a statistical sense, and such a statistical treatment becomes 
more accurate for larger /, J and becomes exact for infinite /, J . 
In this way, we derive a simple estimation as 

iy°*^ = {Wc[U^{s)])/{WcK{s)])ma ^ exp{L(ln(cos^^(s)))MA} (3.21) 

for the contribution of the off-diagonal gluon element to the Wilson loop. From 
this analysis, the contribution of off-diagonal gluons to the Wilson loop is expected 
to obey the perimeter law in the MA gauge for large loops, where the statistical 
treatment would be accurate. 

Now, we study the behavior of the off-diagonal contribution 
W^^ = {Wc[U^{s)])/{Wc[Uf,{s)])MA in the MA gauge using the lattice QCD, con- 
sidering the theoretical estimation Eq.( ^.2lD . As shown in Fig. |3.5| , we find that 



W^^ seems to obey the perimeter law for the Wilson loop with /, J > 2 in the 
MA gauge in the lattice QCD simulation with f3 = 2.4 and 16^. We find also that 
the behavior on W^^ as the function of L is well reproduced by the above ana- 
lytical estimation with microscopic information on the diagonal factor cos6'^(s) as 
(ln{cos^^(s)})MA ^ -0.082 for /3 = 2.4. Thus, the off-diagonal contribution Wf to 
the Wilson loop obeys the perimeter law in the MA gauge, and therefore the abelian 
Wilson loop {Wc[u^{s)])ma should obey the area law as well as the SU(2) Wilson 
loop PVcff^Ml-^)]- From Eq. (|3.21| ), the off-diagonal contribution to the string tension 



vanishes as 

Aa = (Tsu(2) - O'Abel (3.22) 

1 1 

= - lim -—-\n{WRy,T[U^{s)]) + \im -;^ ln(IVRxrK(5)])MA 

R,T~*co til R,T—KX} JrCI 

p _|_ T^ 

~ -2(ln{cos0^(s)})MA^lim ^-^=0. 

Thus, abelian dominance for the string tension, o"su(2) = <^Ahci, can be proved in the 
MA gauge by replacing the off-diagonal angle variable Xt^i^) ^^ ^ random variable. 
The analytical relation in Eq.( p.2lj ) indicates also that the finite size effect on R 
and T in the Wilson loop leads to the deviation between the SU(2) string tension 
crsu(2) and the abelian string tension o"Abei as crsu(2) > ^"Abci in the actual lattice QCD 
simulations. Here, we consider this deviation Act = crsu(2) — crAbei in some detail. 
Similar to the SU(2) inter-quark potential Vsu(2)(^) from ( VFsu(2) ) = {'^[U^i{s)]), 
we define the abelian inter-quark potential VAhei{i^) and the off-diagonal contribution 
Vof[{r) of the potential from ( P^Abei ) = ( W^[w^(s)] ) and Wos, respectively, 

Vsvi2){r) = -^ln(iysu(2)(i?xT)), 
VAUr) = -^ln(H^Abei(i?xr)), 

Ta Ta { WAhei{R x T) ) 

= Vsu(2)(r)-VAbei(r), (3.23) 
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Figure 3.5: The off-diagonal gluon contribution on the Wilson loop, W°^ = 
(w^ \u (s)]) ' ^^ ^^^ function of the perimeter length L = 2(J + J) in the MA gauge 
on 16^ lattice with /3 = 2.4. The thick line denotes the theoretical estimation in 
Eq. (|3l2lD with the microscopic input (Injcos 6'/,(s) })ma ^ -0.082 at /3 = 2.4. The 
data of the Wilson loop with / = 1 or J = 1 are distinguished by the thin cross. 
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Figure 3.6: The inter-quark potential Vir) as the function of the inter-quark 
distance r. The lattice data are obtained from the Wilson loop in the MA gauge on 
16^ lattice with /5 = 2.4 and T = 7. The square, the circle and the rhombus denote 
the full SU(2), the abelian and the off-diagonal contribution of the static potential, 
respectively. The thin line denotes the theoretical estimation in Eg. ( p. 241) . Here, 
the lattice spacing a is determined so as to produce a = 0.89 GeV/fm. Due to the 
artificial finite-size effect of the Wilson loop, the off-diagonal contribution V°^ gets 
a slight linear part. 



where r = Ra denotes the inter-quark distance in the physical unit. We show in 



Fig.|3^ Vsu(2)('^), VAbei(^) and Vos{r) extracted from the Wilson loop with T = 7 in 
the lattice QCD simulation with /3 = 2.4 and 16'^. As shown in Fig. p.6| , the lattice 
result for Kff(^) seems to be reproduced by the theoretical estimation obtained from 
Eq-O, 



Kflf(r) = Vsu(2){r) - VAbel(r) 



2{R + T) 
Ta 



(ln(cos6'^(s)))MA 



(3.24) 



using the microscopic information of (ln(cos0^(s)))MA = —0.082 at /5 = 2.4. From 
the slope of Vos{r) in Eq. (|3.24|) , we can estimate Act = crsu{2) — o"Abei in the physical 
unit as 



Aa = a-su{2) - o-Abei ^ -2(ln ( cos6'^(s) ) ) 



MA 



-{\n{l-sm^9,{s))) 



1 



MA 



Ta^' 



(3.25) 



In the MA gauge, sin^ ^'^(■s) takes a small value and can be treated in a perturbation 
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manner so that one finds 

Aa ^ {sm'e,{s))MA^, = {{Uiis))' + {UlMf)MA^,. (3.26) 

Near the continuum hmit a ^ 0, we find f/" ^ aeA't^/2 (a=l,2,3) from f/^ = 
giaeA^T''/2^ and then we derive the relation between Act and the off-diagonal gluon 
in the MA gauge as 

Aa ^ ^{{eAir + (eAj)2)MA = ^(Mi)' + My')MA, (3.27) 

where t = Ta is the temporal length of the Wilson loop in the physical unit. In 
Eq. (|3.27|) , {{eAj^)^ + (eA^)^)MA is the off-diagonal gluon- field fiuctuation, and is 
strongly suppressed in the MA gauge by its definition. It would be interesting to 
note that microscopic abelian dominance or the suppression of off-diagonal gluons 
in the MA gauge is directly connected to reduction of the deviation Act in Eq.( p.27| ). 
Since ((eA^)^ + (ey4^)^)MA is a local continuum quantity, it is to be independent on 
both a and t. Hence, the deviation Act between the SU(2) string tension (Tsu(2) and 
the abelian string tension o"Abei can be removed by taking the large Wilson loop as 
t —>■ oo OT the small mesh as a — * with fixed t. 



3.3 Gluon Field in the MA Gauge with the U(l)3 
Landau Gauge 

In the MA gauge, the linear confinement potential can be almost reproduced only 
by the abelian degrees of freedom, which is called as abelian dominance on the string 
tension. In this section, we study the probability distribution of the gauge field such 
as the abelian gauge field, the abelian field strength, the off-diagonal gluon in the 
MA gauge. 

From the abelian angle variable Of^{s), the abelian field strength 6^u{s) is defined 
as 

e^,{s) = mod2,(9 A e')^,{s) e (-TT, tt], (3.28) 

where {d A 6''^)^j,(s) is expressed as 

{d A 9'),,{s) ^ 9,,is) ^ 9lis) + 9l{s + /i) - 9^3 + 0) - 9lis). (3.29) 

The abelian field strength 9^u{s) is related to the abelian plaquette as 

°^r (^) = n,{s)u,{s + fi)ul{s + u)ul{s) = e^^"-(^). (3.30) 

In terms of the U(l)3 gauge symmetry remaining in the abelian gauge, the abelian 
field strength 9^y{s) is a U(l)3 gauge-invariant quantity because of the gauge invari- 
ance of the abelian plaquette, while both ^^(s) and {d A 9^)^v{s) are U(l)3 gauge 
variant. 



Figure 3.8: The abelian angle variable Of^{s) in the MA gauge with U(l)3 Landau 
gauge fixing. 
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To investigate the features of U(l)3 gauge-variant quantities hke the abehan 
angle variable 9^{s), it is necessary to fix the residual U(l)3 gauge degrees of freedom 
in addition to the abelian gauge fixing. In this paper, we introduce U(l)3 lattice 
Landau gauge |Q, where the gluon field is mostly continuous under the constraint 



of the MA gauge condition, and the lattice field can be compared with continuum 
field variable more directly. 

The U(l)3 lattice Landau gauge is defined by maximizing 

RL[U^] = J2^eu^{s) (3.31) 

S,fJ. 

using the residual U(l)3 gauge transformation. In the U(l)3 Landau gauge, the 
abelian angle variable 9f^{s) is suppressed as small as possible. In the continuum 
limit a ^ 0, the abelian gauge field A^{x) satisfies the ordinary Lorentz-gauge 
condition as in QED, (9^A^ = 0. 

We show in Figs ^T^ and p.8| the configuration of the abelian angle variable Ol^{s) 



before and after the U(l)3 Landau gauge fixing in the MA gauge at /? = 2.4. The 
magnitude of the angle variable is found to become small and continuous in the 
U(l)3 Landau gauge. We show also in Fig.pTP| the probability distribution P{9^) 
of the abelian angle variable Ol^{s) G (— 7r,7r] in the MA gauge with and without 
U(l)3 Landau gauge fixing. We find that the whole shape of the distribution seems 
Gaussian-type peak around 9^ = in the U(l)3 Landau gauge, while 0^{s) is not 
settled without the U(l)3 gauge fixing. 



We show in Fig. p.lO| the probability distribution P{9^y) of the two form 9^ 



flU 



S = 



{d A9^)^^{s) of the abelian field 9^{s). Without the gauge fixing, there appear three 
peaks around 9^^{s) = — 27r, 0, 27r. In the MA gauge, because of U^ ~ u^, not 

only the SU(2) action S = S//3 = ^^ [^ trn^j,(s)] but also the abelian action 

^ bei ^ ^ [1 — -trn^j'j"^'(s)] is suppressed by the action factor e~^^ in the partition 

S,fJ,>U 

functional. Since the abelian action is written by ^ [1 — cos(6'^i,)], P{9^^) has peaks 

S,fl>U 



around cos{9^,y) = 1. As shown in Fig. |3.10| , most of ^^i,(s) distribute around 9^y{s) 



in the U(l)3 Landau gauge, because the abehan gauge field is mostly continuous 
in the U(l)3 Landau gauge. On the other hand, the probability distribution P{9^^) 
of the abelian field strength i9^i, G (— vr, tt] is U(l)3 gauge invariant. The whole shape 
of P{9fj,i^) is Gaussian-type, as shown in Fig. p.llj . 

Finally in this section, we investigate the off-diagonal phase variable x^l{^) ^^ 
Eq. (|3.5|) . We show in Fig. |3.12| (a) the probability distributions PiXn) ^'^'^ P{9f^) at 



P = 2.4 in the MA gauge with U(l)3 Landau gauge. Unlike P{9^^), P{Xfj.) is flat 
distribution without any structure. This property on the off-diagonal element would 
lead to the validity of the random- variable approximation for x^l{s)> which has been 
used for estimation of the Wilson loop in Eq. (|3.2CI|) . 
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Figure 3.9: The probabihty distribution P{0^) of the abehan angle variable 9^{s) G 
(— TT, tt] in the MA gauge on 16^ lattice with (3 = 2.4. The solid curve denotes P{0^) 
in the U(l)3 Landau gauge fixing after the MA gauge fixing, and the dashed line 
denotes P{0^) without the U(l)3 gauge fixing. Without the U(l)3 gauge fixing, the 
angle variable 6f^{s) is randomly distributed, while Of^{s) has a Gaussian-type peak 
around Ol^{s) = in the U(l)3 Landau gauge. 
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Figure 3.10: The probability distributions PiOfj^^) of the two form ^^j, = {d /\6) ^j^^^s) 
in the MA gauge with and without U(l)3-Landau gauge fixing, which are denoted 
by the solid and dashed curves, respectively. In the U(l)3-Landau gauge, P{9^^) 
has single peak around 6^i, = 0. 
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Figure 3.11: The probabihty distribution P{Ofj_y) 
which is 11(1)3 gauge invariant. 



of the abehan field strength 6^y{s) 



3.4 Randomness of Off-diagonal Gluon Phase and 
Abelian Dominance 



In this section, we reconsider the origin of abelian dominance in the MA gauge in 
terms of the properties of the off-diagonal element 



e^^'*(')sin^^(s) 



(3.32) 



in M^{s) in the link variable U^{s), considering the validity of the random- variable 



approximation for Xf^i^) ^^ the MA gauge with U(l)3 Landau gauge ^7^. In this 



treatment, the contribution of the off-diagonal element in the link variable f//i(s) is 
completely dropped off, and its effect indirectly remains as the appearance of the 
c-number factor cos 6*^(5) in the link variable. Such a reduction of the contribution 
of the off-diagonal elements is brought by the two relevant features on the two local 
variables, 0^{s) and Xfii^)^ ^^ the MA gauge. One is microscopic abelian dominance 
as (cos6'^(s))ma — 1 in the MA gauge, and the other is the randomness of the 
off-diagonal phase variable Xfi{s)- 

1. In the MA gauge, microscopic abehan dominance holds as (cos6'^(s))ma — 
1, and the absolute value of the off-diagonal element |c^(s)| = |sin6'/i(s)| is 
strongly reduced. Such a tendency becomes more significant as (3 increases. 

2. The off-diagonal angle variable Xfii^) is not constrained by the MA gauge-fixing 
condition at all, and tends to be a random variable. In fact, Xfii^) is affected 
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Figure 3.12: The probabihty distributions P{x^i) (sohd hne) and P{du) (dash- 
dotted curve) in the MA gauge with the U(l)3 Landau gauge at /3 = 2.4 on the W^ 
lattice from 40 gauge configurations. 



only by the QCD action factor e~^'^<5CD j^ the QCD generating functional, but 
the effect of the action to x^l{s) is quite weaken due to the small factor sin 6*^(5) 
in the MA gauge. The randomness of Xni^) tends to vanish the contribution 
of the off-diagonal elements. 

Here, the randomness of the off-diagonal angle- variable Xfi{^) is closely related to 
microscopic abelian dominance. In fact, the randomness of Xfi{^) is controlled only 
by the action factor 6"^*^*^^° in the QCD generating functional, however the effect of 
the action to Xij.{^) is quite weaken due to the small factor sin 0^{s) in the MA gauge. 



in the link variable U^{s) 



Near the 



because x^l{^) always accompanies sin 6'^(s) 
strong-coupling limit /5 ~ 0, the action factor e~^'^*3CD brings almost no constraint 
on Xfiis) in the MA gauge. The independence of Xfii^) from the action factor is 
enhanced by the small factor sin 6*^(5) accompanying Xi^i^)- Hence, Xi^i^) behaves 
as a random angle-variable almost exactly, and the contribution of the off-diagonal 
element is expected to disappear in the strong-coupling region. As /? increases, the 
action factor e"^'^'^'^^ becomes relevant and will reduce the randomness of Xt^i^) 
to some extent. Near the continuum limit f3 -^ 00, however, the factor sin^^(s) 
tends to approach in the MA gauge as shown in Fig.^.4|(b), and hence such a 
constraint on Xij.{^) from the action is largely reduced, and the strong randomness 
of Xfi{^) is expected to hold there. Moreover, the reduction of the absolute value 
|c^(s) I = I sin 9fj^{s) I itself further reduces the contribution of the off-diagonal element 
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Figure 3.13: The probabihty distribution P{Ax) of the correlation A^ = 
inodTr{\Xfi{s) — x^J.{s + i>)\) in the same gauge at /3 = (thin line), 1.0 (dotted 
curve), 2.4 (solid curve), 3.0 (dashed curve). 



\Ct.[S 



in the MA gauge. 

Now, we examine the randomness of Xfi{^) using the lattice QCD simulation. 
We calculate the correlation between Xfii^) and Xfii^ + ^) iu the MA gauge with 
the U(l)3-Landau gauge. If x^J.{s) is an exact random angle variable, no correlation 
is observed between Xij.{^) ^-ud XtJ.{^ + ^)- We show in Fig. |3.13| the probability 
distribution P{Ax) of the correlation 



Ax(s) = d{Xf.{s),Xi,is + i))) =mod^|x/.(s) - Xf^i-s + i))\ G [0,7r], 



(3.33) 



which is the difference between two neighboring angle variables, at (3=0, 1.0, 2.4, 3.0. 
In the strong-coupling limit /? = 0, Xm("5) is a completely random variable, and there 
is no correlation between neighboring Xfi- In the strong-coupling region as /? < 1.0, 
almost no correlation is observed between neighboring X/i, which suggests the strong 
randomness of Xfi{^)- As a remarkable feature, the correlation between neighboring 
Xfj, seems weak even in the weak-coupling region as j3 > 2.4, where the action 
factor e~^^^'^° becomes dominant and remaining variables ^j^(s) and 0^{s) behave 
as continuous variables with small difference between their neighbors as A6f^ ~ and 
A6^ ~ 0. Such a weak correlation of neighboring x^l would be originated from the 
reduction of the accompanying factor sin^^(s) in the MA gauge. Moreover, in the 
weak-coupling region, the smallness of sin^^(s) makes c^(s) more irrelevant in the 
MA gauge, which permits some approximation on x^l{^)■ Thus, the random- variable 
approximation for x^l{^) would provide a good approximation in the whole region of 
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(3 in the MA gauge. To conclude, the origin of abehan dominance for confinement 
in the MA gauge is stemming from the strong randomness of the off-diagonal angle 
variable x^l{.s) and the strong reduction of the off-diagonal amplitude | sin6'^(s)| as 
the result of the MA gauge fixing. 

3.5 Comparison with SU(2) Landau Gauge 

In this section, we study the feature of the MA gauge in terms of the concentration 
of the gluon field fluctuation into the U(l)3 sector by comparison with the SU(2) 



Landau gauge |Q. In the MA gauge, off-diagonal gluon components are forced to be 
small by the gauge transformation. Instead, the gluon field fluctuation is maximally 
concentrated into the abelian sector, and monopoles appear in the abelian sector as 
the result of the large fluctuation of the abelian field component. For the qualitative 
argument on the share of the gluon fluctuation into each component, we measure 
{U;{s)) (a =1,2,3) and 

In the MA gauge with the U(l)3 Landau gauge, we find a strong concentration 
of the gluon field fluctuation into the abelian sector as {U^{s)) = {W^X^)) -C {U^{s)) 
and -Rdiag ^ |; as a typical example, we find at /3 = 2.4 (f/^(s)) = {U'^{s)) ^ 0.067, 
{U^{s)) ~ 0.43 and i?diag ^ 0.68. 

For comparison, we consider the SU(2) lattice Landau gauge [Q defined by 
maximizing 

Rl[U,] = j:^TU,is) = 2Y^U'^{s), (3.35) 

where all the lattice gluon components fields become mostly continuous owing to 
the suppression of their fluctuation around Ufj_{s) = 1. In the continuum limit, 
this gauge fixing condition coincides the ordinary SU(2) Landau gauge condition 
dfj,Af^ = 0. In the lattice SU(2) Landau gauge, one finds -Rdiag = | and (f^^(s)) = 
(f/2(s)) = (f/3(s)) < 1, for instance (t/;j(s)) = 0.076 at /3 = 2.4, so that all the 
gluon components are forced to be small equally. 

In the SU(2) Landau gauge, the local symmetry of SU(2) is fixed, and only 
the global SU(2) symmetry remains, because Rl[U^] = ^trf/^(s) is invariant by 

S,fJ. 

any global gauge transformation. In order to compare with the MA gauge, we fix 
SU(2)/(U(1)3 X Weyl) in this global SU(2) symmetry by the additional condition 
so as to maximize RMA[Ufj] = '^^^iU^{s)T3UJ^{s)Ts) by the remaining global SU(2) 

S,fJ, 

gauge transformation. Here, the SU(2) Landau gauge with SU(2)giobai/(U(l)3giobaiX 
Weyl) fixing is regarded as a kind of the abelian gauge. Then, we can extract 
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Figure 3.14: Comparison between the MA gauge and the SU(2) Landau gauge 
with global SU(2)/(U(1) x Weyl) fixed, (a) The probability distribution P(^3) of 
the lattice angle variable 9^{s) in the SU(2) Landau gauge (solid curve) and the 
MA gauge (dotted curve ) with f] = 2.4 on 16^ lattice. The large fiuctuation on 
the link variable Of^{s) disappears in the SU(2) Landau gauge. The abelian part of 
the interquark potential V{r) as the function of interquark distance r in the SU(2) 
Landau gauge and in the MA gauge on 16^ lattice with j3 = 2.4. For comparison, 
the full SU(2) result is added. 
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the abehan variable 0^{s) even for this SU(2) Landau gauge. Figure p.l4| (a) shows 
the probabihty distribution P{9^) of the abehan angle variable 0^{s) in the SU(2) 
Landau gauge and in the MA gauge at /3 = 2.4. We show also in Fig. |3.14| (b) the 



interquark potential in the abelian sector evaluated from the abelian Wilson loop 
in these gauges. Although the global shape of the distribution P{0^) in the SU(2) 
Landau gauge is similar to that in the MA gauge except for the reduction of the 
large fluctuation apart from 6f^{s) = 0, the abelian string tension in the SU(2) 
Landau gauge is much smaller than that in the MA gauge [^. Therefore, the large 
fluctuation ingredient is expected to be responsible for the conflnement property. 

To summarize, in the MA gauge, the fleld fluctuation is maximally concentrated 
into the abelian sector, and hence large fluctuation ingredient appears and the con- 
flnement property is almost reproduced only by the abelian variable. Another clear 
difference between the MA gauge and the SU(2) Landau gauge observed on lattice is 
the density of monopoles appearing in the abelian sector. Indeed, the SU(2) Landau 
gauge includes scarcely monopoles in the abelian sector in comparison with the MA 
gauge, for instance, the ratio on the monopole density is less than 1/10 at /? = 2.4. 
This result seems natural because the SU(2) Landau gauge flxing provides a mostly 
continuous gluon fleld, while the monopole arises from a singular-like large fluctua- 
tion of the abelian fleld as will be shown in Chapter ^ and |^ in detail. In the next 
chapter, we study features of monopole appearing in the MA gauge in relation with 
conflnement and large fleld fluctuation concentrated into the abelian sector. 



Chapter 4 

QCD-Monopole in the Abelian 
Gauge 



In the abelian gauge, QCD is reduced into an abelian gauge theory with QCD- 
monopoles, which appear from the hedgehog-like configuration corresponding to the 
nontrivial homotopy group on the nonabelian gauge manifold, n2(SU(A^c)/U(l)^'=~^) 
_ 2^':"^. The relevant role of the QCD-monopole to the infrared phenomena has 
been studied by using the infrared effective theory and the lattice gauge simula- 
tion |31, 3^. In the dual Ginzburg-Landau(DGL) theory, the linear static quark 



potential, which characterizes quark confinement, is obtained in the monopole con- 



densed vacuum [|33|. In addition, chiral symmetry breaking is also brought from the 
monopole contribution in the DGL theory |^3|, Q. The recent lattice QCD studies 
in the MA gauge suggest monopole condensation in the confinement phase in the 
MA gauge, and show abelian dominance and monopole dominance for nonpertur- 
bative QCD |31|, |3^. Here, monopole dominance means that QCD phenomena are 
described only by the monopole part of the abelian variables in the abelian gauge. 
In this chapter, we study appearance of QCD-monopoles in the abelian sector of 
QCD and clarify the difference between the ordinary QED and abelian projected 
QCD (see Fig. 1.1). 

4.1 Appearance of Monopoles in the SU(2) Sin- 
gular Gauge Transformation 

The abelian gauge fixing, which reduces QCD into an abelian gauge theory, is re- 
alized by the diagonalization of a suitable variable $[A^(x)]. In the continuum 
theory of QCD, the continuous field A^{x) can be taken to be regular everywhere 
in a suitable gauge as the Landau gauge, and then ^[^^(a;)] is expected to be a 
regular function almost everywhere. In the abelian gauge, however, there appears 
the singular point, where the gauge function to diagonalize $[A^(a;)] is not uniquely 

41 
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determined even for the off-diagonal part, and such a singular point leads to the 
appearance of the monopole. 

Here, let us consider the appearance of QCD-monopoles in the abelian gauge 
in terms of the singularity in the gauge transformation [^. For the variable $(x) 
obeying the adjoint transformation, the monopole appears at the "degeneracy point" 
of the diagonal elements of if ■A(a;) = diag{X^{x), \'^{x),- ■ ■ , X^''{x)) after the abehan 
gauge fixing: (i, j)-monopole appears at the point satisfying A*(x) = X^{x). For the 
(z, j)-monopole, the SU(2) subspace relating to i and j is enough to consider, so 
that the essential feature of the monopole can be understood in the SU(2) case 
without loss of generality. Then, we consider the SU(2) case for simplicity. For the 
SU(2) case, the diagonalized element of $(a;) are given by A = ±($? + $i + $1)^''^, 
and hence the "degeneracy point" satisfies the condition $(x) = 0, which is SU(2) 
gauge invariant. This gauge-invariant condition $(x) = can be regarded as the 
singularity condition on $(a:) = $(x)/|$(x)| with |$(x)| = ($"(x)$"(x))^/^. In 
fact, the "degeneracy point" in the abelian gauge appears as the singular point of 
$(x) like the center of the hedgehog configuration as shown in Fig. |4.1| (b) before the 
abelian gauge fixing. 

Since the singular point on ^{x) is to satisfy three conditions $^(a;) = $^(a;) = 
^^{x) = simultaneously, the set of the singular point forms the point-like manifold 
in R'^ or the line-like manifold in R^. We investigate the topological nature near 
the singular point (xo,t) of $(x) for fixed t, i.e., $(xo,t) = |3^. Using the Taylor 
expansion, one finds 

$(x,t) = $'^(x,t)y ^ r"C"^(x-Xo)^ (4.1) 

with C"'^ = ^d^^°'{'Xo, t). In the general case, one can expect detC 7^ 0, i.e., detC > 
or detC < 0, and the fiber-bandle $"(x) can be deformed into the (anti-)hedgehog 
configuration $(x) ~ ±r"x" around the singular point xq by using the continuous 
modification on the spatial coordinate x"^ — i> x" = sgn(detC) ■ C"^(x — xq)^. The 
linear transformation matrix C can be written by a combination of the rotation R 
and the dilatation of each axis A = diag{X^, A^, A^) with A* > as C = sgn(detC)-RA. 
Here, topological nature is never changed by such a continuous deformation. For 
detC > 0, the configuration $(x) can be continuously deformed into the hedgehog 
configuration around xq, $(x) ~ r°x", while, for detC < 0, $(x) can be continuously 
deformed into the anti-hedgehog configuration, $(x) — — r"x'*. Since detC = 
is the exceptionally special case and detC < is similar to detC > 0, we have 
only to consider the hedgehog configuration. This hedgehog configuration around 
the singular point of $(x) corresponds to the simplest nontrivial topology of the 
nontrivial homotopy group H2(SU(2)/U(1)3) = Z^o, and the abelian gauge field has 
the singularity as the monopole appearing from the hedgehog configuration. 

Using the polar coordinate (r, 9, ip) of x, the hedgehog configuration is expressed 



Figure 4.1: Topological structure of variable <l>[A^(a;)] in the abelian gauge fixing in 
the SU(2) QCD. In the abelian gauge, the monopole appears at the singular point 
of ^{x) = $/|$| with |$| = (<|)"$")i/2. (a) For the regular (trivial) configuration 
of $[A^(a;)], no monopole appears in the abelian gauge, (b) For the hedgehog 
configuration of <l>[A^(a;)], the unit-charge monopole appears in the abelian gauge. 



as 

$ = T°'x."' = r sin 6 cos (/? ■ ri + r sin 6 sin ip ■ T2 + r cos 6 ■ t^ 

cos6 e^^'^sin^ \ , . 

e^^sin^ -cos^ j ' ^ ' ' 

and $ can be diagonalized by the gauge transformation with 

H /e^'^cosf sinf \ , . 

^ =Usin| e--cos|j' ^^'^^ 

where 9, ip denote the polar and the azimuthal angles, respectively. Here, on the 
z-axis {9 = or 9 = it), ip is the "fake parameter", and the unique description does 
not allow the (^-dependence on the z-axis. However, at the positive region of z-axis, 
^ = 0, Q^ depends on ip and is multi- valued as 

""-(T A)' (-' 

Such a multi-valuedness of Q^ leads to the divergence in the derivative d^Q^ at 

^ = 0. In fact, d^il^ includes the singular part as cos | {V(p)^ = ^^^j-^ = — —3;, 

which diverges at ^ = 0. By the gauge transformation with Q^ , the variable $ 
becomes $^ = i7$i7''" = rr^, and the gauge field is transformed as 

A^^A'i = niA^ - '-d^)n^. (4.5) 
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For regular A^, the first term QA^j^Q'' is regular, while A?,™^ = —^Qdf^Q^ is singular 
and the monopole appears in the abelian sector originating from the singularity 
of A^™^ [0. To examine the appearance of the monopole at the origin x = 0, 
we consider the magnetic flux $^"^(^) which penetrates the area inside the closed 
contour c{r, 9) = {{r, 9 , (p)\0 < ip < 2ti}. One finds that 



Jc e Jc 



= _lf-,^n^n. = -l^.l±|2i^|, (4.6) 

e Jo dip e 2 2 

which denotes the magnetic flux of the monopole with the unit-magnetic charge 
g = — with the Dirac string [^]. Here, the direction of the Dirac string from the 
monopole can be arbitrary changed by the singular U3(l) gauge transformation, 
which can move e**^ in Q^ from the rs-sector to the off-diagonal sector. In fact, the 
multi-valuedness of ^2 is not necessary to be fixed in r^-direction. Nevertheless, the 
singularity in fldfj^fl^ appears only in the ra-sector, and rs-direction becomes special 
in the abelian gauge fixing. 

The anti-hedgehog configuration of $(x) = — r^'x" provides a monopole with the 
opposite magnetic charge, because anti-hedgehog configuration is transformed to 
the hedgehog configuration by the Weyl transformation. Thus, the only unit-charge 
magnetic monopole appears in the general case of detC 7^ 0. In principle, the 
multi-charge monopole can also appear when detC = 0, however, the condition is 
scarcely satisfied in general, because this exceptional case is realized only when four 
conditions $^ = $2 = $^ = detC = are simultaneously satisfied. To summarize, 
in the abelian gauge, the unit-charge magnetic monopoles appear from the singular 
points of $(a;), however, multi-charge monopoles do not appear in general cases. 

In this way, by the singular SU(2) gauge transformation, there appears the 
monopole with the Dirac string. Here, we consider the role of the off-diagonal 
component in the SU(2) gauge function Q^ to appearance of the monopole, by 
comparing with the U(l)3 gauge transformation. Let us consider the singular gauge 
transformation ii^^^^ = e^'^'^^ E U(l)3 instead of Q^ . This U(l)3 gauge function 
^u(i) jg niulti- valued on the whole region of the z axis (^ = and ^ = vr), and 
^sing ^ _i^u(i)^^^u(i)t g^jgQ j-^g^g g^ siugularlty. The magnetic flux which penetrates 
the area inside the closed contour c{r, 9) = {r, 9 , ip\0 < ip < 27r} is found to be 

$fl"-(0) = /rfx • A'^^^s = - — ^, (4.7) 

Jc e 2 

which corresponds to the endless Dirac string along the z-axis. It is noted that the 
singular U(l)3 gauge transformation can provide the endless Dirac string, however, 
it never creates the monopole. 

The monopole is created not by above singular U(l)3 gauge transformation but 
by a singular SU(2) gauge transformation. Since the multi-valuedness of fl^ is 
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originated from the yj-dependence at ^ = or ^ = tt, we separate the SU(2) gauge 
function (Of) as 



n 

VL = cos-e*'^^^ + ((/^-independent term). 

At ^ = or the positive side of z axis, f2^ coincides with Vl^^^'^ = e^'^'^^ and is 
multi- valued like f2^*^^). Therefore the Dirac string is created at = by the gauge 
transformation Q^ . On the other hand, at ^ = vr or the negative side of 2;-axis, 
(y9-dependent part of Vl vanishes due to cos | = 0, so that the Dirac string never 
appears in Qdfj,fl^ at 6* = vr. Thus, by the SU(2) singular gauge transformation fl^ , 
the Dirac string is generated only on the positive side of the 2;-axis and terminates 
at the origin r = 0, and hence the monopole appears at the end of the Dirac string. 
Around the origin x = 0, the factor cos | varies from unity to zero continuously with 
the polar angle 6, and this makes the Dirac string terminated. Such a variation of 
the norm of the diagonal component cos le*"^ cannot be realized in the U(l)3 gauge 
transformation with r^^*^^-*. In the SU(2) gauge transformation with fi^, the norm 
of the diagonal component can be changed owing to existence of the off-diagonal 
component of fl^ , and the difference of the multi- valuedness between ^ = and 
6 = n leads to the terminated Dirac string and the monopole. In this way, to create 
the monopole in QCD, full SU(2) components of the (singular) gauge transformation 
is necessary, and therefore one can expect a close relation between monopoles and 
the off-diagonal component of the gluon field. 

4.2 Appearance of Monopoles in the Connection 
Formalism 

In this section, we study the appearance of monopoles in the abelian sector of 
QCD in the abelian gauge in detail using the gauge connection formalism. In the 
abelian gauge, the monopole or the Dirac string appears as the result of the SU(A'"c) 
singular gauge transformation from a regular (continuous) gauge configuration. For 
the careful description of the singular gauge transformation, we formulate the gauge 
theory in terms of of the gauge connection, described by the covariant-derivative 
operator Z)^ and D^ = d^ + ieA^j^^x), where 9^ is the derivative operator satisfying 
[d,J{x)]=d,f{x). 

To begin with, let us consider the system holding the local difference of the 
internal-space coordinate frame. We attention the neighbor of the real space-time x^, 
and denote by \q{x)) the basis of the internal-coordinate frame. At the neighboring 
point x^+e^, we express the difference of the internal-coordinate frame as \q{x+e)) = 
Ri;{x)\q{x)) with Re{x) = e^^"^^^ G G being the "rotational matrix" of the internal 
space. We require the "local superposition" on r^ as r-g^+^j = r^^ + r^^ up to 0(e), 
and then we can express re{x) = —ee^A^{x) using a e-independent local variable 
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^4^(2;) G g : \q{x + e)) = e~*^^'''^^*^^^|g(a;)). Then, the "observed difference" of the 
internal space coordinate depends on the real space-time x^, the observed difference 
of the local operator 0{x) between neighboring points, x^ and x^ + e^, is given by 

{q{x + e)\0{x + €)\q{x + £)) - {q{x)\0{x)\q{x)) 

~ e^{q{x)\{df'0{x) + ie[Af'{x),0{x)]}\q{x)) 

= e,{qix)\{[d^ + teA'^ix),Oix)]}\qix))^EM^)\[D'',0{x)]\q{x)). (4.8) 

Here, one finds natural appearance of the covariant derivative operator, D^ = 
d^ + ieA^{x). The gauge transformation is simply defined by the arbitrary internal- 
space rotation as \q{x)) —>■ Q{x)\q{x)) with Q{x) G G, and therefore the covariant 
derivative operator is transformed as D^ -^ D^ = Q{x)D^Q^{x) with Q{x) G G, 
which is consistent with A^ ^ A^ = fl{A^ — ^dfj)Q^ . 

In the general system including singularities such as the Dirac string, the gauge 
field and the field strength are defined as the difference between the gauge connection 
and the derivative connection, 

A, ^ Ud.-B,) (4.9) 

G^, = -{[D^,D,]-[d„d,]). (4.10) 

This expression of G^u is returned to the standard definition Gf^i, = j^[Dfj_,I)u] = 
dfj_Ai^ — di^A^ + ie[A^, Ay] in the regular system. By the general gauge transformation 
with the gauge function (], the field strength G^y is transformed as 

= 9^A^ - 9,< + ie[A% A[?] + -{Sl[d,, dy\n^ - [d,, dy]) 



^- - e 



= {d,A^ - dyA^^) + ze[A^^, A^] + -Q[d„ dyp 

= G^^ + G^-^^ + G^^. (4.11) 

The last term remains only for the singular gauge transformation on ^2^ and Q^^^^ , 
and can provide the Dirac string. 

Figure |3 shows the SU(2) field strength G^^^^\ Gj^jJ''^"'''' and G^^ in the abelian 
gauge provided by f2^ in Eq. (^l3|) . The linear term G^^y^^^ = {dfj,A^ — dyA^) includes 
in the abelian sector the singular gauge configuration of the monopole with the Dirac 
string, which supplies the magnetic flux from infinity. Since each component satisfies 
the Bianchi identity 9"*Gj^^*^'^'' = d"*{d A A^)afj, = 0, the abelian magnetic flux is 
conserved. The abelian part of Gj^^''^^^'" = ie[A^,A^], (Gjj!''^'^^'')^ = -eiA^Al - 
AlA'l), includes the effect of off-diagonal components, and it is dropped by the 
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Abelian Projected QCD 

Figure 4.2: Appearance of monopoles in abelian projected QCD(AP-QCD). After 
the abelian gauge fixing, monopole with the Dirac string appears from G^'^J^'^'^ in 



Eq. (|4.11| ) and the "anti-Dirac string" appears in the singular part G^^^. The off- 



diagonal contribution GJJ^^™*''^'' = ie[A^, A^] forms the anti- monopole configuration 
and compensates to the singularity of the other parts. As the result, the monopole 
without the Dirac string appears in the abelian field strength J^^^ in AP-QCD. 



abelian projection. The last term G^™^ = ^^[d^, dy\Vl) appears from the singularity 
of the gauge function f2, and it plays the important role of the appearance of the 
magnetic monopole in the abelian sector. 

First, we consider the singular part G^™^. In general, G^"^ disappears in the 



regular point in Vt. It is to be noted that G^™^ is found to be diagonal from the 
direct calculation with ^2^ in Eq. 



O C- ^ 

= - {9^^l9u2 - 9fM29ui ) ^ [c^i ,02\'^- n 

= — {9^ll9u2 - 9^i29ui)0{x3)6{xi)6{x2) ■ — , 



(4.12) 



where we have used relations, 
[di,d2](p = -2Tr6{xi)6{x2) 



1 + cos6' 



5{xi)5{x2) = 9{xs)5{x,)5ix2). (4.13) 



The off-diagonal component of Q^[d^,du]Q^'' disappears, since the singularity ap- 
pears only from (/^-dependent term. As a remarkable fact, the last expression in 
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Eg . (14.121) shows the terminated Dirac string, which is placed along the positive z- 
axis with the end at the origin. Hence, in the abelian part of the SU(2) field strength, 
G^'^s leads to the breaking of the U(l)3 Bianchi identity. 



2 



K = d''*Gt,^ = ;j^a/^9"G^7 = —e^^ud''{6{x,)6ix2)eix3)}^ 



= —9,oS{xi)6{x2)6{x3)^, (4.14) 

which is the expression for the static monopole with the magnetic charge g = — 
at the origin. Thus, the magnetic current /c^ is induced in the abelian sector by 
the singular gauge transformation with Q^ and the Dirac condition eg = 4n is 
automatically derived in this gauge-connection formalism. 

In the covariant manner, G^™^ is expressed as G^™^ = ^*(^ A k)^i, using the 



f^i^ n-d 

monopole current k^ in Eg. ([4. 14 ) and a constant 4- vector n^. Actually, for the above 
case, one finds for ra^ = (7^3 

—^*{nAk)f,i, = j dx'^{x'i\^--\x.^)e^ymn^ — 5{xi)5{x2)5{x'^)^ 

4vr T3 

= — {9^ll9u2 - gt,2gui)0{xz)K^i)Kx2) — 

= ^fi^[9^,9,]fi^t = G5g^ (4.15) 

using the relation (a;„|;;^|x„) = 9{xn — x^). 

Thus, the last term Gf^™^ corresponds to the Dirac string terminated at the origin. 
Since Gj^^*^'^'' shows the configuration of the monopole together with the Dirac string, 
the sum of Gj^^*^'^'' + G^™^ provides the gauge configuration of the monopole without 
the Dirac string in the abelian sector. Thus, by dropping the off-diagonal gluon 
element, Gjj'j™''^'' vanishes and the remaining part (G!^°'^'*'' + Gf,"s)^ describing abelian 
projected QCD includes the field strength of monopoles. 

Next, we consider the role of off- diagonal gluon components for appearance of the 
monopole. The gluon field is divided into the regular part QA^Q'' and the singular 
part —-fld^Q'^. Since we are interested in the behavior of the singularity, we neglect 
the regular part of the gluon field. Then, Gjj^^™'"^'' is written as 

= -{d,A^-d,A]l)-'-n[d,,d,p, (4.16) 
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where the last term appears as the breaking of the Maurer-Cartan equation. In 
the abelian gauge, the singularity of the monopole appearing in G^™^^^ + G^'^s jg 
exactly canceled by that of G|^'^™°^''. Thus, in the abelian gauge, the off-diagonal 
gluon combination {G)i^'''''''f = -e{ {A^Y {A^f - {A^)\A^f } includes the field 
strength of the anti-monopole, and hence the off-diagonal gluons {A^Y and {A^Y 
have to include some singular structure around the monopole. 

The abelian projection is defined by dropping the off-diagonal component of the 
gluon field A^, 

Accordingly, the SU(2) field strength G^^ is projected to the abelian field strength 
T = F — 

G"",. = {G'iy-^ = {d,A^,-d^A^^)+ie[A^^,A^] + -n[d,,d,]n^ 

i 
_i 

e 

9,,A - 9,^^ - -F^r, (4.18) 



^^,u = df,A^ - djyAf, + -fi[<9^, di,]n^ 



where J^^™^ = -F^"^-f- = —lfl[dfj_,d,y]Q'' is diagonal and remains. Here, the bilinear 
term 2e[A|^,y4^] vanishes in AP-QCD because it is projected to ie[A^,Au] = by 
the abelian projection. The appearance of J^f^^^^ leads to the breaking of the abelian 
Bianchi identity in the U(l)3 sector, 

a"*^.^ = -9"* J-:^ = d''*{'-n%, d,,]n^} = k^, (4.19) 

where Eg. ( [4. 151) is used. Thus, the magnetic current fc^ is induced into the abelian 
gauge theory through the singularity of the SU(2) gauge transformation. 

Here, we compare AP-QCD and QCD in terms of the field strength. The SU(A^c) 
field strength G^u is controlled by the QCD action, Sqcd = J d'^x{ — ^tTGf^^G^'^}, so 
that each component G^^ cannot diverge. On the other hand, the field strength jF^j, 
in AP-QCD is not directly controlled by Sqcd, since the QCD action includes also 
off-diagonal components. It should be noted that the point-like monopole appearing 
in AP-QCD makes the U(l)3 action S'Abei = ! d'^x{ — \iTTfj,uT^^} divergent around 
the monopole, such a divergence in JF should cancel exactly with the remaining 
off-diagonal contribution from G^^"^^^^^ to keep the total QCD action finite. Thus, 
the appearance of monopoles in AP-QCD is supported by the singular contribution 
of off-diagonal gluons. In this way, abelian projected QCD includes monopoles 
generally. 
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4.3 Monopole Current in the Lattice Formalism 

In this section, we show the extraction of the monopole current in the lattice for- 
malism [Q. The monopole in lattice QCD is defined in the same manner as in the 



continuum theory. 

The abehan field strength 9^u{s) is defined as Ofj_^{s) = mod27r(9 A 9^)^iy{s) G 
(— 7r,7r], which is U(l)3 gauge invariant. In general, the two form of the abelian 
angle variable 9^{s) is divided as 

V(s) = (9 A ^')^.(s) = V(^) + 27rn,..(s), (4-20) 

where n^j^^s) G Z corresponds to the quantized magnetic flux of the "Dirac string" 
penetrating through the plaquette. Although n^^ ^ provides the infinite magnetic 
field is the continuum limit as 27m^i,/a, the term 2Tin^y{s) does not contribute 
to the abehan plaquette n^^^^(s), and it is changed by the singular U(l)3 gauge- 
transformation as 9^{s) -^ 9f^{s) + df^ip{s) with (p{s) being the azimuthal angle. 
Thus, 27m^i, corresponds to the Dirac string as an unphysical object. 
The monopole fc'"*(s) is defined on the dual link as ||5^ , 



1 



kl^\s) = —d^*9^^{s) = -d^*n^^{s), (4.21) 

using the abelian field strength 9^^{s). Here, /i;^"*(s) is defined such that the topolog- 
ical quantization is manifest, fci"*(s) G Z. In this definition, for instance, one finds 
^cT* — 2^ijkdinjk and k'f* = {i = 1,2,3) for the static monopole. The magnetic 
charge of the monopole on the dual lattice is determined by the total magnetic flux 
of the Dirac strings entering the cube around the monopole (see Fig. [4.3| (a).) 

We show in Fig. ^4.4| a typical example of the monopole current at a time slice 
in the lattice QCD at /? = 2.4 in the maximally abelian (MA) gauge. In each 
gauge configuration, the monopole current appears as a distinct line-like object, 
and the neighbor of the monopole can be defined on the lattice. However, taking 
the temporal direction into account, the monopole current forms a global network 
covering over R^. 

Here, we summarize several relevant properties of fc^(s). 

1. The monopole current fc^ is topologically quantized and /c!f*(s) takes an integer 
A;^"*(s) G Z in the definition of Eq. ([4.21|) . As the result, A;^"*(s) forms a line-like 
object in the space-time R'^, since fc|f* is a conserved current as d^k^°^^ = 0. 
These features of fc!f*(s) G Z are quite unique and different from the electric 
current j^(s) G R, which can spread as a continuous field. 

2. In the lattice formalism, /c^"* = ^d^9a^ is defined as a three-form on the dual 
link. For the use of the forward derivative, k^^*{s) is to be defined on the dual 
hnk between s'^^^ = s + | + | + | + |±|. For instance, k\^*{s) is placed on the 



Figure 4.3: (a) The (static) monopole defined on the dual lattice is equivalent to 
the total magnetic flux of the Dirac string, (b) The neighboring links/plaquettes 
around the dual link. 

dual link between {s^ + l,Sy + |, s^ + |, St) and {s^ + l,Sy + |, Sz + |, St + 1). 
Thus, the monopole is deflned to appear at the center of the 3-dimensional 
cube perpendicular to the monopole-current direction as shown in Fig.[4.3|(b). 



3. Because of k^ = d^Faf^ = —\s^iai3')daFp.y, k^ only affects the perpendicular 
components to the /i-direction for the "electric variable" as Fq,^ in a direct 
manner. For instance, the static monopole with /cq 7^ creates the magnetic 
field Fij (i, j=l,2,3) around it, but does not bring the electric field Fq^. Hence, 
in testing the field around the monopole in the next chapter, one has to con- 
sider the difference between such perpendicular components and others. 

We now consider the relationship between the lattice variable and the field vari- 
able in the continuum theory. The continuous abelian field Afj,{x) = A^{x)y is 
expressed as 

eAl ^el-- (4.22) 

with the gauge coupling constant e and the lattice spacing a. The abelian field 

3 

strength T^y{x) = F^y{x)^ in the continuum theory is written as 



eF^, = mod2.(V) ■ ^ = V ■ ^. (4-23) 



9 A"^ 9 



and Ffj_,y is composed of two parts according to the decomposition ( |4.2CI| ) 



F^, = (c?AAV-^ur- (4.24) 



Figure 4.4: The typical example of the S-dimensional time-slice of the monopole 
current in the MA gauge in the lattice QCD with (3 = 2.4 on 16^. 

Thus, in the SU(A^c)-lattice formalism, the difference between the field strength 
Ff^i, and two-form {d A A)^^,^ arises from the periodicity of the angle variable in the 
compact subgroup 11(1)^'=^^ embedded in SU(A^c)- Here, the singular Dirac-string 
part -F^™^ is directly related to 2Tin^i, and is written by 

2 Am 1 
F^:^ = 27rn^, ■ — = —n,,-. (4.25) 



Owing to existence of Ff^^^ in Eq.( [4.24 ), the monopole current k^[x) = k 

T -7-2 

An An 



da*Fafi\ appears in the continuum theory and is written as 



^'x}^ 



3 



K = K ■ — = d^*ri^,-, (4.26) 

where the magnetic-charge unit g = — naturally appears in fc„. 



Chapter 5 

Large Field Fluctuation around 
Monopole 



In this chapter, we study the QCD-monopole appearing in the abehan gauge in 
terms of the gluon field fiuctuation [^3|, ^. For simphcity, we take Nc = 2. In the 
static frame of the QCD-monopole with the magnetic charge g, a spherical "magnetic 
field" is created around the monopole in the abelian sector of QCD as 

H(r) = -^r (5.1) 

with Hj = EijkdjAl. Thus, the QCD-monopole inevitably accompanies a large 
fiuctuation of the abelian gluon component At around it. As was discussed in 
the previous chapter, in the abelian gauge, the formal action of abelian projected 
QCD or the abehan part of the QCD action is given by S*^^^^ = — ^ / d'^x{{df^Al — 
di^A^)"^ — Ff^^^^}, where —Ff^^^ appears and eliminates the Dirac-string contribution. 
In the abelian part, the field energy created around the monopole is estimated as 
the ordinary electro-magnetic energy, 

£{a)= d^x-n{rf = ^, (5.2) 

where a is an ultraviolet cutoff like a lattice mesh. As the "mesh" a goes to 0, 
the monopole inevitably accompanies an infinitely large energy-fluctuation in the 
abelian part and makes 5"^^*^^ divergent. 

Since there seems no plausible reason to eliminate such a divergence via renormal- 
ization, the monopole seems difficult to appear in the abelian gauge theory controlled 
by 5'^^'^'. This is the reason why QED does not have the point-like Dirac monopole. 
Then, why can the QCD-monopole appear in abelian projected QCD ? To answer 
it, let us consider the division of the total QCD action 5"^'"° into the abelian part 
^Abei g^j^^ ^]^g remaining part 5'°''^ = S'^'^^ — S'^^'^', which is contribution from the 
off-diagonal gluon component. While S**^*"^ and 5'^'^'^' are positive definite in the 
Euclidean metric, 5*°*^ is not positive definite and can take a negative value. Then, 
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Figure 5.1: The local correlation between the abelian angle variable Of^{s) and the 
monopole current k^ in the MA gauge with U(l)3 Landau gauge at /? = 2.4. The 
closed symbol denotes the monopole current on the dual link. 

around the QCD-monopole, the abelian action 5'^tiei gj^guld be partially canceled by 
the remaining contribution 5'°''^ from the off-diagonal gluon component, so as to keep 
the total QCD action S^^^ finite even for a ^ 0. Similar cancellation between the 
gauge field and the Higgs field fluctuation is also found around the GUT monopole. 
Thus, we expect large off-diagonal gluon components around the QCD-monopole 
for its existence as well as a large field fluctuation in the abelian part. Based on this 
analytical consideration, we study the field fluctuation and monopoles in the MA 
gauge using the lattice QCD. 



5.1 Gluon Field Configuration around Monopoles 

We study the properties of monopole in terms of the gluon configuration in the 
MA gauge. In particular, we investigate the correlation between monopoles and the 
abelian angle variable Of^^s) and abehan projection rate -RAbei- 

First, we show the local correlation between the abelian angle variable and the 
monopole using the gauge configuration in the MA gauge with U(l)3 Landau gauge 
in Fig. |5.1| . The closed symbol denotes the monopole current on the dual link. 
Around the monopoles, the abelian angle variable tends to fluctuate largely. 
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Next, we consider the correlation between the field variables and the monopole 
quantitatively in the lattice QCD. For this argument, one has to recall the property 
of the monopole current shown in Section |4.3| . In particular, one should note that 
fc^(s) is defined on the dual link and only affects the perpendicular components to the 
/i-direction for the electric variable as Fa(3 because of k^ = d*^Fa^ = —^Sfxa/s-ydaFis.y. 
Taking account of these properties, we study the local correlation between the field 
variables and the monopole current k^{s) in the MA gauge with the U(l)3 Landau 
gauge. We first measure the average of the abelian angle variable Of^{s) over the 
neighboring links around the dual link (see Fig. [4.3| (a)), 

1 1 1 

l^'(5,A)l = T^E E 7;\^>.ap^\-\eUs + m(3 + n^)l (5.3) 

al3-ym,n=0 ^ 

which only consists of the perpendicular components considering the above monopole 
property. Here, the index fi denotes the direction of the dual link, and \9^{s,fi)\ 
corresponds to the average over the 12 sides of the 3-dimensional cube perpendicular 
to the /i-direction [45|. We show in Fig. |5.2| the probability distribution P(|^'^|) of 
\9^{s, fi)\ in the MA gauge with the U(l)3 Landau gauge at /3 = 2.4. The solid curve 



denotes P(|^^|) around the monopole current, while the dashed curve denotes the 
total distribution on the whole lattice. The abelian angle variable |^j^(s)| takes a 
large value around the monopole. In other words, the monopole provides the large 
fiuctuation of the abelian gauge field, which would enhance the randomness of the 
abelian link variable. 

Similar to |^^(s, fi)\, we measure the average -RAbei of the abelian projection rate 
i?Abei(s,/i) = cos6'^(s) over the neighboring links around the dual hnk, 

i?Abei(s,/i) = — E E T^l^fia^-rlcoseais + mP + n^) (5.4) 

in the MA gauge to investigate the correlation between off-diagonal gluons and 
monopoles. As shown in Fig.|5.3|(a), -RAbei around the monopole current becomes 



smaller than the total average of -RAbei and therefore the magnitude of the off- 
diagonal gluon component becomes larger around the monopole. The (5 dependence 
of the abelian projection rate (-RAbei) is shown in Fig.p73|(b). Although (-RAbei) on the 
whole lattice approaches to unity as /3 — >■ oo, (-RAbei) around the monopole is about 
0.88 and is not changed even in the large (5 region. Thus, the monopole provides 
the large fluctuation both for the abelian field and for the off-diagonal gluon. 

5.2 Plaquette Action Density around Monopoles 

We next study monopoles in terms of the plaquette action density. We define the 
SU(2), abelian and "off-diagonal" plaquette action densities as 

S'^'^^is) ^ l-^trnj^(^)(.), (5.5) 
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Figure 5.2: The solid curve denotes the probability distribution P(|^^|) of the aver- 
aged abelian angle variable \0^{s, jl)\ around the monopole current in the MA gauge 
with the 11(1)3 Landau gauge fixing. Here, \9^{s,jl)\ is the average of l^^l-^)! o'^^r 
the neighboring links around the dual link, For comparison, the total distribution P 
on the whole lattice is also added by the dashed curve. Around the monopole, |^^| 
corresponding to the abelian gluon component takes a large value. 



S: 



Abel/ 



sfM) 



1 - ^trnAbei(,) 



S']]^'\s)-S, 



AbcW 



(5.6) 

(5.7) 



where n^J/(^)(s) and n^^'^^(s) denote the SU(2) and the abelian plaquette variables, 
respectively; 



= u^{s)uu{s + fi)ul{s + i')ul{s). 






(5.8) 
(5.9) 



S*^^, instead of the Lorentz 



Here, all of S^^ are defined as symmetric tensors, S^^^ 
scalar, considering the above property of the monopole current. In the continuum 
limit a — > 0, S?^^'^\s) and Sf^^^\s) are related to the SU(2) and the abelian action 
densities as S^]^^'^\s) -^ ^a'^e'^tiGl,, and S{^}^'^^(s) -^ \a'^eHiJ^l„, and then we call 



HU 



S^y as the action density, in spite of the lack of the summation on the Lorentz 
indices. Here, S'^^ corresponds to the contribution of the off-diagonal gluon. While 
^'f^^'^^ s-^d S^^^^ are positive-definite, 5'°^ is not positive-definite and can take a 



negative value. 
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Figure 5.3: (a) The solid curve denotes the probabihty distribution P(^Abci) of 
the averaged abehan projection rate -RAbci(s, /i) around the monopole current in the 
MA gauge in the SU(2) lattice QCD with /3 = 2.4 on 16*^. For comparison, the total 
distribution P on the whole lattice is also added by the dashed curve, (b) The solid 
curve denotes abelian projection rate (-RAbei) around the monopole current in the 
MA gauge as the function of (5. The dashed curve denotes (RAhei) on the whole 
lattice. 
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In order to examine the correlation between the action densities and the monopole 
current defined on the dual link, we measure the average of the action density S{s) 
over the neighboring plaquettes around the dual link, 

1 ^ 1 

^is,!^) = -^Y^Y. 7;\^f^ap^\Sap{s + mj). (5.10) 

" a/37 ™=o 

Here, ft appearing in S{s,fi) denotes the direction of the dual link, and S{s,fi) 
corresponds to the average over 6 faces of the 3-dimensional cube perpendicular to 
the /i-direction. 

We show in Fig.^]^ the probability distribution P{S) of the action densities 
S{s,jl) in the SU(2), the abelian and the off-diagonal parts. Before the argument 
around the monopole current, we show the action densities on the whole lattice in 
Fig. |5.4| (a). On the whole lattice, most 5*°^ are positive, and both S^^'^^ and 5*°^ 
tend to take smaller values than S^^^"^^ = S^^'^^ + 5*°^. In other words, S^^^^ and 
positive 5'°''^ additionally contribute to S^^^'^\ 

However, such a tendency of the action densities is drastically changed around 
the monopole as shown in Fig. ^.4| (b). We find remarkable features of the action 
densities around the monopole as follows. 

1. Around monopoles, most 5*°^ take negative values, and S^^'^^ is larger than 

cSU(2) _ cAbel _|_ coff 

2. Due to the cancellation between S^^*^^ and 5*°^, S*^^^^-* does not take an ex- 
tremely large value around the monopole. 

Thus, the large abelian action density S^^"^^ around the monopole is strongly canceled 
by the off-diagonal contribution S°^ to keep the total QCD action S^'^^ = S^^'^^ + 
^off g]2iall. Here, different from S^^^'^\ ^Abei j^ggjf (Joes not control the system 
directly, and hence there is no severe constraint from S^^"^^. However, large S^^^^ 
is still not preferable, because the large-cancellation requirement between 5'^'^''^ and 
5*°^ leads to a strong constraint on the off-diagonal gluon and brings the strong 
reduction of the configuration number. 

Around the monopole, the abelian action density S^^^^ takes a large value, and 
this value can be estimated from a following simple calculation. Without loss of 
generality, the monopole- current direction is locally set to be parallel to the temporal 
direction as A;q"*(s) = ^da*Oao{s) = ±1. Here, kJ-^^^s) is expressed as the sum of 
six plaquette variables 6ij (i, j=l,2,3) around the monopole, because of k^Q^{s) = 
-j^eijkdi9jk{s) = -^ Ei Ei<fc £ijk{Ojk{s + i)- Ojk{s)}. Hence, the total sum of six 
\9ij{s)\{i < j) is to exceed 27r to realize ko{s) = ±1. Since large |%(s)| accompanying 
large S*^^"^' is not preferable, the magnetic field \9ij\ around the monopole is estimated 
as \6ij\ ~ 27r/6 = vr/S on the average, using the spherical symmetry of the magnetic 
field in the vicinity of the monopole. Accordingly, we estimate as Sff"^^ = 1 — 
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Figure 5.4: (a) The probability distribution P{S) of density S{s,fi) on the whole 
lattice in the MA gauge at /? = 2.4 on 16'^ lattice. (b)The probability distribution 
Pk{S) of the action density S{s,fi) around the monopole current k^. The dotted 
and the solid curves denote P(S'^^*-^-') and P{S^^'^^), respectively. The dashed curve 
denotes P{S°^) for the off-diagonal part S°^ of the action density. 
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Figure 5.5: The action density as the function of /5 in the MA gauge in the 
SU(2) lattice QCD. The closed symbols denote the action densities {S) around the 
monopole current, while the open symbols denote those on the whole lattice. The 
square, circle and rhombus denote (S*^^^^-*), (S'^*''^') and {S°^), respectively. The 
monopole accompanies a large U(l)3 plaquette action, however, such a large U(l)3 
action is strongly canceled by the off-diagonal part. 



cos(|6'jj|) ~ 1 — cos f = ^ around the monopole on the average. The above argument 
can be easily generalized to the case with arbitrary monopole-current direction. 

Then, existence of monopoles brings a peak around 5'^^'^' = | in the distribu- 
tion P{S^^^^). In fact, the abelian action density S^^"^^ has two ingredients; one is 
nontrivial large fluctuation about S^^"^^ = 1/2 originated from the monopole, and 
the other is remaining small fluctuations, which is expected to vanish as 5"^*^*^' — >■ 
as a — >^ 0. As shown in Fig. ^.5|, the peak originated from the monopole is almost 
f3 independent, while the other fluctuation becomes small for large (3. At a glance 
from this result, the monopole seems hard to exist at the small mesh a, since the 
monopole needs a large abelian action 5''^'''^'. Nevertheless, the monopole can exist in 
QCD even in the large (3 region owing to the contribution of the off-diagonal gluon. 
As shown in Fig. |5.4| (b), the off-diagonal part 5*°^ of the action density around the 
monopole tends to take a large negative value, and strongly cancels with the large 
abelian action 5'^^''^ to keep the total SU(2) action S'^'^^ finite. 

Here, we consider the angle variable Xm(^) ^^ the off-diagonal gluons A^(x) 
around the monopole. In the MA gauge, the amplitude of A'j^lx) is strongly reduced, 
and x^{x) can be approximated as a random variable on the whole, because Xm(^) 
is free from the MA gauge condition entirely and is less constrained from the QCD 
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action due to the small |A^(a;)|. However, around the monopole, the off-diagonal 
gluon v4j^(x) inevitably has a large amplitude even in the MA gauge to cancel the 
large abelian action density. This requirement on the reduction of the total action 
density severely constrains the randomness of the angle variable Xfj,{x) of the off- 
diagonal gluon A^(a:) around the monopole. As the result, the randomness of X^{x) 
is weaken, and continuity of Xni^) or ^t{^) becomes clear in the vicinity of the 
monopole even in the MA gauge. This continuity of A^{x) around the monopole 
ensures the topological stability of the monopole itself as n2(SU(2)/U(l)) = Zoo. 

To summarize, existence of the monopole inevitably accompanies a large abelian 
plaquette action 5*^^*^' around it, however, the off-diagonal part 5*°^ takes a large 
negative value around the monopole and strongly cancels with 5''^'''^' to keep S^^^ 
not so large. Due to this strong cancellation between S^^'^^ and 5'°''^, monopoles 
can appear in the abelian sector in QCD without large cost of the QCD action 
^QCD^ which controls the generating probability of the gluon configuration. The 
extension of the off-diagonal rich region around the monopole can be interpreted as 
the effective size or the structure of the monopole, because the abelian gauge theory 
is largely modified inside the QCD-monopole like the 't Hooft-Polyakov monopole. 

Finally, in this section, let us consider the correlation between monopoles and 
instantons [^] in terms of the gluon-field fluctuation. The instanton is a nontrivial 



classical solution of the Euclidean Yang- Mills theory, corresponding to the homotopy 
group Ha (SU(A^c)) = Zqo [0,11, Hi- For the instanton, the SU(2) structure of the gluon 
field is necessary at least. In spite of the difference on the topological origin, recent 
studies indicate the strong correlation between monopoles and instantons in the 
QCD vacuum in the MA gauge |3^, ^ . What is the origin of the relation between 
two different topological objects, monopoles and instantons ? In the MA gauge, off- 
diagonal components are forced to be small, and the gluon field configuration seems 
abelian on the whole. However, even in the MA gauge, off-diagonal gluons largely 
remain around the QCD-monopole. The concentration of off-diagonal gluons around 
monopoles leads to the local correlation between monopoles and instantons: instan- 
tons appear around the monopole world-line in the MA gauge, because instantons 
need full SU(2) gluon components for existence. 



Chapter 6 

Monopole Projection and Scaling 
Properties in the MA Gauge 



6.1 Decomposition into Monopole and Photon Angle- 
Variables 

Abelian projected QCD is obtained by neglecting the off-diagonal gluon component 
in the abelian gauge, and it includes both the electric current j^ and the monopole 
current A;^. Here, j^ is generated by charged gluons, and k^ is generated by the 
singular SU(2) gauge transformation in the abelian gauge. In the lattice formalism, 
the U(l)3-gauge field can be separated into the "monopole part" and the "photon 
part" using the Coulomb propagator. Here, the monopole part only includes the 
monopole current fc^ and reproduces the purely linear part of the static quark po- 
tential ||5^. On the other hand, the photon part includes the electric current j^ only, 
and provides the Coulomb potential between the quark and the anti-quark, similarly 
in the ordinary QED. In this section, we extract the monopole contribution from 



abelian projected QCD using the lattice Coulomb propagator |^, and study the 
role of the monopole for confinement. 

The abelian gauge field 0^{s) can be decomposed into the monopole part 9^°{s) 
and the photon part ^^^(s) 



e;f°(s) = 2nY^{s\d-'\s')d^n^,{s'), (6.1) 

s' 

C(^) = T.{^\d'"W)dJa,{s% (6.2) 

s' 

where {s\d~'^\s') = —^jjz^ is the Coulomb propagator in the 4-dimensional Eu- 
clidean space. As for the U(l)3 gauge invariance, 9^°{s) is gauge- variant, since 
27171^1, is gauge- variant. On the other hand, Oj^^{s) is U(l)3 gauge-invariant, since it 
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is composed of gauge- invariant 9^^. In the Landau gauge <9^6'^(s) = 0, one finds 

Of''is) + 0l\s)=el{s). (6.3) 

Here, we investigate properties of the monopole and the photon parts in terms 
of the electric and the magnetic currents. The two form of these angle variables are 
written as 

(SA^^V =C + 2™;t° with ^Jt° = niod2.(9A^^Ve(-vr,7r] 

(9aOm- =C + 2™|:' with 6l^ = mod2AdAe''''),,e{-n,7i], 

using n^°, n^^ G Z. The monopole current k^{s) and the electric current j^{s) in 
the monopole part are defined as 

> Mo ^ J_g*mo ^ _A*^Mo -Mo ^ Q qMo /^ .^ 

and those in the photon part are defined as 

C - ^^:c; = -^>S' jT - ^"C- (6-5) 

From the physical point of view, the decomposition of the abelian sector into the 
photon and monopole sectors directly corresponds to the separation of the electric 
current j^ and the monopole current k^ near the continuum limit. 

In the actual lattice simulation, we can observe k^ ~ 0, j^ = j^ and fcjf" ~ 0, 
k^° = kfj_ in the MA gauge within a few percent error for large (3. For instance, 
we show in Fig. |6.1| (a) the lattice result of the monopole current density pj^c — 
Es^. \k^is)\/ EsM 1 in the MA gauge. One fined p^^j' - Pmc and p^c - 0- 

Existence of the monopole accompanies the Dirac string. We show also the 
Dirac-sheet density pas = J2sfiu \^ij.{.^)\/ J2siiu 1 ^^ the MA gauge with U(l)3-Landau 
gauge in Fig.|6.1|(b). Owing to the U(l)3-Landau gauge fixing, the abelian angle 
variable ^j^ becomes mostly continuous and |(<9 A 9^)fj,i^\ scarcely exceeds vr, so that 
redundant Dirac sheets are eliminated and the correlation between p'^Q and p^s 
appears. More quantitative argument of p^^ and pjfg will be done in the physical 
unit in Section p^ . 

For comparison, we investigate the monopole density pj^^ in the SU(2) Lan- 
dau gauge with SU(2)/U(1)3 fixed, which is a kind of the abelian gauge. In the 
SU(2) Landau gauge, each component of the gluon field is maximally smooth, and 
monopoles scarcely appear because they are generated from the singular-like large 
fiuctuation of the abelian field. For instance, p'^Q in the SU(2) Landau gauge in less 
than 1/10 of pj^^ in the MA gauge at /3 = 2.4. In fact, in the regular SU(2) gauge 
configuration, the gluon fiuctuation is shared into each component almost equiva- 
lently, and there appears no singularity as the monopole due to the "wise sharing" 
of the fluctuation. On the other hand, in the MA gauge, large gluon fluctuations are 
concentrated into the abelian sector by the suppression of off-diagonal gluon com- 
ponents, and therefore the monopole appears as the singularity or the topological 
defect from the large fleld fluctuation in the abelian sector. 
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Figure 6.1: (a) The nionopole-current density pj^^ in the MA gauge as the function 
of p. The circle, square and triangle denote p^^Q in the abelian, monopole and 
photon sectors, respectively, (b) The Dirac-sheet density p^g in the MA gauge with 
U(l)3 Landau-gauge as the function of p. The circle, square and triangle denote 
Pos in the abelian, monopole and photon sectors, respectively. 
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6.2 Link and Plaquette Variables in Monopole and 
Photon Sectors 

In general, confinement observed as the area-low reduction of the Wilson loop is 
closely related to the large gluon fiuctuation as is suggested by the strong coupling 
QCD. Hence, from monopole dominance for the confinement force observed on the 
lattice, one may expect large fiuctuation of the gluon field in the monopole sector 
and small ones in the photon sector. This naive expectation on the large fiuctuation 
9^° may be also suggested from the definition of 9^° and 9j^^ in Eq. ( |6.3|) , since the 
monopole angle variable 6'ff° consists of the integer part 27m G 27rZ, while the photon 
angle variable 9^^ consists of the fractional part ^^ G [— tt, tt). We measure the lattice 
angle variables 9^{s), 9^°{s) and 9^^{s) in the MA gauge with U(l)3 Landau gauge 
as shown in Fig. |6.2| . On the whole shape of the probability distribution P{9^) of link 
variable, the distribution of monopole part P{9^") is smaller than that of the photon 
part P{9^)^ and P{9f^) is similar to P{9f^^) rather than P{9j^°). These results are in 
contradiction to the above expectation on the large fiuctuation of P{9^°), and seem 
to be surprising in terms of the confinement properties, because the string tension 
in the abelian sector resembles that in the monopole sector and the string tension 
in the photon sector is almost zero. 

Quantitatively, the whole shape of probability distributions P{9^), P{9^°) and 
P{9f^^) can be fitted with the Gaussian curve around 6* = as 

P{9) = ^e-^^^ (6.6) 

with a =0.167 for Pi9f^), a =0.038 for P{9f°) and a =0.114 for P{9l'^). 

Now let us reconsider the fiuctuation of the monopole part. In the U(l)3 Landau 
gauge, the appearance of the Dirac sheet with n^j, 7^ is strongly suppressed, 
because the field is forced to be continuous there. Therefore, almost all of the two 
form 9^° = d^9^° - d^9^° satisfies -vr < 9ji° < vr, i.e., one finds n]f°{s) = almost 
everywhere on the lattice. This is the reason why the fiuctuation of the monopole 
angle variable 9j^° becomes small. However, there appears the relic of the large 
fiuctuation of the abelian angle variable in the monopole angle variable, as shown 
in Fig. |6.3[ Similar to P{9l), the distribution P{9^°) of the monopole angle variable 
has a non-Gaussian large tail even for \9^°\ < |, while P{9j^^) has the Gaussian 
tail only and mainly distributes within — f < 9^^ < | as shown in Fig.|6.3|. Such a 
large fiuctuation of 0ff° is responsible for the appearance of the monopole, and it 
also would play a relevant role for confinement properties, because a large reduction 
of the string tension is observed in the lattice QCD by the artificial elimination of 
the large fiuctuation. 

Next, we consider the plaquette action densities of the monopole and the photon 
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Figure 6.2: The probabihty distribution P{6^) of lattice angle variables 6*^^(3), 
6'Jf°(s), 0^{s) G (— 7r,7r] in the MA gauge with the U(l)3 Landau gauge fixing at 



j3 = 2.4 on 16^ lattice, (a) the abelian angle variable P 
angle variable P{6^°). (c) the photon angle variable P{Of^^ 
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Figure 6.3: The tail of P{dl), P{ef°) and P{ 
a large fluctuation in^ 
fluctuation ingredient. 



Both P(^3) and P{9^°) include 
a large fluctuation ingredient as 6^ > |, while P{6^^) does not have such a large 



^,^'^ 



parts, 



C(^) 



S^^is) 



^ l-^trnMo(,) 



(6.7) 
(6.8) 



similar to the SU(2) and abelian action densities in Eqs.( |5.5|) and ( ^.6|) . Her" '~'^° 
and n^^(^s) are the plaquette variables in the monopole and photon parts 



:6D. Here, n^^-ls) 









S + /iX°t(5 + ^)^Mot(^, 



S + fl)u 



Phf, 



s + z>)mP^^V5) 



(6.9) 
(6.10) 



respectively. To see the monopole effect on the action density Sf^^{s), we measure 
the average S{s,fi) over the neighboring plaquettes around the dual link, 



1 ^ 1 

" a/37 '"=0 



(6.11) 



similar to Eq.( |5.10|) . As shown in Fig.|6.4|, both action densities P{S^°) and P{S^^) 
have a peak near S* = 0. For the monopole parts, there are two ingredients in the 
probability distribution P{S^'^) of the monopole action density S^°{s,fi): one is 
the large fluctuation around 5*^° 



I and the other is the small fluctuation near 
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^Mo ^ g_ rjj^g large fluctuation ingredient corresponds to the average action-density 
over six plaquettes of the cube including the monopole inside 

C = 1 - f^'^'i: = 1 - ^°« 'C ^ 1 - ^°^ Y = ^' (6-12) 

which is similar to the abelian action density P[S^^'^^) around the monopole, as 



shown in Fig. |5.4| . Here, 27r is the total magnetic flux created by the unite charge 
monopole. On the other hand, the small fluctuation ingredient corresponds to the 
plaquette apart from the monopole. For the photon part, the global shape of the 
probability distribution P{S^^) of the photon action density S^^{s^fi) resembles 
p^^Abei^ as shown in Fig.|5.4|, however, P{S^^) does not have the large fluctua- 



tion peak around S^^ = |, corresponding to the absence of the monopole unlike 
pj^^Abei-j_ jj^ ^Yie distribution P{S^°), the peak position around 5"^° = ^ is almost 
/^-independent, since it is geometrically determined as Eq.( |6.12|) . On the other hand, 
the peak near S = becomes narrow both in P{S^°) and in P{S^^) as (3 increases. 
Finally in this section, we consider also the monopole effect on the discontinu- 
ity of the angle variables. As the result of the large fluctuation of the monopole 
plaquette variable around the monopole as 5''^° ~ |, there appears the discon- 
tinuity of the link variable 6^°. We measure the correlation neighboring links 
u^{s) = e*^^*-* and u^(s + 0) = e^^^'^~^'^\ We show in Fig. |6.5| the angle difference 
A = |mod27r(6'^(s) — 9^{s + z>))| with ji ^ u. The photon angle variable 0^^{s) tends 
to be continuous with the only small difference between the neighboring links as 
A < |. This tendency on the continuity of 0^^{s) becomes clear as [3 increases. On 
the other hand, similar to the abelian sector, the monopole angle variable 0^°{s) 
includes the discontinuity as A > |. Thus, the appearance of point-like monopoles 
in the MA gauge brings the large fluctuation and the discontinuity of the angle 
variable O^^is) around them. 



6.3 Dual Field Formalism 

In the presence of magnetic monopoles, the abelian gauge fleld A^{x) inevitably 
includes the singularity as the Dirac string, which leads to some difficulties in the 
fleld theoretical treatment. Here, we investigate the description with the dual fleld 
formalism, which is useful to describe the monopole sector in the QCD vacuum. In 
the dual formalism, the dual gauge fleld B^{x) is introduced as F"^^^^ = d^B^ — 
dyB^ with -F^"'^^ =* F^y = ^e^„ai3F°''^ , and the interchange between A^ and B^ 
corresponds to the electro-magnetic duality transformation, F^^, ^^ F'^^'^^ or H <->■ E 
||33|| . The monopole sector, which carries essence of the nonperturbative QCD, 
includes the magnetic current only and does not include the electric current: fc^ 7^ 
and j^ = 0. Therefore, the monopole sector is the dual version of the ordinary QED 
system with j^ 7^ and A;^ = 0. The dual gauge fleld B^{x) can be introduced 
without the singularity like the Dirac string owing to the dual Bianchi identity. 
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Figure 6.4: (a) The probability distribution P{S ) of the monopole part of the 
action density S^^°{s,fi) around the monopole. The dotted curve denotes S^°{s,fi) 
on the whole lattice. The value of S^°{s,fi) is the average of the action density 
over the neighboring plaquettes around the dual link, (b) The probabil- 
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ity distribution P{S^^) of the photon part of action density S^^{s,fi) around the 
monopole. The dotted curve denotes S^^{s,jl) on the whole lattice. 
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-d^T, 



dual 



—jy = 0, corresponding to the absence of the electric 



current. 

Let us consider the derivation of the dual gauge field B^{x) from the monopole 
current k^{x), taking the dual Landau gauge, (9^-B^ = 0. Then, the relation 



dT, 



dual 

fjLV 



d By — dy{d^B^) = ky becomes a simple form d B^ = k^. Therefore 



starting from the monopole current configuration k^[x), we derive the dual gauge 



field B^{x) as 



B,{x) = / dS{x\d-^\y)k,{y) 



47r2 



d^y- ^^^y^ 



{x - yf 



(6.13) 



Using the dual gauge field 5^, the Wilson loop in the monopole sector is expressed 
as 



iy^° = exp{ieidx''Af°} = exp{ie f fda''''*Fjfj'} 



exp[—te 



da'"'*FI^;^^^) 



exp(-ie f f da^"'*{dAB] 



^Ujl 



(6.14) 



where the Stokes theorem is applicable because of abelian nature. 

Now, we apply this dual field formalism to the monopole sector in the lattice 
QCD in the Euclidean metric. The lattice dual gauge field 6*^"^' is defined as 6*^^ 
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Figure 6.6: The comparison of the probability distribution P{9'^^^) with P{9^°) in 
the monopole sector in the MA gauge. The dual gauge variable O"^^"^ is defined in 
the dual Landau gauge d^Of^'^^ = 0. The angle variable ^Jf° in the monopole sector 
satisfies the Landau gauge condition dn9^° 



^t^" fj. 



0. 



aeB^/2 similar to 6*^ = aeA^/2. In the dual Landau gauge, 9'^^^^ is obtained from 
the monopole current fc^ as 



e^-% + f,) = 2nY^{s\d''\s')kl^\s') 



(6.15) 



using the lattice Coulomb propagator {s\d ^|s') = {s\{d^d'^) ^\s'), where (9^ and d'^ 
denote the forward and backward derivatives, respectively: 



a;/(x) = /(.)-/(. -A). 

The two-form of 6'^]^^' is defined by 



oiT\s) 



I fldual/ \ o/ fldual/ 



&A 



'v^i 



et^^ + 27rnf"^i 



[iV 



iLV 1 



where the dual abelian field strength Q^J^"^ = ea?F^1^'^^ (2 is defined as 
d^;f{s) = mod2.(9 A e^"^%.(s) e (-TT, tt] 

in the dual gauge invariant manner. 
], w 



(6.16) 
(6.17) 

(6.18) 
(6.19) 



In Fig.|6l6|, we compare the probability distribution P{0'^^^) of dual gauge variable 



^duai ^^^ p^0Mo^ q£ ^YiQ gauge variable ^ff° in the monopole sector in the MA gauge. 



72 



6 Monopole Projection and Scaling Properties in the MA Gauge 



§. 



CD 




Figure 6.7: The probabihty distribution P{Or^ ) of the dual abelian field-strength 
^j]"^' in the monopole sector in the MA gauge. Here, PiOf^^"^) exactly coincides with 



dual 



0. 



and 6'J^° is defined in 



Eq. (|6.1| ) satisfies the Landau gauge d^6^° = 0. These two distributions seem to have 



Here, 9'^ is defined in the dual Landau gauge 9^6' 

almost the same widths. This is because the Dirac sheet 71^,^ is strongly suppressed 
in the Landau gauge, while there is no dual Dirac sheets with n^^'^^ 7^ in the 
dual Landau gauge owing to the absence of j^. In fact, only the essential minimal 
fiuctuation remains both for the monopole variable 6^° in the Landau gauge and 
for the dual variable 0'^^'^^ in the dual Landau gauge. Different from 6^°, however, 
^duai Y^^^ ^^Q bumps around 6'^'^^' = =1=^ originated from the monopole current. 

We show in Fig.^]^the probability distribution P{9'^^^^) of the dual abelian field- 
strength ^j]"'^'. There appear also some bumps originated from the quantization 
of the monopole current and lattice discretization. Here, P(^^"^') is the same as 
PiOp) because of Of^'^''^ = *ef° in the continuum limit 



fii/ 



between n'^^^ and 



However, the two-forms 
^^ =* (9 A 9^")fj,^, according to the difference 
"nj^". In the description by 6}^°, there inevitably appears the 
singularity as the Dirac sheet with n^jy 7^ to generate the monopole current through 
the breaking of the Bianchi identity, although the appearance of the Dirac sheet is 
strongly suppressed in the Landau gauge. In the dual formalism, on the other hand, 
the monopole sector becomes regular and does not include the dual Dirac sheet, i.e. 
n^"*^' = 0, in the dual Landau gauge owing to the dual Bianchi identity. Actually in 
the monopole sector in the dual Landau gauge, 6*1]"^' = (d A 6''^"^')^,^ distributes only 
in the region of -n < 6'^"^' < vr, and hence one fines i9^]^^' = mod27r((?A6''^''^')„„ = 9'^''^^ 



fiu 



Hiy 
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and n*^^^ = everywhere on the lattice. Thus, the monopole sector in the lattice 
QCD is described by the dual gauge field without the singularity as the Dirac string 
in the dual Landau gauge. 

Using the dual gauge variable 6*^"^', the static potential between the monopole 
and the anti-monopole is obtained from the dual Wilson loop ti Yl]^^iu'^^^{si} with 



.dual/ 



M"""'(s) = exp{ir30']"^'(s)} [^. The large dual Wilson loop obeys the perimeter law, 



and the inter-monopole potential can be fitted by the Yukawa potential including 
the monopole size effect |5^. Thus, the dual gauge field 5^ seems to acquire its 



mass in the infrared region, which would be a direct evidence on the dual Higgs 
mechanism by monopole condensation occurring in the QCD vacuum. 

6.4 Scaling Properties on Monopole Current and 
Dirac Sheet 

In this section, we study the /5-dependence of the gluon field and the field strength in 
the abelian, monopole and photon sectors in the MA gauge with the U(l)3 Landau 
gauge. We study also the monopole-current density pmc and the Dirac-sheet density 
Pds using the lattice QCD simulation with various (3. We argue these quantities both 
in the lattice unit and in the physical unit, and examine their scaling property. 

In the MA gauge with the U(l)3 Landau gauge, the probability distribution of 
abelian angle variables, Of^, 6^° and 6j^^, also exhibit a peak around 6'^ = as was 



shown in Fig.|6.2|. We show in Fig.|6.8|(a) the /3-dependence of (16*1^1), (I^Jj'^l) and 



{\0jj^°\) in the MA gauge with the U(l)3 Landau gauge, and find that the fiuctuation 
of the abelian angle variables decreases with /3 on the average. On the continuum 
abelian fields, eA. = ^6^, eA^J" = ^6^"^ and eA^° = ^6^°, we show in Fig.p|(b) 
the /3-dependence on (leA^I), (leAj^^l) and (|ey4|f°|). While there appears a /?- 
dependence on eA^ and eA^^, the continuum field eA^° in the monopole sector is 
also found to be less /3-independent. 

In the MA gauge, the probability distribution of the abelian plaquette variable 



9fj_iy, has a peak around 6'^,^ = as was shown in Fig. |3.11| , and hence the abelian 
plaquette action seems to control the abelian system like the SU(2) total action. 
Such a tendency on the gathering around ^^i^ = becomes clearer as (3 increases. 
Quantitatively, we show in Fig.pl9|(a) the /3-dependence of (|^^,y|), {\0j^}^\) and (|^^°|) 
in the MA gauge, and find decreasing of all the plaquette angles with (5. From 
these lattice data, we try to extract the continuum abelian field strength, eF^,^ = 
^e^,, eFj:^ = ^ep and eF^J^ = ^9^°. We show the /3-dependence on (|eF^,|), 
{\eF^^\) and (|eF^°|) in Fig.|S.9|(b), where the lattice constant a used is determined 
so as to reproduce the string tension cr=0.89GeV/fm in the SU(2) sector. While 
there appears a /3-dependence on eF^,^ and eFJ^^, the continuum field eF^° in the 
monopole sector is almost /?- independent. In other words, we find a good scaling 
property even for the local variable as eFj^° in the monopole sector. 
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Figure 6.8: (a) The absolute value of the lattice angle variables 9'^{s), 9 

in the MA gauge with U(l)3 Landau gauge fixing on 16'' lattice as the function of p. 
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in the MA gauge on 16^ lattice as the function of (3. These are obtained from 



431 



9^ = Af^- Y with lattice spacing a. The lattice spacing a is defined so as to produce 



the string tension a =0.89 GeV/fm. 



6.4 Scaling Properties on Monopole Current and Dirac Sheet 75 



CD 



0.8 



0.6 



0.4 



0.2 



0.0" L. 



(a) 



A- - _ _ 



-'-—--■::-:: 



2.2 2.3 2.4 2.5 



2- 



> 

o 

> 

tin 



1 - 















.••■;-^ 




.•■■••■"- 




..•■••■"■- 


:•'*'' 


.;-;f' 




• 


■■,-A- 












if ' 














- 












■ 














(b) 



2.2 



2.3 



2.4 



2.5 



P 



Figure 6.9: (a) The absolute value of the abelian field strength ^p,v\ in the MA 
gauge as the function of /5. The circle, square and triangle denote ^^v\ in the abelian, 
monopole and photon sectors, respectively, (b) The absolute value of the continuum 
abelian field strength \eFy,J\ = |^^,^| ■ ^ in the physical scale as the function of fi. 
The circle, square and triangle denote |e-F^jy| in the abelian, monopole and photon 
sectors, respectively. 
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In order to understand the good scahng property on the abehan continuum 
variables in the monopole sector, we examine the /3-dependence of the monopole- 
current density pmc and the Dirac-sheet density p^s in the MA gauge with the U(l)3 
Landau gauge. In the physical unit, the monopole- current density pmc is expressed 
as the ratio of the total monopole-current length L\„^'^ to 4-dimensional physical 
volume V^ot^* of the system, 

dat ^ I MV /I jvhys 



j)hys _ PmC _ f^H — tot 

,3 - ^^4 - yphys- 



vnys _ rivnj _ ^^H _ t°t /^^ r,p,N 

PmC = -^ - ^^"1 - T rvhvs ■ [O.ZU) 



S,fJ, 

Similarly, the Dirac-sheet density pDS in the physical unit is expressed as the ratio 
of the total Dirac-sheet area S^^t'^ to 4-dimensional physical volume Vtot^'* oi the 
system, 

lat ^-^ \ ii.v\ )\ qphys 

Pds - ,, - J, ^4 -y,n,. (6.21) 



Jihys _ PdS _ ^./^.'^ _ '^tot' 

d 



If the monopole-current length L^ot^ and the Dirac-sheet area Slot'^ are physical 
values in the continuum limit, Pmc** and p^^'^ are to be almost /5-independent as long 
as the mesh a is small enough. We show in Fig.|6.10| the /^-dependence of Pmc'* and 



Pj^g''' in the MA gauge with the U(l)3 Landau gauge. As for the monopole density 
pyi(^ -I there is a /5-dependence to some extent. On the other hand, the Dirac-sheet 
density p^s^* i^ almost /5- independent, and exhibit a clear scaling property. 

Thus, n^i^^^s) = n^u{s)/a? is expected to be almost /5-independent on the av- 
erage over the whole system. Near the continuum limit, we can then expect a- 
independence of eA^° as 



e/i;f°(s) = -^^(s)^° = -Y,D'^\s - s')d':W,{s') ^ jd'x'D{s - x')d^nl 

Hi Hi / fJ 



where the last expression is to be a-independent. Thus, the good scaling prop- 
erty of the Dirac-sheet density is considered as the origin of the approximate (3- 
independence of the abelian continuum variables eA^° and eF^° in the monopole 
sector. 

It is noted that the Dirac-sheet area S'^t^* seems to have a continuum limit in 
a fixed 4-dimensional physical volume Kot^'*) although the monopole-current length 
Lfojf' tends to increase as the mesh a decreases. Here, let us consider the reason of 
the good scahng for the Dirac-sheet area in the U(l)3 Landau gauge. In the U(l)3 
Landau gauge, the U(l)3 gauge variable ^^(s) is forced to be mostly continuous as 
^j^ ~ on the lattice, which suppresses the large 6^^, fluctuation which generates 



rinu 7^ as was shown in Fig. 3. 10. Thus, the generation of the Dirac sheet is 
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Figure 6.10: 



The monopole-current density (p^ ^^^ 



^MC 



and the Dirac-sheet density 



(Pds'*)^ as the function of (3 in the MA gauge on 16'* lattice. The area of the Dirac 
sheets is /^-independent, while the monopole-current length slightly increases with 

p. 



strongly suppressed in the U(l)3 Landau gauge, although the Dirac sheet induced 
by the monopole current must appear. Geometrically, the area of the Dirac sheet is 
minimized in the U(l)3 Landau gauge, within the constraint that the boundary of 
the Dirac sheet is fixed as the monopole current. This situation resembles the "soap 
bubble" with the boundary fixed. Thus, the Dirac-sheet area is insensitive to the 
small fluctuation of its boundary or the monopole current due to the minimizing 
condition of the area in the U(l)3 Landau gauge. 

On the lattice, as the mesh a becomes small, the monopole current may be 
fractal, and the length of the monopole current may become large with small a. 
Nevertheless, the Dirac-sheet area remains almost unchanged in the U(l)3 Landau 
gauge owing to the minimizing condition of the area, even when its boundary or the 
monopole current becomes fractal. This would be the advanced point of the use of 
the Dirac sheet in the U(l)3 Landau gauge for the scaling property rather than the 
monopole current. 



6.5 Estimation of the Wilson Loop from the Monopole 

Up to now, we have found that the monopoles appearing in the MA gauge brings 
the large fluctuation of the abelian field variable around them, which would be 
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the origin of the string tension. In this section, we consider the relation between 
the monopole density pmo and the string tension o in order to clarify the role of 
monopole for the confinement force in the analytical manner. Here, we idealizes 
the large field fluctuation created by the monopole, and pay attention to the area- 
law reduction of the Wilson loop (VT) . For the system with the monopole density 
Pmo, one (anti-) monopole is expected to occupy the area of the cube with the 
length ^Mo = Pmo ' which would be a typical scale of the confinement physics |]58 |. 



Here, tmo = ^f^ corresponds to the size of the QCD monopole, and is estimated 
as rMo = 0.25 ~ 0.3fm from the lattice QCD with (5 = 2.4 as shown in Fig.|0 



As an important feature, the unit-charge monopole provides a minus factor to the 
Wilson. Let us consider the static monopole with the unite charge and the "spatial" 
Wilson loop surrounding the monopole as shown in Fig. |6.11| (a). Taking the plane 
sheet S as the surface area inside the loop, one can use the simple estimation of the 
monopole effect as a "minus factor" , 



-3 ^ ^ _ _3 



expjiey rfs^Ajy} = exp{iejjda^uF^^^} = exp{i2nx—} = exp{inT^} = -1, 

which corresponds to the half oi the quantized magnetic flux between the monopole 
and anti-monopole as shown in Fig. |6.11| (b). In general, unit-charge monopole cur- 
rent which perpendicularly penetrates to the Wilson loop W provides the "minus 
factor" to W. 

On the other hand, the effect of a monopole apart from the Wilson loop is 
expected largely canceled by the nearest anti-monopole, and is not expected to 
provide a large contribution to the Wilson loop in the dense monopole system like 
the AP-QCD vacuum. Such a screening effect is also suggested by recent lattice 
simulations on the inter-monopole potential ^^. Here, we try to estimate the 



contribution of monopoles to the Wilson loop of the R x T rectangular, where R 
and T are physical lengths. Defining the plane area S inside the loop, we expect 
that only nearest monopoles to 5* can contribute to the Wilson loop without the 
screening effect by other monopoles. Then, only monopoles inside RxT x duo plate 
area near S are provided the minus factor to the Wilson loop. The total monopole 
number in the RxT x dyio plate area is estimated as 

RxT 
Nmo = {RxT X duo) X pmo = —72 — • (6.22) 

Introducing the lattice with a fine mesh a <^ duo, we try to estimate the Wilson 
loop as the function of duo- Since the lattice plaquette number on S is A^'"* = ■^^^, 
we divide the plate area V into A^'"* "cells" as a x a x ^mo boxes. When a is taken 
small enough, each cell includes at most one (anti-)monopole, and the existence 
probability q of the monopole inside one cell is obtained by g = ^^ = -^ <C 1. 



Mo 



Next, we apply the statistical consideration for the estimation of the Wilson loop 
{W). The probability that n cells include the monopole among A^'"^* cells is simply 
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Figure 6.11: (a) The monopole in the plaquette. In the lattice unite, the magnetic 



field surrounding the monopole is d 



total 



27r, and therefore the magnetic field pene- 



trating the north or south half sphere is 9^i, = vr. The contribution of the monopole 
to the Wilson loop is e*" '^•'"'^''^ = e*^^'''^^ = e*''^^ = — 1. (b) The magnetic flux, 
which has positive and negative charges at the ends, penetrated in the plaquette. 
The quantization condition of the flux requires to 6^i, 

2pf 12. 
^Jm" 9. 



contribution of this flux to the Wilson loop is e*'' '^■''"'' 2 = q^i^^'^s 



21111 with n G Z. The 



given ]\fiatCnq"'{l — Q')^" " in a statistical source. Since n monopoles provide (— ) 
factor to the Wilson loop W, the expectation value {W) is estimated as 



(W) 



Aflat 

^ ^latCnq'^il — q) 

n=0 



N'-°-*-n (V"- 



Ajlat 

22 Af;«tC'„g"(l — q) 

ra=0 



:i - 2g) 



Ajlat 



ATlat 



exp{A^'"*ln(l - 2g)} = exp{:5-^ln(l - 2g)}. 



(6.23) 



Therefore, the Wilson loop obeys the area law, and the string tension a is obtained 
as 



a 



RxT 



^ ln(W^) = -^ln(l-2g) = -^ln(l-2 ° 



"Mo 



"Mo 



(6.24) 



because of g = -^ <^ 1. Thus, we derive a simple relation between the string 
tension a and the inter-monopole distance duo, 

V^-V2dul (6.25) 
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As shown in Fig. |6.1(J| , we have found duo = Puo = 0-5 ~ 0.6 fm from the lattice 
result with /5 = 2.2 ~ 2.5. Then, the string tension is evaluated as y/a ~ ^/2d^^ = 
(0.42 ~ 0.57) GeV from the monopole density p using this estimation. 



Chapter 7 

Monopole Current Dynamics 



In the infrared region, color confinement can be understood by the dual Higgs the- 
ory with monopole condensation. The string tension is almost reproduced by the 
monopole sector of the abelian gluon component, which is called as monopole dom- 
inance for the color confinement force. Thus, the monopole degrees of freedom 
would be the key variable for the confinement physics. In the dual Higgs theory, 
the monopole is assumed to be condensed, which has been suggested by the forma- 
tion of global network of the monopole current in the lattice QCD. In this chapter, 
we examine the realization of monopole condensation in the QCD-vacuum at the 
infrared scale. However, QCD is described by the gluon field not by the monopole 
current, and therefore it is difficult to clarify monopole condensation only with the 
lattice QCD simulation. To this end, we generate the monopole-current system 
on the lattice using a simple monopole current action and study the role of the 
monopole current to the color confinement |p9| . Comparing the QCD vacuum with 
the monopole current system, we try to clarify the monopole condensation occurred 
in the non-perturbative QCD. 

7.1 Monopole Current Dynamics and Kosterlitz- 
Thouless-type Transition 

To begin with, we consider the monopole-current action S[kfj_{s)] which simulates the 
monopole-current system observed in the lattice QCD in the MA gauge. In general, 
S[k^{s)] includes the nonlocal Coulomb interaction as Sc = J d'^xd'^yk^{x)D{x — 
y)k^{y) with the Coulomb propagator D{x). In the dual Higgs phase, however, the 
effective interaction between the monopole currents would be short-range due to the 
screening effecthj the dual Higgs mechanism |^, 5^ similar to the Debye screening 



60|| . Then, the partition functional of the monopole current A;^(x) = kl^{x) ■ t^/2 
would be written as 



/ Dk^exp{-a f d^xtrkl{x)}5{d''k^), (7.1) 



J Ja 

81 
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where a is the energy per unit length of the monopole current. Here, the (5-function 
is necessary to ensure the current conservation ^^k^J, = 0, and a is an ultraviolet 
cutoff larger than the screening length. To perform the path-integral (|7.1| ), we put 
the system on the 4-dimensional lattice with the lattice spacing a. In the lattice 
formalism, the monopole currents are defined on the dual lattice, 

kl^\s)^e/{An).a%{s). (7.2) 

The partition function is given as 

z = Eexp{-«'"*E^r(^)'}'^(^M^r(^))' (7-3) 

where a'"* = a/2 ■ (47r/ea)^. 

The lattice QCD simulation shows that one long monopole loop and many short 
loops appear in the confinement phase |^. Only the long loop becomes important 



for the properties of the QCD vacuum, while many short loops are originated from 
the fluctuation in the ultraviolet region and therefore can be neglected. In this 
system, the partition function can be approximated as the single monopole loop 
ensemble with the length L, 

Z = j:p{L)e~-^ (7.4) 

L 

where L and p{L) are length of the monopole loop and its configuration number, 
respectively. The monopole current with length L is regarded as the L step self- 
avoiding random walk, where 2D — 1 = 7 direction is possible in each step in the 
D dimensional space-time. Therefore, p{L) is roughly estimated as {2D — 1)^ = 7^. 
Using this partition function, the expectation value of the monopole current length 
is found to be 

(L) ^^^p(L)Le-'^ = I {" - '"P^ - '»" '!" > !"P^ - ;' (7.5) 

Z "Y \ oo if a < ln(2Z:' - 1). ^ ^ 

When the energy a is larger than the entropy \n{2D — 1), the monopole loop length 
is finite. However, when a is smaller than the entropy, the monopole loop length 
becomes infinite, which corresponds to monopole condensation in the current repre- 



sentation ||2^. Here, the critical value on monopole condensation Oc — \n(2D — 1) ~ 
ln7 = 1.945, which corresponds to the "entropy" of the self-avoiding random walk. 
Such a transition is quite similar to the Kosterlitz-Thouless type transition in (1+2)- 



dimensional superconductor |2T|], where vortex condensation plays an important role 
to the transition. 

In performing the simulation, we construct the monopole current as the sum of 
plaquettes of the monopole current because of the current conservation condition, 
d^kn = 0. As the initial current configuration, we prepare a random monopole 
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Figure 7.1: The monopole current in the monopole- current system in R^ at a fixed 



time : (a) the "critical phase'' 
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Figure 7.2: The histograms of the monopole- current length L of each monopole 
cluster in the monopole-current system. The data at each a are taken from 40 
current configurations, (a) In the "critical phase", only short monopole loops ap- 
pear, (b) In the " monopole condensed phase" near the critical point, monopole 
currents become in a cluster, (c) In the "confinement phase", there appears one 
large monopole cluster in each current configuration. 
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u 




Figure 7.3: The monopole density pmc as the function of a. The monopole density 
becomes small for large a and almost vanishes for a> ac = ln7 = 1.945. 



current system (hot start) or no monopole current system (cold start). Then, we 
update the link of the monopole current using the Metropolis method ||42|| . 

We generate the monopole current system using the monopole current action 
7r3| . Fig. [7.1| shows the monopole current on 8^ lattices for the typical cases, (a = 
1.5, 1.8, 1.9) at a fixed time. For these cases, we show in Fig. |7.2| the histograms of 
monopole loop length [^ of each monopole cluster with 40 current configurations. 
For large a, only small loops appear and monopole density is small. On the other 
hand, for small a, there appears a global network of one large monopole cluster, and 
the monopole currents are complicated and dense. The lattice monopole density 



Pmc 






^m(^) 



(7.6) 



and the clustering parameter 



v = 



{J:^Li 



63|| are shown in Fig. |7.3| . and Fig. |7.4| , respectively. Here, 
the i-th monopole cluster. For the extreme limit of r^ = 1 



(7.7) 

Li is the loop length of 
, all the monopole loops 
are linked in a cluster, while many monopole loops are isolated for small rj. As a 
decreases, the monopole density pmc becomes larger gradually for a ~ ac- However, 
the clustering parameter is drastically changed at a = 1.8 closed to ln7. Thus, these 
monopole current simulations clearly show the Kosterlitz-Thouless-type transition 
around ac = ln7 in agreement with the theoretical consideration on monopole (- 
current) condensation. 
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Figure 7.4: The clustering parameter i] = J2i Lf/{J2i LiY as the function of a. Near 
the critical point Oc = ln7 = 1.945, r] is drastically changed from unity to zero. 

7.2 Role of Monopoles for Confinement 



In this section, we study how these monopole currents contribute to the color con- 
finement properties. Quark confinement is characterized by the linear inter-quark 
potential, which can be obtained from the area-law behavior of the Wilson loop, 
{W) = {Pexp{ie § Afj_dx^)). Therefore, it is desired to extract the gauge variable 
from the monopole current k^. We now derive the abelian gauge variable in stead 
of the non-abelian gauge field, because the lattice QCD results show that the color 
confinement phenomena can be discussed only with abelian part to some extent. 
However, in the presence of the magnetic monopoles, the ordinary abelian gauge 
field A^{x) inevitably includes the singularity as the Dirac string, which leads to 
some difficulties in the field theoretical treatment. Instead, the dual field formalism 
is much useful to describe the monopole current system, because the dual gauge 
field B^{x) can be introduced without the singularity for such a system with k^ j^ 
and j^ = 0. This is the dual version of the ordinary gauge theory with Af^{x) for 
the QED system with j^ 7^ and /c^ = 0. 

Now, we apply the dual field formalism to the monopole current system discussed 
in Section [rT| . Since the monopole current k'-^^{s) is generated on the lattice with 
the mesh a. 



the continuous dual field B^{x) 



is derived from kfj_{x] 



using Eq.( |6.13 



in principle. In estimating the integral in Eq. ( |6.13D numerically, we use a fine lattice 
with a small mesh c. To extract B^{x) correctly, the mesh c is to be taken enough 
small. However, too fine mesh is not necessary because the original current k^!^^{s) 
includes the error in the order of a. Numerical analyses show that the use of c ~ a 
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Figure 7.5: The expectation value of the Wilson loop {W{I x J)) on the fine mesh 
c = a/2 for a = 1.7, 1.8, 1.9. The area-law behavior of (W) indicates the linear 
inter-quark potential, although the string tension is zero for a = 1.9 > etc- 



is too crude for the correct estimation of the integral in Eq. (|6.13|) . Instead, the 
calculation with c < a/2 is good enough for the estimation of B^{x), and hence we 
take c = a/2 hereafter. 

The expectation value of the Wilson loop (W) is shown in Fig. |7.5| . The Wilson 
loop exhibits the area-law and the linear confinement potential: In(iy) decreases 
linearly with the inter-quark distance, where its slope corresponds to the string 
tension. Quantitatively, the string tension is measured by the Creutz ratio, and we 
show in Fig.^]6| x(3, 3) as a typical example. For the monopole condensed phase as 
a < ac, the string tension gets a finite value, while it vanishes for the non-condensed 
phase of monopole as a < ac- Thus, the confinement phase directly corresponds to 
the monopole condensed phase and therefore monopole condensation is essential as 
essence of the confinement mechanism. 



7.3 Monopole Condensation in the QCD Vacuum 

In this section, we compare the lattice QCD with the monopole-current system in 
terms of monopole condensation and confinement properties, as shown in Fig. |7.7| . 
The lattice QCD simulation shows that the QCD vacuum in the MA gauge holds 
the global network of monopole current as shown in Section |4.3| . Considering the 
similarity on the monopole clustering, the QCD vacuum can be regarded as the 
monopole-condensed phase with a < a^ in the monopole current system, as shown 
in Fig.lO(c). Such identification of the QCD vacuum with the monopole-condensed 
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Figure 7.6: The Creutz ratio as the function of a in the monopole current system. 
The dotted hue denotes the Creutz ratio as the function of /? = 1.25a in the lattice 
QCD. 



phase is also suggested in terms of the confinement properties, because the con- 
finement phase corresponds to the monopole condensed phase as shown in Section 



7.2. 



7.4 Comparison with Vortex Condensation in 1+2 
Superconductor 

In this section, we summarize in Table 7.1 the correspondence between monopole 
condensation in the QCD vacuum and vortex condensation in the 1+2 dimensional 
superconductor |21[] in terms of the topological object and the phase transition. For 
these two systems, condensation of the line-like topological object plays the relevant 
role for the determination of the phase. The phase structure is controlled by the 
balance of power between "entropy" (configuration number of the topological object) 
and "energy" (self-energy of the topological object). 

In the 1+2 dimensional superconductor, the Abrikosov vortex is the important 
topological object. While the vortex scarcely appears at the low T, such topological 
excitations frequently occur at the high T, which enhances the entropy factor in 
the free energy. Then, at the critical temperature, vortex condensation occurs due 
to "entropy dominance" and the system goes to the normal phase. This phase 



Figure 7.7: Comparison between QCD in the MA gauge and the monopole current 
system. The QCD vacuum corresponds to the monopole-condensed phase in the 
monopole current system because of the similarity on confinement and the monopole 
clustering. 

transition is known as the Kosterlitz-Thouless transition. 

Similarly in the QCD vacuum, the monopole current plays the relevant role as 
the line-like topological object in the 4-dimensional space-time. In the monopole- 
current system, the control parameter is the self-energy a of the monopole current. 
While only the local fiuctuation of the monopole current appears for a > «£ ~ ln7, 
monopole condensation occurs for a < ac as the result of entropy dominance on the 
monopole-current configuration. Moreover, this monopole condensation leads to the 
electric confinement. 



System 


1+2 dim 
Superconductor 


4-dim 
Monopole system 


Line-like 
topological object 


Vortex 


Monopole current 


Control parameter 


Temperature T 


Monopole self-energy a 


Condensed phase 
(Entropy dominance) 


High temperature 
Normal phase 

T 

Vortex condensation 


Small a 
Confinement phase 

t 
Monopole condensation 


Non-condensed phase 
(Energy dominance) 


Low temperature 
Superconducting phase 


Large a 
Deconfinement phase 
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Table 7.1 Correspondence between the monopole current system and the 1+2 dimensional 
superconductor 

7.5 Monopole Size and Critical Scale in QCD 

Up to now, we have argued about the correspondence between the lattice QCD and 
the monopole current system described by { \l.l\j in terms of the qualitative aspects as 
monopole condensation and confinement. In the final section of this chapter, we at- 
tempt further consideration on the correspondence between the coupling constants, 
a and /3 = ^^. Here, we note also the limitation of the simple "monopole current 
approach" in the ultraviolet region and the critical scale of the dual superconductor 
picture of QCD. 

To begin with, let us consider one magnetic monopole with an intrinsic radius 
R in the multi-monopole system. In the static frame of the monopole, it creates a 
spherical magnetic field, H(r) = |^r = ^^ ioi r > R and H(r) = |J^r = j^ 
for r < R, where the QCD running gauge coupling e(r) is used to include the vacuum 
polarization effect. Here, the effective magnetic-charge distribution is assumed to 
be constant inside the monopole for the simplicity. 

Now, we consider the lattice formalism with a large mesh a > R. The electro- 
magnetic energy observed on the lattice around a monopole is roughly estimated 
as 

M{a) ~ r d^xluir)' ~ ^ = -^, (7.8) 



a 



2 Sna e'^{a)a' 



which is largely changed depending on the lattice mesh a. This simple estimation 
neglects the possible reduction of g{r) in the infrared region due to the asymptotic 
freedom of QCD. The screening effect of the magnetic field by other monopoles also 
reduces g{r) effectively in a dense monopole system. However, M{a) is modified by 
at most factor 2 {M{a) = -^rnx^) even for the screening case as g{r) = g{a)-6{2a — r). 
Then, even in the multi-monopole system, M{a) would provide an approximate value 
for the electromagnetic energy on lattices created by one monopole, and we call M{a) 
as "lattice monopole mass". 

For the large mesh a > R, the monopole contribution to the lattice action reads 
S = M[a)a-L, where L denotes the length of the monopole current measured in the 
lattice unit a. Therefore, M{a) is closely related to the monopole- current coupling 



a 



and (3 = 2Nc/e , a ~ M{a)a ~ -Jf^ = | ■ /3su{2)- For the above screening 
case, this relation becomes a'"* ~ M{a)a ~ -^^ = | ■ /?su(2), which is consistent 
with the numerical result, a = 0.8/?, discussed in Section |7]^. Here, as long as the 
mesh is large as a > R, the lattice monopole action would not need modification 
by the monopole size effect, and the current coupling a'"* is proportional to /?. 
Quantum mechanically, there is the energy fiuctuation about a~^ at the scale a, and 
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therefore monopole excitation occurs very often at the long-distance scale satisfying 
M{a) ~ a~^. Thus, M{a)a ~ a'°* is the control parameter for monopole excitation 
at the scale a(> R), and we can obtain the quantitative criterion for "monopole 
condensation" as M{a)a ~ \n{2D — 1) from the analysis using the current dynamics 



in Section 7.1 



Second, we discuss the ultraviolet region with a < R. In the current dynamics, 
there exists a critical coupling ac — In (2D — 1) as shown in Section [7.1| and |7.2| . 
Above ac, the lattice current action provides no monopole condensation and no con- 
finement, while /3 — > cxD can be taken in the original QCD keeping the confinement 



property shown in Fig. |7.6| . Such a discrepancy between P and a can be naturally in- 
terpreted by introducing the monopole size effect. Obviously, the monopole- current 
theory should be drastically changed in the ultraviolet region as a < R, if the QCD- 
monopole has its peculiar size R. 

Here, let us reconsider the relation between a and a'"* in the lattice current 
action. Similarly in the lattice QCD, the action has no definite scale except for the 
lattice mesh a(> R), and therefore the scale unit a would be determined so as to 
reproduce a suitable dimensional variable, e.g. the string tension a ~ IGeV/fm, 
in the monopole current dynamics. For instance, a is determined as a function 
of a'"* using the Creutz ratio x — (^(^^ i^^ Fig. |7.6| . Therefore, a should reach R 
before realizing a'"* — > a^*^*, and the framework of the current theory is to be largely 
modified due to the monopole size effect for a < R. 

In conclusion, the QCD-monopole size R provides a critical scale Re for the 
description of QCD in terms of the dual Higgs mechanism. Quantitatively, the 
difference between the monopole current system and QCD appear larger than (3 = 
1.25a, which corresponds to Re — 0.25fm. In the infrared region as a > Re, QCD 
can be approximated as a local monopole- current action PD[ , and the QCD vacuum 



can be regarded as the dual superconductor. On the other hand, in the ultraviolet 
region as a < Re, the monopole theory becomes nonlocal and complicated due to 
the monopole size effect, and the perturbative QCD would be applicable instead. 



Chapter 8 

Dual Ginzburg-Landau Theory 



The lattice QCD simulation shows abelian dominance and monopole dominance for 
the string tension in the maximally abelian (MA) gauge p3| , p^ p6| , and hence the 
monopole is considered to be the essential degrees of freedom for the confinement 
properties. In the confinement phase in the MA gauge, there appears the long 
and complicated monopole current covering the whole lattice space, which suggests 
monopole condensation. Thus, monopole condensation scheme is one of the realistic 
candidates of the physical interpretation for the confinement mechanism. In this 
chapter, we study the confinement phenomena using the dual Ginzburg-Landau 



(DGL) theory ||22|, ^, where the confinement mechanism is described by monopole 



condensation [|1^, |1^, [T^ . 

The DGL theory [^2|, ^] is the infrared effective theory of QCD based on the dual 



superconductor picture. In the ordinary Ginzburg-Landau theory, condensation of 
the Cooper-pair field leads to exclusion or squeezing of the magnetic field in the 
superconductor. Here, the U(l)e-gauge symmetry is spontaneously broken and the 
abelian gauge field, the photon, becomes massive as the result of the Higgs mech- 
anism. The Cooper-pair field plays a role of the Higgs field. In the QCD vacuum, 
condensation of the "monopole field" would provide the exclusion or squeezing of 
the color-electric field through the dual Meissner effect. In the monopole-condensed 
system, the dual gauge symmetry is spontaneously broken and the "dual gauge field" 
becomes massive. 

Let us formulate the DGL theory starting from QCD. The QCD Lagrangian has 
the SU(A^c) symmetry and is written by the quark field q{x) and the gluon field 
Af,{x) as 

£qcd = -^tr {G,Mn +q{il^D^^ - m) q, (8.1) 

where G^i. denotes the SU(A^c) field strength, 

G^, = ^{[D^,D,]-[d^,d,]) with D^ = d^ + ieA^. (8.2) 
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In terms of the Cartan decomposition, the gluon field A^{x) = y4"(x)T" G su(A^c) 
with Nc=3 is decomposed into the diagonal component A^ = (A^, A^) and the 
off-diagonal component C^ G C, 



A, = A,. H + J2iC;'^E, + C;E.,). 

a=l 

Here, the Cartan subalgebra H is defined as 

and the raising and the lowering operators E±a{a = 1,2,3) are defined as 



(8.3) 



U) 



L5) 



H,Ea 



daEa and 



H,E_a 



-daE_a, where «„ denotes 



One finds the relations 

the root vector of the SU(3) algebra: ai = (1, 0), 02 = (— |, ^ 2 ) ^^'^ "^3 ^ (~h 2 
The off-diagonal gluon component C^ is written as 

Cl = ^{Al + zAl), Cl = ^{Al-zAl), Cl = ^{Al + zAl). (8.6) 



According to the Cartan decomposition, the QCD Lagrangian is expressed as 

^/'Abcl I /'off 
QCD — ''- + i- , 

where £^^^1 is diagonal part of Cqcb, 

Z:^^^' = -^{d,X-d,A,y+q{ij,d^-ej^A^-H-m)q, 
and C"^ is remaining part including the off-diagonal gluons. 



^.7) 



C 



off 



le 



E 

a=l 
3 



{{d, - ^6(5-^ ■ 4)) A C^)^^ - ^a,tcC;'C, 






-iY.iU^-^a){c*''^cT''+- 



3 



^a,(C*'^AC"; 



^v 



.a=\ 



-e E {CT^^l^^aq) + C^"(g7^i5-ag) 



a=l 



L9) 



Pl^V 



with (a Ab)f,^ = a^h^ - 6 

In the abelian gauge defined by diagonalizing a gauge dependent operator $, the 
SU(3) gauge theory is reduced into [ U(l)3 x U(l)8]'^ gluon gauge theory p2|, |3^ 



and monopoles appear. Particularly, in the MA gauge, off-diagonal components are 
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forced to be small and behave as massive fields, while diagonal gluon components 



remain largely propagate over the long range |^. In fact, in the MA gauge, off- 
diagonal gluons are infrared-irrelevant and can be neglected at the long distance. 
Therefore, after MA gauge fixing with Q, we construct the abelian projected QCD 
by dropping off the charged part and define the abelian gauge field Afj, as 

A, = ^tT{A,H)-H = A^-H. (8.10) 

Accordingly, the SU(3) field strength 

G% = {d^A,~d,A,)+ze[A^,A,] + '-Q[d„d,p (8.11) 

is modified as 



where JP^™^ = -VL[d^^dy\VL^ . This corresponds to Eg. ( [4. 241) in the lattice formalism. 



Here, the breaking of the abelian Bianchi identity is brought by the appearance of 
T^uu^- In the SU(3) case, there are three kinds of magnetic charges of the QCD- 
monopole corresponding to the subspace on degeneracy of the diagonalized operator 
$diag; for ^diag = diag($i, $2, $3), $1 = $2, '^2 = -^s and $3 = $1. For example, 
the magnetic charge is proportional to the root vector cJi in the case of the $1 = $2- 
Thus, the magnetic charge of the monopole current is discretized as 

3 

fh = g'^^aOa, ^a e Z. (8.13) 

a=l 

For A^ monopole system, the magnetic current k^ = k/j^ ■ H is expressed as 

k^{x) = -gY.m dTi^^5\x-xi{Ti)). (8.14) 

In the DGL theory, the monopole field is introduced by integrating over all trajec- 
tories of the monopole particle as will be shown in the later section. 



8.1 Dual Gauge Field 

The extended Maxwell equation is described by the field strength F^y as 

d^F^i^=f, d;F''^' = k^ (8.15) 

with the electric current j^ and the magnetic current k^, and the electro- magnetic 
duality {F^^ ^^ *Ffj_y, j^ <->■ fc^) is manifest in this formalism. In the ordinary QED, 
the description by the gauge field A^ is useful, since QED includes only the electric 
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current. In this description by A^, however, this electro-magnetic duality in the 
Maxwell equation is not manifest and the introduction of the magnetic monopole 
leads to the singularity in the gauge field A^. As shown in Section |673| , for "the dual 
system" including only the magnetic current, it is rather useful to take the dual 
description with the dual gauge field B^ instead of ^4^. Here, the dual gauge field 
B^ is introduced so as to satisfy 

*F^'' = ^e^"'''^F^p = d^'B'' -d^B^", (8.16) 

or equivalently, the electro-magnetic field is written as 

W = -*F°* = -a°5* + a^fi°. (8.18) 

Now, we adopt here the Zwanziger formalism [Q described both with A^ and 



Bfj^ in order to keep the electro-magnetic duality in the extended Maxwell equation 
manifest. In the general case with j^ ^ and /c^ 7^ 0, we rewrite the field strength 
to keep the duality manifest in the following way. For an arbitrary anti-symmetric 
tensor G^i, and an arbitrary constant four- vector n^, there is an identity 

G = ^{[nA{n- G)] -* [n A {n ■* G)]} (8.19) 

with (n ■ G)^ = njG'^^. Substituting G = F = (5 A A) and *G = *F = (9 A B), we 
get 

F = -^iinAln-idAA)]}-* {nA\n-{dAB)]}), (8.20) 

*F = ^{*{nA[n-{dAA)]} + {nA\n-{dAB)]}), (8.21) 

and then the extended Maxwell equations d^F'^^^ = j^ and d*F'^^ = k^ are written 
as 

d^F''^ = ^(^[{n-dyA^'-{n-d)d^'{n-A)-n^{n-d){d-A)+n^'d\n-A)} 

-{n-d)e>',^^n^d-B^)=j^, 
dlF""^ = ^[{{n-dyB^'-{n-d)d^'{n-B)-n^'{n-d){d-B)+n^'d\n-B)} 

+ {n ■ d)e^'^^xn''d^A^) = k^". 
Therefore, the Lagrangian which leads to these field equations is derived as 

C = -^[n-{dAA)][n-*{dAB)] + ^[n-{dAB)][n-*idAA)] 

^ [n-{dAA)f-^[n-{dAB)f~j^A'^-k^B'^. (8.22) 



2n2 ' ' '' 2n2 
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In the Zwanziger formahsm, the duahty transformation F^^ ^^ *F^u is given by 
the replacement A^ ^->- B^ and the electro- magnetic duality is manifest, and the 
system can be described without singularities. Particularly in the absence of j^ and 
kfj,, the theory can be written using "only A^" , "only B^" , or "both A^ and B^" in 
the Zwanziger formalism as 



/'pv4e^/'^'^^''l^l = fvBe'-f'^'''^°^^^= j VAVBe'^ '^''"^^^^'^^ 



where 



CM] 
Co[A,B] 

cm 



-\{^^A)\ 



^.23) 



(8.24) 



-^^[n.{d ^A)][n.* {d ^B)] + ^^[n.{^ ^B)][n.* {d ^A)] 



^ 'n-{dAA)f 



2n2 
1 



^[n-(9A5)f 



4 



{^^BY. 



(8.25) 
(8.26) 



In the Zwanziger formalism, the Lagrangian can be written with holding the du- 
ality between the gauge field A^ and the dual gauge field B^. However, occasionally 
it is more useful to describe only with A^ or B^ for practical calculations. Here, 
we try to describe the Lagrangian with only A^ or B^ by the functional integration 
over B^ or A^ in the partition functional |^ . 



Since the dual gauge field B^ is bilinear in the Lagrangian ( |8.25| ), B^ can be 
analytically integrated in the partition functional in the Zwanziger formalism. 



fvAVBe'^'^'^^^'^'^^ = fvAe'-f'^'''^''^^^ 



(8.27) 



with 



Ck[A] 



-{{dAA) 



/.w 



psingia 
IJ,u J • 



f8.28) 



Here, the singular term -F^™^ = £^iuai3:^n^k^ appears due to existence of the mag- 



netic current fc^, and includes the nonlocal operator (a;| ;^ | 



{,Xn-yn)r{,x^-y^). 



Thus, the field strength F^^ includes the nonlocality in this description only by 
A^j, in the presence of the magnetic current k^. Similarly, the theory can be described 
only with the dual gauge field B^ after the functional integration on A^. In this 
case, the theory is replaced by {B^^j^) in the above argument, and the Lagrangian 
reads 



CAB] 



-}^{*{dAB),,-e,..p^n^3^}\ 



(8.29) 



where the nonlocal part appears *F^y reflecting existence of the electric current j^. 
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In this way, the system including both the electric and the monopole currents 
can be described without any singularities using the Zwanziger formalism, where 
the dual field B^ is introduced as well as the gauge field A^. On the other hand, the 
description only with A^ or B^ leads to the appearance of the nonlocal operator in 
the field strength in the general case. In constructing the DGL theory from abelian 
projected QCD, the above formalism can be applied by the simple replacement of 
Af,-^ Af, = Af,- H, B^^ B^ = B^- H and so on. 



8.2 Monopole Field from the Monopole Particle 

In this section, we derive the Lagrangian for the monopole field from the monopole- 
particle trajectories following Stone and Thomas |^ ^ ^. At the infrared scale 
in the MA gauge, the monopole current becomes dense and complicated, and the 
monopole currents interact each other only in the short distance due to the screening 
effect ||5^ corresponding to monopole current. The infrared effective action fc^ is 
expected to be described by the local action of /c^ 

rffc^e""^'. (8.30) 

Here, a is the monopole self-energy and the "mesh" a is introduced as an ultravi- 



olet cut-off larger than the screening length |^. As shown Section |4.1| , almost all 
monopoles have the unit charge in the abelian gauge, and therefore the partition 
functional is approximately described by the ensemble of the monopole loops 

^N CO 

Z = Y.^ with Zioop ^ ^ p(L)e-"^ (8.31) 

where p{L) is the number of closed loop configurations with the length L. Here, we 
have taken the lattice unit a = 1. Each monopole loop with the length L can be 
approximated as L step random walk |6^ because of the absence of the nonlocal 



interaction. 

Let us consider the grand canonical ensemble for closed-loops on the D-dimensional 
lattice in the unit of the lattice spacing 1. At the infrared scale, the partition func- 
tion for a single monopole loop is written as 

Zioop = / dTp{r)e-''\ (8.32) 

Jo 

The probability distribution of starting at x and ending at x' after a walk of r steps 
is given as 

p{x,x',t) = / (i[x(-)]exp{— / -— x(r')^(ir'}, (8.33) 

Jx Jo 2 



98 8 Dual Ginzburg-Landau Theory 



which satisfies 

dp{x,x',T) 1 d'^p{x,x',T) 



dr 2D dxl 



(8.34) 



The total number of paths from x to x' with r steps is estimated as 

r(x,x',r) = {2DYp{x,x\t) 

= N [ d[x{-)]exp{- r dT'—x'^-\n{2D)}, (8.35) 

J Jo 2 

where the factor 2D is the configuration number for a walk of one step. Using the 
eigen-function ^rn{x) satisfying 

{-^dl-\n{2D)}^U^)=LU^^^ and Jd''x^lix) = l, (8.36) 

r(x, x', t) is written by 

r(x, x', r)=Y, <l>„(x)<l>„(a;')e--"^ (8.37) 

n 

For the closed loop, all r starting points on the loop define the same configuration, 
and hence r(x, x',r) satisfies 

I d^xT{x, X, r) = tp{t). (8.38) 

Then, we get 

- ^§^ = r drrpir)e-''^ = T dr [ d^'xTix, x, r)e-"^ = T drV] e-(-"+")^ 
oa Jo Jo J Jo n 

= E^— = Tr(-^52-ln(2D) + a)'\ (8.39) 

„ cj„ + a \ 2D ^ J 

Integrating on a, we obtain 

Zioop = -Trln(-^9j-ln(2D) + a) 

= -lnDet(-^9j-ln(2Z}) + a). (8.40) 

Thus, the partition function for the whole system is expressed as 

7^ 1 

Z = E^ = exp[^ioop] = Det-^(-^5; + m^) 

= J V^ J Vv*exp[- J d^x^*{-^dl + m')^]. (8.41) 

Hence, noninteracting loops can be written as a free-field theory of the complex 
scalar field with the mass m^ = a — ln(2D) 
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Here, considering the probability weight of e ""^5 the "expectation value" of the 
number of the walk which starts at x and ends at x' is proportional to 



dTT{x,x' , r)e 



^ $n(x)$„(x') 1 I \/ I /\ 

= > => {x\n)(n\ — i— — -\n){n\x) 

= ^ J V^ J V^*ip{x)ip*{x')exp{~ J dxip*{-^dl + m')ip} 

= {ipix)ip*ix')). (8.42) 

In particular for \x — x'\ -^ 00, {{p(x)ip*{x')) indicates the density of loops whose 
length is infinity and becomes {{p){ip*) = Kv^)^. Thus, the expectation value \{(p)\ 
of the free massive scalar field is directly related to the loop density. When the self- 
energy a is larger than "entropy" ln(2D), i.e., m^ > 0, the loop density is suppressed 
by the action factor and the system is described by the free scalar field theory. On 
the other hand, when a is less than ln(2D), i.e., m^ < 0, infinite long loop appears 
and the expectation value of the scalar field takes a non-zero value, i.e, the scalar 
field (f condenses in the description of the scalar field theory. In order to stop 
the loop density being infinite, we need — A|v9|^ with positive A in the Lagrangian, 
corresponding to a short-distance repulsion in Eq.( ^.4lD . 



In the monopole-current system (|8.3CI| ) appearing in the QCD vacuum, this com- 
plex scalar field (p corresponds to the "monopole field". Here, the monopole current 
system can be regarded as a self-avoiding random walk, and then the factor "2Z)" in 
this argument is expected to be replaced by "2D — 1" . Hence, whether the monopole 
condenses or not is determined by the balance of the "self energy" a and the "en- 
tropy" \n{2D — 1), as is already shown in the Chapter |^. In the presence of the dual 
gauge field i?^, the derivative 9^ is simply replaced by the dual covariant derivative 
dfj, + igB^ with the dual gauge coupling constant g. Finally, the Lagrangian for the 
monopole field x — T^V i^ written as 

with 2At;2 = 2Dm^. 

In the SU(3) QCD in the abelian gauge, there appear three kinds of monopoles 
with the unit charge m^ = icJa (a = 1, 2, 3) corresponding to the nontrivial root 
vector [^. Among the three monopole currents fcJJ (a = 1, 2, 3), only two monopole- 
current degrees of freedom are independent due to the relation Yl\=i <3a = 0. In fact, 
/c^ can be expressed as — (fc^ + fc^)- However, the description by the two independent 
monopole currents makes the global Weyl symmetry unclear and becomes compli- 
cated in the presence of the interaction with B^ = B^-H. Therefore, in order to treat 
A;^ (a = 1, 2, 3) as the independent variable, we introduce the Lagrange-multiplier 



100 8 Dual Ginzburg-Landau Theory 



field B^ and rewrite the sum over the monopole current as 

E- E kJ:K)= E /l)i?Je^/^^^^^"S-^^> (8.44) 

:?-,.fc2 kl.kl.kl a=l kl.kl.kl 



Jul h.2 ul ^2 l3 rt = l ^1 ^2 l3 



in the partition function of the monopole current. This prescription corresponds to 
the introduction of the unphysical U(l)™ magnetic charge coupled with 5°. Now, 
we can apply the above formalism for each /c" (a = 1, 2, 3) independently, and then 
we obtain the partition functional, 

a=l 

r ^ 
•exp{z / d^x J2{\i9, + ^gaa ■ B, + 2(?5j)xa|' - AdxaP - t;2)'}},(8.45) 

■^ a=l 

where «„ denotes the non-trivial root vector; ai = (1,0), 02 = (~2'~2) ^^"^ "^3 ~ 
(— |, ^). Since the partition functional Z is invariant under the transformation 

^^^^) ^ e^'(^^Xa{x), Bl{x) ^ Bl{x) + -9^^(x), (8.46) 

the theory holds the extra dual gauge symmetry U(l)™ as well asU(l)^ x U(l)™. 
Here, the global Weyl symmetry is kept manifest in Z. Due to the extra local U(l)[}* 
symmetry, the overall phase degrees of freedom of Xa{x) in Eq.( ^.45| ) is absorbed 
into B^ and becomes unphysical. In fact, without loss of generality, we can set the 
local condition as 

3 

EargXa = 0, (8.47) 

a=l 

by using a suitable U(l)[}^ dual gauge transformation and the shift of the integral 
variable B^ . In principle, B^ can be integrated out in Z. After the integration on 
5° there appear some interaction terms as the functional determinant, however, 
those contribution is assumed to be included into the self- interaction term of Xa- 
Thus, the final expression of Z is obtained as 

Z = /'DS^e^/'^'^^-^(^^^)'^[n Dxa] exp{i f d'^xC"'"''}, (8.48) 

a=l 

where £™°° denotes the monopole part of the Lagrangian in the SU(3) case, 

^mon ^ J2{\{d, + igaa ' B ^)Xa? " K\Xa? ' v'f}, (8.49) 

a=l 
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with the constraint J2 argx^ = 0. 

a=l 



In the Zwanziger form (^.25), we obtain the dual Ginzburg-Landau Lagrangian 



as 



£ 



DGL 



2^2 
1 



[n-idAA)Y[n-*{dAB)l 



2^2 



[n-{dAB)Y[n-*{dAA)], 



2^2 

3 r 

+ E 

a=l 



n-{dAA)f 



2^2 



[n-{dA B)Y+q{i-f^d'' - e-f^A" ■ H - m)q 



\{id^- gda- B^)Xa -X{\Xaf-v'^y 



^.50) 



with J2 argXa = 0. 

a=l 

In terms of the Cartan decomposition, this DGL Lagrangian is also expressed as 
Cdgl= tr{ -hn- {d AAY[n-*{d AB)l + hn- {d AB)Y[n-*{d AA)l 
-\[n ■ {d A A)f -hn-{dA B)f}+q{i^^d>' - e^^A^ - m)q 

^z ^z 



+ [DT\xnDt''''\x]-Kx^X-vY}, 



(8.51) 



where A^, S^, Z)^""^ and x denote ^^ = A^-H, B^ = B^^H, D;1"'^i = d^ + igB^ and 
X = I]a=i ^XaEai respectively. The kinetic term of the monopole field leads to the 
original expression in Eq. (|8.5CI| ) 

ir{[bf-\xnDf-\x]) 

= tT{[d^ + igB^,x^][d^, + igBf„x\) 

= 2ti{{d^x:E-a + igB^xllH, E_a\)id^XbEt + tgB^Xb[H, E,])} 

= (d^ + igda ■ B^)xa ' . (8.52) 



Here, the DGL Lagrangian ( ^.50|) has [U(l)g]2 = U(l)g x U(l)g gauge symmetry and 
[U(l)^]2 = U(l)^ X U(l)^ dual gauge symmetry, since the Lagrangian is invariant 



under the gauge transformation by Vt^ = e 



-iOe 



-iOe-H 



G U(l)^ X U(l)^ 



q -> Qe{x)q, Af^ -^ Qeix) [A^ df, ) Qeix^ = Af, + -df,e, 



(8.53) 



and also invariant under the dual gauge transformation by Qm = e ■""" = e 
U(l)^ X U(l)i 



B,. 



X 



Qmix) i B^ d^ j Qmi^y = B^ + -d^9m, 

^^mX^^m / , Xa^ ^a- 



(8.54) 
(8.55) 
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In terms of the Weyl symmetry, it is useful to rewrite the field variable 

'A^-(^i \ /A^ 



A^ = A^-H=\ A^-uJ2 \ = \ S Al |, (8.56) 
l.-cDs, 



V A^ 



where cJa(a =1,2,3) denotes the SU(3) weight vector. Here, A^, A^ and A'^ satisfies 
the relations, 

E ^m(^) = 0' (8-57) 

a&{R,B,G} 

{Ay + {Alf + {A^f = l{{Aiy + {Ain (8.58) 



and 



{dAA)l = {dAAX^+{dAA')l 



= 2{(aAA% + (9AA% + (9AA%}. (8.59) 

Similarly, the dual gauge field B^j, and the quark field q are written as 
(B^ \ (q^\ 

B,^\ Bl , g = g^ . (8.60) 

V B^J [q^J 



Here, -B" also satisfies the similar relations to Eqs. (|8.57| )- (|8.59| ). The DGL La- 
grangian in the gauge sector is expressed as 

/'gluon \ ^ /'gluon 

a£{R,B,G} 

/:f°" = -^[n-{dAA)X[n-*{dAB,)l + ^[n-{dAB)r[n-*{dAA,)l 
_ J_[^ . (5 A A,)]' -L[n.{dA B,)]'. (8.61) 

The quark sector in the DGL theory is written as 

ae{-R,B,G} 

In this expression, the U(1)rX U(1)bX U(1)g symmetry seems to hold, however, 
one of U(l) symmetries is fixed by the constraint as ( ^.57|) . Under the Weyl trans- 



formation as A^ ^-> Aj^, B^^ <-^ B^^ and g^ <-*> qt with a, b & {R,B,G}, the Lagrangian 
J2a&{R,B,G} l^a IS manifestly invariant. Since the Lagrangian including the monopole 
field X is also invariant under the Weyl transformation as x ~^ x' = WxW^ with 
W G Weyl, the DGL Lagrangian has the Weyl global symmetry as the relic of the 
SU(3) symmetry. In this way, the dual Ginzburg-Landau theory has the [U(l)3 x 
U(i)g]|^oca/ ^ [u(l)3 X U(1)8]J^'="' X Weyl symmetry. 
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8.3 Dual Meissner Effect by Monopole Conden- 
sation 

Based on the dual superconductor picture, monopole condensation leads to the color 
confinement, which is brought as the result that the color-electric field is excluded 
from the QCD vacuum through dual Meissner effect. In this section, we show the 
dual Meissner effect caused by monopole condensation in term of the DGL theory. 
We separate the monopole field Xa into the mean field {xa) = v and its fluctuation 

Xa = {v + Xa)e'^'', (8.63) 

where the monopole condensate does not depend on a because of the Weyl symmetry. 
Due to the constraint Yl\=i argXa = 0, there are only two independent degrees of 



freedom among the three phase variables ^a {ct =1,2,3). Using this expression ( ^.631) 
the Lagrangian (|8.50|) is expressed as 



-DGL 



^ [n-{dA A)Y[n ■ *(9 A B)], + -l^[n ■ {d A fi)]^[n • *(9 A A)], 



2^2 L V n I y J,^ 2n2' 

^ [n- {d A A)f - ^{n- {d AB)f 



2n^' ' '' 2^2' 



+qit^,d'^ - e^,A^ ■H-m)q + ]-mlB^ + ^ [{d.Xaf - -'^~' 






+ Y. [/(«a • 4)'(Xa + '^Xavf - X{4vxl + xt)] , (8-64) 



a=l 



where me = V^gv and m^ = 2^/Xv denote masses of the dual gauge field S^ and 
the monopole Xa, respectively. Thus, monopole condensation makes the dual gauge 
field massive and the two independent phase variables of the monopole field are 
changed into the longitudinal degrees of freedom of the dual gauge field, which is 
the dual Higgs mechanism. Here, the color electric field E cannot propagate over 
the longer distance than l/rriB and the dual Meissner effect occurs. In this way, the 
color-electric flux is confined into the QCD vacuum by the similar mechanism on 
the magnetic flux in the superconductivity. 

8.4 Color-flux-tube in the Monopole-Condensed 
Vacuum 

As a result of the dual Meissner effect, the color electric flux is squeezed like a tube in 
the monopole condensed vacuum. In this section, we now investigate the structure 
of flux-tube with a quark and an antiquark at the both ends. In the standard 



Figure 8.1: The electric charge of quarks and the magnetic charge of monopoles. 
Since the quark color Qi is vertical to the magnetic charge gai of the monopole field 
Xi, condensation of xi does not contribute to confinement of the color Qi. 



notation, [^, ^, the quark charges are Qa = eWa, where Wa{a = 
weight vectors of the SU(3) algebra, wi = (0, —-j£), W2 = (— |, -^h^ 
for the three color states 



= 1, 2, 3) are the 
75;, W2 = l-f, ^), W3 = (|, ^), 
, red(i?), blue(-B) and green(G), respectively. Using 
the Gauss law, one finds the color electric field E and then the dual gauge field B^, 
obeying F^,^ = *{d A B)^^, are proportional to the color charge Q. The monopole 
Xa couples with the quark charge Qb in the form of Oq • QbXa- We note the algebraic 
relation between the root vector and the weight vector as. 



2wa ■ab = J2 ^<^hc e {-1, 0, 1}, 

c=l 



(8.65) 



and therefore one kind of monopole Xa couples with not three but two of the quark 
charges as is shown in Fig. ^.l| . 

For the example, in the case of (R-R) system, |xi| is never affected; \xi\ = v. 
For this case, we can rewrite Eq.(|8.5Cl|) as. 



1 
3 



id,B^-d^By + 2 



^DGL - 

where B^ = WiB^ and |xi| = v, X2 = 
invariant under the transformation, X2 
lagrangian in the simple GL form; 



Rn,,R 



id, + -gB^)x 



2A(x^ -v'YiS.QQ) 



X"" and xs = X^*; because the system is 
In this case, we can rewrite the DGL 



,R 
X3- 



£ 



DGL 



-\id,B^ 



duB,y+ {d, + tgB,)x -Xi\x\-v') 



-,2\2 



(8.67) 
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Figure 8.2: The color-electric field E{r) the monopole condensate x{f) ^i^re plotted 
as the function of r in the cylindrical coordinate. 



where the field variables and coupling constants are redefined as 



S„ 






9 



V3 



g, x = V2x'', 



V2i 



A 



A 



f8.68) 



We get the same expression for the other two cases {B-B, G-G). Hereafter, we will 
drop the notation " since there is no confusion. The field equations for 5^ and x are 
derived by the extreme condition, 



d'x + 2tgB^{d,x) + ^9{d,B^)x - g'B^x + 2A(|xr - v'Yx = 



(8.69) 



d,{d>'B^ - d'^B^) + ig{{d''x*)x - {d''x)x'} + 2(?^5lxr = 



^.70) 



For these equations, there is a static solution of vortex, which is known as the 
Nielsen-Olesen vortex |7y]. Using the cylindrical coordinate(r, 6', z), we consider the 
static solution satisfying 



B = Beee, 



E=-_[r5(r)]e, 
r dr 



X = X{r)e 



ind 



f8.711 



under the suitable boundary condition. 

We show in Fig. |8.2| the profile of the color electric-field E{r) and the monopole 
field xi''^)- Apart from r = 0, the monopole condenses, while the monopole conden- 
sate disappears near the center of fiux-tube. Accordingly, the color electric-field is 
squeezed like a one-dimensional tube, whose radius is about l/niB- In this way, in 
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the QCD vacuum, there are three types of the flux-tubes according to the color of 
quarks, being different from the superconductivity. However, the DGL Lagrangian 
in the all color system becomes a simple GL type and there is a solution where the 
color-electric field is squeezed like a one-dimensional tube as the Abrikosov vortex. 
This result is consistent with the Regge trajectory of hadrons and the lattice QCD 
results, which indicates that the QCD vacuum can be regarded as the dual version 
of the superconductor. 

8.5 Quark Confinement Potential 

In the DGL theory, the color-electric flux between quarks is squeezed into one- 
dimensional tube as the result of monopole condensation. This indicates that the 
static potential between quark and anti-quark is linear confinement potential. Here, 
we investigate the interquark potential using the gluon propagator including the 
nonperturbative effect [Q. 

The static potential between heavy quarks can be obtained from the vacuum 
energy where the static quark and antiquark exist. Here, we take a quench approx- 
imation, i.e. we neglect the quantum effect of the hght quarks. As a first step, 
we approximate the monopole filed Xa as the mean field \xa\ = v neglecting the 
effect of the quark sources on the condensed monopole. Later, we include the cor- 
rection in order to eliminate the divergence originated from this approximation. In 
this scheme, information on confinement is included in the gluon propagator, which 
leads to the strong interaction in the infrared region. 

The vacuum energy V{j) in the presence of the static quark sources j is obtained 
from 

Z =< 0\e'I^>^^''^'''\0 >= N fvA^VB^e'I^^-^^^^"^'^''' = e-'^^^^f'^K (8.72) 

The Lagrangian on the mean field level for the monopole field is written as 



Cdgl-mf = -^n- (d A A)Y[n-*id AB)], + -^[n- {d A B)Y[n-*id A A)l 



-^.[^■idAA)r-^,[n.{dAB)r 

+g(z7^9^ - e^^A^ . H - m)q + ^m^S^. (8.73) 



Integrating out the dual gauge field Bn, it becomes 



Cdgl-mf = -\Uur + ^A^K^.A^+qii^^d^ - e7^A^ ■ H - m)q (8.74) 



with 



1 Tn 



n m^ 
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— [-n^d'^g'''' + {n ■ dfg^"" + n^'rfd'^ - {n ■ d){n^'d'' + rfd") + n^d^'d" 



where X^y satisfies the relation, X^^ = X^^ and X^yd^ = X^yn" = 0. When we 
introduce the external source j in stead of dynamical quark part, we get 

Cdgl-mf = ]^A^^D-^A^ +]^A^ (8.76) 

= i(l'^ + /,D"^')D;J(1- + D^P]p) - \],Dnu. (8.77) 

Integrating out A^, we obtain the non-local current-current correlation as 

C, = —l^D^u. (8.78) 

where the propagator of the diagonal gluons is written as 



^;^' = Um.^^ - (1 - -)dM + . _ ""^^"^^ X^y (8.79) 



ae ^ '^J {n ■ dy + m\n^ 
1 / .d^du\ 1 m% n^ 

D,. = ^ [9>^^ + ("^ - ^)^^j - 9^92 + ^|(n . oy^^-" ^^-^^^ 

in the Lorentz gauge df^A'^ = 0. If the monopole does not condense f = or 
TJiB = 0, the second term in (|8.80|) disappears, and this propagator returns to 
familiar propagator derived in the perturbative sense. The nonperturbative effect 
is included in the second term. Indeed, the second term leads to the confinement 
potential as is shown from now on. 

The generating functional is given as 

Z = e*/^^'^'^ = e'^' = e-*^(j)/'^*. (8.81) 

Now we try to estimate 

S, = J qd'x = —Jj^D^-^Wx. (8.82) 

The static source of the quark with charge Qa located at a and the anti-quark Qb 
at b is given by 

U^) = -gAQaS'i^ - a) + 4<^'(x - b)} (8.83) 

and its Fourier transformation leads to 

Uk) = -g^Qae-^''-^ + Qbe-''^''). (8.84) 
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Figure 8.3: Confinement potential V^(r) as the function of the distance r between 
a static quark-antiquark. 



Substituting the propagator ( p. 801) and the source ( ^.841 ) into (|8.82| ), we get using 
the unit vector n 

S. = -I f S^M-k)D^''Uk) (8.85) 



' 2 J (27r)4 

1 r d^k 



2 J (27r)4 
1 r d^k 



J,{-k)[ 



k"^ — m\' 



flU 



—nil 



k^ — m?^ {n ■ k) 



" (9--^)liW 



■(4 + «l,e*T)(4 + 4e-*-i-)x 



n^ 

1 m 

+ 



2J (271)^'^'' ' ^"^ ;v^a , -^^^ ^'~^k2 + m| k' + m|(n-k)2J 

where r =b— a is relative vector between quark and anti-quark, and 27t6{0) = J dt 
is used. Here, we take n//r, because of the axial symmetry of the system and the 
energy minimum condition. The inter-quark potential is obtained from ( |8.81| ) as 



V{r; n) = 14ukawa(^) + V^mcar(^; n), 

where the yukawa term is 



(8.86) 



K 



yukawa V J 



d^k 

2J (2^ 



■tk-r 



iQa + Qbe'^''){Qa + Qbe 



ik-T\ 



k^ + m| 



Qa ■ Qb 



d^k 



e'^-^ Qa ■ Qb e-^^' 



(27r)3k2 + 



rrii 



Att r 



(8.87) 



and the linear term is 



linear (r;^) 



2 J (27r 



'^''' {Qa + Qbe^-'){Qa + Qbe~'^-'' ""^ 



k^ + m|(n-k)^ 



dt. 
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l/|^(gl + Q| + 24.4oos(U...))i^^ 



Introducing the ultraviolet cutoff m^ corresponding to the core of the flux-tube, we 
get the expression, 



Minear(r) = -^J ^{QI + QI + 2Q. ■ Q, cos(A;.r)) ln( ^; - , % 89) 



where kr is parallel component of k to the Dirac-string direction n. Since Qa • Qb is 



— \e^ for any color system, we finally get 



^ -ln( ^ „ 

127r r 247r ^ m^ 



H^) = -TTTT ■ — + ^^ H^^-T^y. (8.90) 



The string tension is given from the slope of the linear confinement potential as 

247r Y m^ y 

which corresponds to the vortex energy for the unit length in the typell supercon- 
ductor. Here, the gauge field mass m\ is neglected in the numerator in the logarithm 
function, because it is much smaller than the mass of the Higgs field (Cooper-pair 
field), m\ <C m\. Thus, the confinement potential arises from the second term of 
the gluon propagator ( p.SOj ). 

The parameters in the dual Ginzburg-Landau theory is determined so as to 
reproduce the interquark potential of heavy quarknium or results of the lattice QCD 
simulation; e = 5.5, monopole condensate v = 0.126GeV and interaction strength 
between the monopoles A =25. These parameters lead to the unit magnetic charge 
g = 2.3, the mass of the dual gauge field uib = 0.5GeV, monopole field mass 
m^ = 1.26GeV and the string tension a = l.OGeV/fm. Using these parameters, we 
compare in Fig.^]3| the phenomelogical potential like Cornell potential |J71| 



e^ 1 

l^Cornell = -7r^-+aer. (8.92) 

Sir r 

In the short range r < 0.2fm= (lGeV)~^, the Coulomb term dominates, while the 
linear potential dominates in the longer distance. 

Thus, in the DGL theory, the gluon propagator includes the long range cor- 
relation between quarks through the non-local operator , Lg • Such a long range 
correlation along the direction of the Dirac string leads to the linear confinement 
potential. 






DD 



Figure 8.4: The lowest-order polarization diagrams of the dual gauge field 5^, 
which is denoted by the wavy line, in the DGL theory. The dotted line denotes the 
monopole field. 

The DGL theory in the pure gauge is renormalizable similar to the abelian Higgs 



model [^, and is not asymptotically free on g in view of the renormalization group: 
g{p'^) become large as p^ increases. Hence, asymptotic freedom is expected for the 
QCD gauge coupling constant e owing to the Dirac condition: eijP') defined by 
e{p^)g{jp') = 477 become small as p^ increases. Thus, the DGL theory qualitatively 
shows asymptotic freedom on the QCD gauge coupling e p3| , [73[| . This asymptotic 
behavior in the DGL theory is consistent with QCD qualitatively, and seems a 
desirable feature for an effective theory of QCD. 

Next, we attempt to calculate the /3- function and the running coupling constant 
from the polarization tensor 11"^ (p) of the dual gauge field -B^ in the DGL the- 
ory. In particular, we are interested in the infrared region {p < IGeV), where the 
perturbative QCD calculation is not reliable. 

With the dimensional regularization |jl|, ^ by shifting the space-time dimension 
as d = 4 — e, we calculate the simplest nontrivial radiative correction from the 
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Figure 8.5: The running couphng constant g{fi) as the function of the renormahza- 
tion point /i. The sohd curve denotes the result in the DGL theory, which is directly 
calculated. The dashed curve denotes the leading-order perturbative QCD result, 
where g{fi) is obtained using the Dirac condition eg = Att. 
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Figure 8.6: The running coupling constant e{fi) as the function of the renormaliza- 
tion point fi. The solid curve denotes the result in the DGL theory, where e{fi) is 
obtained using the Dirac condition eg = An. The dashed curve denotes the leading- 
order perturbative QCD result. 
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monopole loop diagrams as shown in Fig.p^ 



™fe 



5''\gi^-r{p'g,.-P,Pu)- + 0{e') 



(8.93) 



where g is the bare dual-gauge coupling and /i the renormalization point. In the 
minimum subtraction scheme 0], the 0(l/e) divergence is eliminated by the counter 
term contribution, 



TrCaf) 



^u 



(P) 



Z3-1)<5<^''(pV-PmP. 



(8.94) 



where Z^ is the wave-function renormalization constant [|74| of the dual gauge field 



B 



A" 



ZM = 1 - 



(^/^^ 



-t\2 



327r2 



(8.95) 



Because of the Ward identity Zi = Z2 [|7^, the renormalized coupling constant is 
simply given by g{iJ?) = Z^d^Y^'^g. The /3- function |jl|, |7^ is then expressed as 



/? = /^:i-c/(/^) 
an 






-g{^^f + 0[g{fiY 



which determines the behavior of the running coupling g{fi) as 

ln(/iV/io) 



(8.96) 



(8.97) 



g^ifi) g^ifio) 327r2 

within the leading order. 

By summation of the multi-polarization diagrams, one obtains the final formula 
for the running coupling g{fi) including the higher order correction, 

1 



9 if^)= 9 (/^o) 



327r- 



ln(^ /^o)- 



8.98 



In the DGL theory, the QCD gauge coupling e{fi) defined by e{fi)g{fi) = An is simply 
expressed as 

I- -. ^ -1 



e\fi) = e2(/io) <^ 1 + 



e2(/xo) 



ln(/i//io) 



(8.99) 



We show in Fig. |8.5| the running coupling constants, g{jj,) and e(/i), in the DGL 



theory with the parameter set in Ref. [^: m;^=1.67GeV. Similarly in QED or 
the abelian Higgs model, we have imposed the renormalization condition as g{fi = 
2m^) = 7.9 [ e(/i = ^m^) = 1.59 ], which is taken to be consistent with the parameter 
g = 6.28 (e = 2.0) in Ref. [|75| in the infrared region (see Fig. ^.5D . This result shows 
that the gauge coupling e(/i) behaves as "walking coupling constant", which means 
the coupling varies slowly, even in the infrared region as /i < IGeV. 

Thus, owing to the Dirac condition eg = An, the DGL theory has asymptotic 
freedom nature on the gauge coupling constant e, where the "walking coupling 
constant" is predicted for e even in the infrared region. 



Chapter 9 

QCD Phase Transition at Finite 
Temperature 



The QCD vacuum is non-trivial at zero temperature. In this vacuum, quarks and 
gluons are confined in hadrons and the chiral symmetry is spontaneously broken. 
As the temperature increases, however, the color degrees of freedom in hadrons are 
defrozen. Above the critical temperature, the QCD vacuum is in the quark-gluon 
plasma (QGP) phase, where quarks and gluons move almost freely. 

The QCD phase transition is investigated in the various field of physics. As for 
the lattice QCD, the QCD phase transition is one of the most important subject in 
the computer science and studied using the Monte Carlo simulation. The simulation 
demonstrates that such a phase transition happens at about 280 MeV for the pure 
gauge case |7^ and at about 100~200MeV for the full QCD case [^]. In the early 



Universe, the QCD phase transition is considered to have occurred, i.e, the system 
was changed from the QGP phase into the hadron phase around the 200-300 MeV as 
the temperature decreases. As for the actual experience, in the recent years, some 
experimental groups are trying to create QGP in the laboratory using high-energy 
heavy-ion collisions. The RHIC (Relativistic Heavy Ion Collier) project is aimed at 
forming QGP and at studying its properties. 

In this chapter, we investigate the behavior of the color- confinement at high 
temperature by studying the change of the properties in the QCD vacuum with 
temperatures in terms of the dual Ginzburg-Landau theory |^, ^]. To this end, 
we concentrate on the pure-gauge QCD case, where glueballs appear as the physical 
excitation. Although such a pure gauge system is different from the real world, it is 
regarded as a proto-type of real QCD and is well studied by using the lattice QCD 



simulation ||4^. It is worth mentioning that our framework based on the DGL theory 
can be extended to include the dynamical quarks straightforwardly [^, |33[ keeping 
the chiral symmetry of the system, which is explicitly broken in the color-dielectric 



model 1 79 1 or in the lattice QCD with the Wilson fermion |42 
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9.1 Effective Potential Formalism at Finite Tem- 
perature 

In order to find the stable vacuum in the field theory, the effective potential for- 
malism is useful pO|]. The effective potential indicates the vacuum energy at zero 



temperature and corresponds to the thermodynamical potential at finite tempera- 
ture. Here, we investigate the effective potential in the path integral formalism in 
the DGL theory [^. The partition functional is written as 



Z[J] = JvxaVB,,exp Ujd^x{CT,GL - JE IXarij, (9.1) 



where we take the quadratic source term [^ instead of the standard linear source 
term [Q, ^ . As is well-known in the 0^ theory |1], |8y] , the use of the linear source 



term leads to an imaginary mass of the scalar field Xa in the negative-curvature 
region of the classical potential, and therefore the effective action cannot be obtained 
there due to the appearance of "tachyons". In this respect, there is an extremely 
advanced point in the use of the quadratic source term [^ , because the mass of the 



scalar field Xa is always real even in the negative-curvature region of the classical 
potential owing to the contribution of the source J to the scalar mass (see Eq.( p.5|) . 
Then, one obtains the effective action for the whole region of the order parameter 
without any difficulty of the imaginary- mass problem. Moreover, the effective action 
with the quadratic source can be formulated keeping the symmetry of the classical 
potential. Since this method with the quadratic source term is quite general, it 
is convenient to formulate the non-convex effective potential in the 0^ theory, the 
linear a model or the Higgs sector in the unified theory |jl|. 

The vacuum expectation value of Xa («=1,2,3) is the same value x due to the 



Weyl symmetry ||22[ , and therefore we separate the monopole field Xa into its mean 
field X s-nd its fiuctuation Xa as 

Xa = ix + Xa) exp (i^a)- (9.2) 

Here, the phase variables C,a have a constraint, J2a=i ^a = 0, where two independent 
degrees of freedom remain corresponding to the dual gauge symmetry [11(1)3 x 
U(l)g]m [23, 0- When monopoles condense, the phase variables ^a turn into the 
longitudinal degrees of the dual gauge field B^, which is the dual Higgs mechanism. 
Since we are interested in the translational-invariant system as the QCD vacuum, 
we consider the x-independent constant source J, which leads to a homogeneous 
monopole condensate. In the unitary gauge, the Lagrangian with the source term 
is rewritten as 

/:dgl - JELi \Xa\' = -CcKx) - 3 Jx^ _ 2x[2A(x' - v') + J] ELi Xa 

-\{d,B, - dj,f + \mlBl + Yi=A{d,Xa? - nilxl] (9.3) 

+ ELi{^'(4 ■ B,f{xl + nxa) - K^xxl + xi)}, 
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where C-c\{x) is the classical part, 

£d(x) = -3A(x'-t;2)2. (9.4) 

Here, the masses of Xa and B^ are given by 

ml = 2A(3x' -v^) + J = A\x\ m| = ?>g^x\ (9.5) 

where we have used the relation between the mean field x and the source J, 

J = -2\{x^-v'). (9.6) 

This relation is obtained by the condition that the linear term of Xa vanishes. It is 
remarkable that the scalar mass m^ is always real owing to the source J even in the 
negative- curvature region of the classical potential, x < w/v^- 

Integrating over B^ and Xa by neglecting the higher order terms of the fluctua- 
tions, we obtain the partition functional, 

Z[J] = exp [i j d'x{C,,{x) - 3Jxn)[Det(^D^i)]-MDet(^D-i)]-3/2, (9.7) 

where the exponents, —1 and —3/2, originate from the numbers of the internal 
degrees of freedom. Here, Db and D^^ are the propagators of B^ and Xa in the 
monopole condensed vacuum, 

V ""^B J k^ — m% + le k^ — m^ + le 

in the momentum representation. Hence, the effective action is given by the Legen- 
dre transformation |l[], 

r(x) = -t\nZ[J]+ fd^x3Jf 

= I (fxCdix) + i In Det(iD5^) + -i In Dei{iD-^). (9.9) 

The functional determinants are easily calculable in the momentum space, and we 
obtain the formal expression of the effective potential ||l| , 



Kff (X) = -r(x)/ / d'x = 3A(x^ - v'Y +3 / -0^ ln(m| - P - 



le] 



+ |/^M^^-A:'-^e). (9.10) 



In the finite-temperature system [Q, the partition functional Z is described by 
the Euclidean variables; Xq = —ir, and the upper bound of the r integration is 
/3 = 1/T with T being the temperature. Then, the /co-integration in the functional 
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determinant becomes the infinite sum over the Matsubara frequency fSQl- The effec- 
tive potential at finite temperatures physically corresponds to the thermodynamical 
potential, and is given by 

Kff (x; T) = 3X{f - v'f +3T E^=_oo / (0 ln{(2nvrT)2 + P + m%} 

+ |T E^=-oo / Hi2nnTr + P + m^ } (9.11) 

in the DGL theory. Performing the summation over n and the angular integration, 
we obtain the final expression of the effective potential at finite temperatures. 



Kff (x; T) = 3X{f - vr +3^ r dke In 1 - e-V^"^/^ 



+|^/o°^rfA:Pln l-e-v/^"^/^ , (9.12) 



where tub and m^ are functions of x as shown in Eq.( |9.5| ). Here, we have dropped 
the T-independent part (quantum fiuctuation), because we are interested in the 
thermal contribution to the QCD vacuum p2| . 

Before the numerical calculation, we examine the outline of the phase transition 
using the high-temperature (high-T) expansion. The effective potential is well ap- 
proximated by the high-T expansion, when the particle mass is much smaller than 
the temperature. Since the particle masses are almost zero for |x| ~ as shown 
in Eq. (|9.5|) , the high-T expansion is applicable for the effective potential around 

1x1 = 0. 

The temperature-dependent part of the effective potential for bosons is expressed 
as 

Vt{x]T) = — dke\n{l-e-^^^^^/^) 

Tl^ Jo 

= lljdyy'Hl~e-V^^), (9.13) 

with y"^ = -^ and a^ = ^. In the high-T expansion, Vt{x', T) is expanded in powers 
of a, 



Vt{x:T) = V{a\x)) 



dV 

a?=0 + 



/■oo 1 1 

/ , 9 J- J- 2 

a2=oa 



da? 
T^ /-"^ . o, . ... T^ f^ . o 1 1 



— / dyyHn{l - e"^) + -^ / dyy^^===^== 

vr^ Jo 277^ Jo vy + « gV?' +'^^ _ 1 

T^ 7r\ T^ 7r^m,2 vr^T^ T\m,^ , ,, 

■n^ 45 27r^ 6 i 45 12 i 



Thus, KflF in Eq. (|9.12|) is expressed as 



.2n2 ^ ^ , ^ /3 2 , o_2 



Kff(x; T) ^ 3A(x' - t;^)^ - ^^ + ^^(3^ + 3^b) (9-15) 
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in the high-T expansion. As the temperature increases, the effective potential at 
Ixl = is changed from the minimum point to maximum point at the lower critical 
temperature 

which satisfies the flat curvature condition as 

' 0. (9.17) 



^\x=o 



In fact, Ves{XiT) at x = is a local maximum for T < Tiow, while it becomes a 
(local) minimum and the system is (meta-) stable for T > Tiow For the 2nd-order 
phase transition, V^g has only one absolute minimum at each T because of continuous 
variation of the order parameter, and therefore Ti^w coincides the thermodynamical 
critical temperature Tc. On the other hand, for the Ist-order phase transition with 
a jump of the order parameter, there are two absolute minima at Tc-, and the second 
minimum at x = appears at the lower critical temperature Tiow, so that Tiow is 
lower than T^. In any case, using the high-T expansion, we derive the analytical 
expression for Tiow, which provides a lower bound for the critical temperature. 

9.2 Numerical Results on QCD Phase Transition 

Here, we show the results on the QCD phase transition which obtained from the 



numerical calculation 47 



We first show in Fig. p.l| the effective potential at various temperatures (ther- 
modynamical potential), K;ff(x; T), as a function of the monopole condensate x, an 
order parameter for the color confinement. The parameters, A = 25, v = 0.126GeV 
and g = 2.3, are extracted by fitting the static potential in the DGL theory to the 
Cornell potential |5 These values provide uib = O.SGeV and m^ = 1.26GeV at T = 
0. The (local-)minimum point of Ves{Xi T) corresponds to the physical (meta-)stable 
vacuum state. As the temperature increases, the broken dual gauge symmetry tends 
to be restored, and the monopole condensate in the physical vacuum, XphysiT), is 
decreased. A first order phase transition is found at the thermodynamical critical 
temperature, Tc — 0.49 GeV, and the QCD vacuum becomes trivial, XphysiT) = 0, 
for T > Tc- This phase transition is regarded as the deconfinement phase transition, 
because there is no confining force among colored particles in the QCD vacuum with 

Xphys(T) = 0. 

We show the behavior of the monopole condensate in the physical vacuum. 



Xphys(^), as a function of the temperature T in Fig.|9.2|. One finds Xphys= 0.126 



^ We examined several possible parameter sets, and found a small parameter dependence on 
our results shown in this paper. 
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Figure 9.1: The effective potentials at various temperatures as functions of the 
monopole condensate x- The numbers beside each curve are the temperatures. The 
absolute minimum points of the effective potentials are shown by crosses. 



GeV at T = 0, and the monopole condensate decreases monotonously up to Xphys 
= 0.07 GeV at the upper critical temperature Tup = 0.51 GeV, where the minimum 
at finite x disappears in V^ff(x; T). On the other hand, the local minimum is devel- 
oped at X = in V^fj(x;T) above the lower critical temperature Tiow = 0.38 GeV, 
which is analytically obtained by using the high-temperature expansion in previous 
section |8^, ^. The minimum value of Kff(x;T) at x = becomes deeper than 
that at finite x above the thermodynamical critical temperature Tc = 0.49 GeV. 
Here, we get the first-order phase transition because we have considered full orders 



in x^ as shown in Eg. ( p. 121) . On the other hand, Monden et al |82[ did not get the 
first-order phase transition due to the use of only the lowest order in x^ in the high- 



temperature expansion [Q, and therefore they had to introduce the cubic term in 
Xa in the Lagrangian. 

Here, we consider the possibility of the temperature dependence on the parame- 
ters (A,f ) in the DGL theory. The critical temperature, Tc = 0.49 GeV, seems much 
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Figure 9.2: The monopole condensate Xphys(^) at minima of the effective potential 
as a function of the temperature T. The solid curve denotes XphysiT) corresponding 
to the confinement phase, which is the absolute minimum up to Tq =0.49 GeV and 
becomes a local minimum up to Tup =0.51 GeV. A minimum appears at x = 
above Tiow =0.38 GeV and becomes the absolute minimum above Tq =0.49 GeV. 
The dot-dashed curve denotes the value of x at the local maximum. 



larger than one of the recent lattice QCD prediction, Tq = 0.26 ~ 0.28GeV, which 
is, for instance, estimated from the relation, Tcj \fo ~ 0.62 |Q and the string ten- 
sion a = 0.89 ~ l.OGeV/fm. However, we should remember that the self-interaction 
term of Xa has been introduced phenomenologically in the DGL Lagrangian. In par- 
ticular, the asymptotic freedom behavior of QCD leads to a possible reduction of 
the self-interaction among monopoles at high temperatures. Hence, we use a simple 
ansatz for the temperature dependence on A, 



A(T) = A (: 



Tc-aT 



T, 



c 



(9.18) 



keeping the other parameter v constant. Here, the constant a is determined as 
a = 0.89 so as to reproduce Tc = 0.28GeV. (We take A(T) = for T > Tc/a.) The 



results for the monopole condensate Xphys(7') are shown in Fig.p.3|. The qualitative 
behavior is the same as in the above argument with a constant A. We find a weak 
first-order phase transition in this case also. Here, we find a large reduction of the 
self-interaction of the monopoles near the critical temperature Tc'- A(T ~ Tc) — 2.7 
is considerably smaller than A(T = 0) = 25. 

Next, we investigate the variation of the masses of the dual gauge field B^ and the 
monopole field Xa at finite temperatures. Here, Xa would appear as the color-singlet 
glueball field with 0+ [E2|, ^. These masses uib and rriy, at the finite temperature 
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Figure 9.3: The nionopole condensate Xphys(^) at minima of the effective potential 
as a function of the temperature in the case of variable A(T) so as to reproduce 
Tc =0.28GeV. The meanings of the curves are the same as in Fig.p]2. 



T are given by 

iriBiT) = VSgXptysiT), m^{T) = 2^fMT)xphys{T) (9.19) 

as shown in Eq. (|9.5| ). In Fig. |9.4| , we show msiT) and m^{T) as functions of the 



temperature T using variable A(T) in Eq. (|9.18|) . (Their behaviors are almost the 
same as the case of a constant A except for the difference of the value of T^.) It is 
worth mentioning that mB{T) and m^{T) drop down to mB,m^ ~ Tc(=0.28GeV) 
from mB,Tn^ ~ 1 GeV near the critical temperature Tq- In other words, the QCD 
phase transition occurs at the temperature satisfying mB,Tn^ ~ T, which seems 
quite natural because the thermodynamical factor l/{exp(v^P~+l7?/T) ±1} be- 
comes relevant only for ?7i ~ T. Thus, our result predicts a large reduction of the 
dual gauge field mass tub and the monopole field mass m^ near the critical tem- 
perature Tc- It is desirable to study the change of the scalar glueball mass at finite 
temperatures, especially near Tc, in the lattice QCD simulation with the larger 
lattice size and the higher accuracy. 

We investigate the string tension a at finite temperatures, since a is one of 
the most important variables for the color confinement, and controls the hadron 
properties through the inter-quark potential. We use the expression of the string 
tension cr(T) provided by SST ||33| , 



a{T) 



e^m|(T), (ml{T)+ml{Ty 



247r 



In 



m' 



m 



(9.20) 



where the masses mB{T) and m^{T) are given by Eq. (|9.19| ). The results are shown 



in Fig. |9.5| as a function of the temperature T. In the case of constant A, the string 
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Figure 9.4: The dual gauge field mass mB{T) and the monopole field mass m^{T). 
The solid lines denote the case of variable A(T) with a constant v. The dashed lines 
denote the case of variable v{T) with a constant A. A large reduction of these masses 
is found near the critical temperature. The dotted line denotes m = T. The phase 
transition occurs at the temperature satisfying mBiTTty ~ T. 



tension a{T) decreases very gradually up to the temperature, Tup= 0.51 GeV. On the 
other hand, in the case of variable A(T), the string tension cr(T) decreases rapidly 
with temperature, and cr(T) drops down to zero around Tc = 0.28 GeV. Hence, 
one expects a rapid change of the masses and the sizes of the quarkonia according 
to the large reduction of (j{T) at high temperatures. We plot also the pure-gauge 
lattice QCD results for the temperature dependence of the string tension by black 

^ . We find our results with variable A(T) agree 



dots 



with T( 



c 



0.28GeV 



with the lattice QCD data. 

We discuss further the temperature dependence of the parameters (A,f ) in the 
DGL theory. Definitely, we should follow the lattice QCD data for this determina- 
tion as the case of the Ginzburg-Landau theory of superconductors extracting the 
temperature dependence from experiments. Since there exists the lattice QCD data 



on the string tension a [^, we try to reproduce a by taking a simple ansatz on A 
and V. We try the following ansatz, 

'Tc - aT' 



BiT) = 3X{T)v\T) = 3Xv' 



T< 



c 



(9.21) 



where the constant a is determined so as to reproduce Tq = 0.28GeV. (We take 
B(T) = for T > Tc/a.) The variable B{T) corresponds to the bag constant, the 
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Figure 9.5: The string tension cr{T) as a function of the temperature T. The sohd 
and dashed hues correspond to the variable A(T) case with a constant v and the 
variable v{T) case with a constant A, respectively. The constant (A, v) case is also 
shown by the thin line. The lattice QCD results in the pure gauge are shown by 
black dots near the critical temperature. 



energy-density difference between the nonperturbative vacuum {\xa\ y^ 0) and the 
perturbative vacuum {\xa\ = 0) in the DGL theory; see Eq.( |9.12|) . The ansatz (|9.21 



suggests the reduction of the bag constant at high temperatures, which provides the 
swelling of hot hadrons by way of the bag-model picture. Since we have already 
examined a typical case for variable A(T) keeping v constant, we show here another 
typical case for variable v{T) keeping A constant. The string tension cr(T) in the 
variable v{T) case with a = 0.91 is shown by the dashed line in Fig. |9.5| . We find 
almost an identical result and find again a good agreement with the lattice QCD 
data. Other combinations on A(T) and v{T) under the relation (|9.21| ) also provide 
equally good results on cr(T). 

Finally, we investigate the relation between the scalar glueball mass m^lT) and 
the string tension (t{T). For variable A(T) keeping v constant, one finds from 
Eq.(p.20|) an approximate relation. 






1.6, 



(9.22) 



near the critical temperature Tc- On the other hand, for variable v{T) keeping A 
constant, the scalar glueball mass at finite temperatures, m^{T), is shown by the 
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dashed line in Fig.|9.4|, and Eq.( |9.20|) leads to a simple relation 




for the whole region of T. Thus, the DGL theory suggests a proportional relation 
between the scalar glueball mass and the square root of the string tension at least 



near Tc- It is worth mentioning that Engels et al. |83| obtained a similar relation 



mQB{T) = (1.7 ± 0.5)wcr(T), for the lowest scalar glueball at finite temperatures 



from the thermodynamical study on the SU(2) lattice gauge theory. Eqs. (|9.22|) and 
( |9.23[ ) can be examined from the thermodynamical study on the glueball mass in the 
lattice QCD, which may also reveal T-dependence on the parameters in the DGL 
theory. 

9.3 Hadron Bubble Formation in Early Universe 

The remaining part of this section, we consider the application of the DGL theory 
to big bang [Q. As Witten proposed |8^, if the QCD phase transition is of first 



order, the hadron and the QGP phase should coexist in early Universe. Such a 
mixed phase may cause the inhomogeneity of the Universe in the baryon number 
distribution. This inhomogeneity affects the primordial nucleo-synthesis [^ . 



As a result of the 1st order phase transition, hadron bubbles appear in the 
QGP phase near the critical temperature. We now consider how hadron bubbles 
are formed in the DGL theory. In the supercooling system, the free energy of the 
hadron bubble with radius R profile x(^; R) is written using the effective potential 
at finite temperature, 

E[x{r-R)]=47iJ^^drr'{3{^^^r + V,six;T)}. (9.24) 

We use the sine-Gordon kink ansatz for the profile of the monopole condensate, 

X{r; R) = xh tan"^ e^R-)l^ / tan^^ e^'\ (9.25) 

where the thickness of the surface 5 is determined by the free energy minimum 
conditions. The result is shown in Fig. |9.6| . The monopole condensate x('"; R) is 
connected smoothly between inside and outside the bubble. The energy density of 
the hadron bubble is shown in Fig.|9r^. It is negative inside and positive near the 



boundary surface. The total energy is roughly estimated as the sum of the surface 
term (corresponding to the positive region near the surface) and the volume term 
(corresponding to the negative region inside the bubble). 

The energy of the hadron bubble with radius R is shown in Fig. |9.^ . The bub- 
ble whose radius is smaller than critical radius Re collapses. Only larger bubbles ( 
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Figure 9.6: The profile of the monopole condensate in the hadron bubble. There 
is the QGP phase without monopole condensation outside the bubble, while hadron 
phase with monopole condensation remains inside the bubble. 
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Figure 9.7: The energy density of the hadron bubble. It is negative inside and 
positive near the boundary surface. 
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Figure 9.8: The total energy of the bubble is plotted as a function of the hadron 
bubble radius/?. The total energy is roughly estimated as the sum of the volume 
term and the surface term. 
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Figure 9.9: The critical radius Re, corresponding to maximum of the energy in 



Fig.|9^, is plotted as a function of temperature. 



Figure 9.10: The scenario of the QCD phase transition in the early universe. The 
shaded and white regions denote the hadron and QGP phases, respectively, (a) 
Slightly below Tc, only large hadron bubbles appear with very small creation rate, 
(b) Hadron bubbles expand with radiating shock wave, (c) Near Tiow, many small 
hadron bubbles are created, (d) The QGP phase pressured by the hadron phase is 
isolated. 



R > Re) are found to grow up from the energetical argument. However, the cre- 
ation of large bubbles is suppressed because of formation probability. In the bubble 
formation process, there exists a large barrier height h of the effective potential and 
therefore the creation of large bubbles needs the large energy fluctuation above the 
barrier height. Such a process is suppressed because of the thermal dynamical fac- 
tor (proportional to bubble formation rate), exp(— |7r_R|?/i/T). Thus, the only small 
bubbles are created practically, although its radius should be larger than R^ energit- 
ically. The temperature dependence of the critical radius is shown in Fig.|9.9|. In the 



temperature region of the supercooling state, i.e, Tiow < T < T^ the hadron bub- 
bles are created. As the temperature decreases, the size of hadron bubble becomes 
smaller, however the bubble formation rate becomes larger. 

From these results, we can imagine how the QCD phase transition happens in 
the big bang scenario [0], as shown in Fig. p.lO] . At the first stage slightly below 



Tc, only large bubbles are created but its rate is quite small. As temperature is 
lowered, smaller bubbles are created with much formation rate. During this process, 
the created hadron bubbles expand with radiating shock wave which reheats QGP 



phase [87]. Near Tiow niany small bubbles are violently created. Finally QGP phase 



is isolated like the bubble 187 . Such an evolution of the hadron bubble can be 



obtained from the numerical simulation using the DGL theory. 
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Thus, using the quadratic source term instead of the linear source term, we have 
obtained the effective potential in the DGL theory in the all |x| region. Thermal 
effects reduce the monopole condensate and lead to deconfinement phase transi- 
tion. Since the critical temperature is found very large, Tc ~ O.SGeV, with the 
parameters unchanged at finite temperature, we introduce a large reduction of the 
self-interaction between monopoles so as to make the phase transition temperature 
about 0.28GeV, which is suggested in the lattice QCD. We find large reduction of 
string tension with the temperature in accordance with the lattice QCD results. We 
predict that the glueball mass decreases considerably near the critical temperature, 
which is to be checked by lattice QCD simulation and experiment. Based on this 
effective potential at finite temperature, we further investigated properties of hadron 
bubbles created in the early Universe and discussed the hadron bubble formation 
process from estimation of the size of hadron bubble at various temperatures. 



Chapter 10 

Application to Quark Gluon 
Plasma 



In the recent years, some experimental groups are trying to create quark gluon 
plasma (QGP) as the new form of matter in the laboratory by high-energy heavy-ion 
collisions. The RHIC (Relativistic Heavy Ion Collier) project will start in the next 



year. The scenario of producing QGP is based on Bjorken's picture [^. Just after 
heavy ions pass through each other, many color-fiux-tubes are produced between the 
projectile and the target, and pulled by them as shown in Fig. |10.1| (a). Usually, it 
is guessed that these flux-tubes are cut into several pieces through quark- ant iquark 
pair creations, and these short fiux-tubes, which behave as excited 'mesons', are 
thermalized by stochastic collisions among themselves. If the energy deposition is 
larger than a critical value, the thermalized system becomes the QGP phase, whereas 
if it is lower, the system remains to be the hadron phase. 

The features of the multi color-flux-tube system strongly depend on their den- 
sity of the fiux-tubes created by hard process in early stage. When the fiux-tube 
number density is low enough, the system is approximated as the incoherent sum 
of the individual fiux-tube. Its evolution is expected to be superposition of random 
multiple hadron creations of many color-fiux-tubes produced between many nucleon- 
nucleon pairs. On the other hand, when the fiux-tube number density is sufficiently 
high, many flux-tubes overlap each other and would be melted into a big flux-tube. 
During this process, each flux-tube loses its individuality and the whole system can 



be regarded as a huge flux-tube between heavy-ions like a condenser ||8^ 



In this chapter, we would like to study the properties of the multi-flux-tube 
system [Q. QCD is very hard to deal with in the infrared region analytically due 
to the breakdown of the perturbation technique. Moreover, for such a large scale 
system there is a severe limit on the computational power even in the lattice QCD 
simulation. For the study of the QGP formation, the DGL theory would provide 
a useful method, because it describes the properties of the color-electric flux tube, 
which are important in the pre-equillibrium system just after the ultrarelativistic 

128 



Figure 10.1: The scenario of the QGP formation in high-energy heavy-ion colh- 
sions. (a) There appear many color-flux-tubes between the projectile nucleus and 
the target nucleus just after the collision, (b) When the distance between the two 
nuclei becomes large, There appears the pair creation of quark and anti-quark, (c) 
Many created quarks and anti-quarks make frequent collisions to form a thermal 
equilibrium and form QGP, when the energy density is larger than a critical value. 

heavy-ion collisions. 

10.1 Formalism on Multi-flux-tube System 

In this section, we formulate the multi-flux-tube system in the DGL theory. In 
this theory, the color-flux-tube is described as the dual version of the Abrikosov 
vortex |]70[. Different from the multi- vortex system in the superconductor, however. 



there are some kinds of the color-fiux-tubes corresponding to the kinds of color. 
Furthermore, in the superconductor, the direction of the fiux-tube is all same and 
the system can be described only by the ground state such as triangle lattice system, 
while in the QCD vacuum many color fiux-tubes distribute randomly and the system 
includes the highly excited states. Here, we simplify such a complex system and 
discuss qualitatively. For simplicity, we consider a single color charge system using 
the dual Ginzburg-Landau Lagrangian Eq.( |8.67| ). 

We consider two ideal cases of multi-fiux-tube penetrating on a two dimensional 
plane. The directions of all the color-fiux-tubes are the same [Fig. |10.2| (a)] in one case 



or alternative[Fig.[TO]^(b)] in the other case. When the fiux-tubes are long enough, 
the effect of the flux edges is negligible. Hence, taking the direction of the fiux- 
tubes as the 2;-axis, the system is translationally invariant in the 2;-direction and is 
essentially described only with two spatial coordinates {x,y). For the periodic case 
in the {x, y) coordinate, we can regard the system as two fiux-tubes going through 
two poles (north and south poles) of the S^ sphere. For the system of fiux-tubes 
with all the directions being the same, we take two flux-tubes passing through the 
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two poles on S"^ sphere as shown in Fig. l0.2| (a) (which we call the two flux-tubes 



system). For the alternative case, on the other hand, we take a flux-tube coming in 
from the south pole and the other leaving out through the north pole (which we call 
flux-tube and anti-flux-tube system). Such a prescription leads the exact solution 
for the periodic crystal of the sine-Gordon kinks, and also provides a simple but 
good description for the flnite density Skyrmion system studied by Manton [Q. 

The two color-flux-tube system on the sphere S"^ with radius R corresponds to 
the multi-flux-tube system with the density p = 1/{2ttR'^). Introducing the polar 
coordinates {R, 9, ip) on S'^, we consider the static solution satisfying 

Eo = 0, B = B{9)e^ ^ #^e^, X = x(^)e^"^ (10.1) 

rism u 

where we have used the axial symmetry of the system. Here the electric fleld pene- 
trates vertically on the sphere surface, E//ej.; 

E = V X B = Eer, (10.2) 

which corresponds to the electric fleld penetrating vertically on the plane, E/ /e^. 
The fleld equations are given by 

^ "^ ^sin^^) - [i;7i27^(^ - 9mr + 2X{f - v')]x = 0, (10.3) 



R'^sinOdO^ dO' 'R'^sin^O 



d ^ 1 d f^f^^s^ , Sg- ^ ^f)/n\\-2 



Consider the closed loop C on S"^ with a constant polar angle 9 = a and G [0, 2tt), 
the electric flux penetrating the area surrounded by the loop C is given by 

$(a) = f E-dS= fv xB-dS= iB-d\ = 2TxB{a). (10.5) 

The boundary conditions for the two flux-tubes system as shown in Fig. |10.2| (a) are 
given as 

$(a) = 2TiB{a) = 0, x(") = as a ^ 0, (10.6) 

~ 1 /" ^TTTl 71 

$(a) = 2'KB{a) = - ER^dVl = ± as a -> - ± e, (10.7) 

Here, n corresponds to the topological number of the flux-tube, which appears also in 
the vortex solution in the superconductivity. This boundary condition at a ^ ^ ± e 
has a discontinuity for the dual gauge fleld B. Because the electric flux leaves out 
from the two poles, there should be some sources to provide the electric flux. In 
this case, the Dirac-string like singularity appears on the Q = 7r/2 line, through 
which the electric flux comes into the sphere from long distance. For the flux-tube 
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(a) 





(b) 



Figure 10.2: (a) A multi-flux-tube system with the same direction of the flux-tubes 
is approximated by the two color-flux-tubes going out from the north and the south 
poles on a sphere S"^. (b) A multi- flux-tube system, where flux-tube direction is 
alternative, is approximated by the flux-tube and 'anti-flux-tube' system penetrating 
on S'^ with the color-flux going in from the south pole and leaving out from the north 
pole. 
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and anti- flux-tube system as shown in Fig. l0.2| (b), the boundary condition around 
^ = 0, TT is given as 



^(a) = 2TTB(a) = 0, 



X{a) = 



as a ^ 0, IT. 



10.8) 



In this case, there does not appear the Dirac-string like singularity, since the electric 
flux is conserved. The free energy for the unit length of the color-flux-tube is written 
as 

2 / , , _ \ 2" 



+ 



R^dn 
R^dVL 



2 yR^sinede^^^^) ^\Rde) 



i?2 sin^ e 



[n-gB{e)Yx' + ^[r-v 



27r sin^d^ 
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;io.9) 



First, we consider a limit of i? -^ oo, which corresponds to the ordinary single 
vortex solution. Introducing a new variable p = i?sin6', the free energy in the limit 
R ^ p{6 ~ 0) is written as 
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2npdp 
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and the fleld equations are 
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;f -;f {pB{p)) + ^(n - 9pB{p))t = 0, 
dp pdp p 



with the boundary condition. 



$ = 2npB{p)\ 



27m 



as p ^ oo. 



(10.12) 
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Above equations, (|10.10[ - p~0.12|) , coincide exactly with those of ordinary single vortex 
solution in the cylindrical coordinate. Thus, we get the desired results. 
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One can analytically investigate the dependence of the proflle functions(i?(^), x(^)) 
on the flux-tube number density. For this purpose, we express the free energy as 



fo + Ie- ^ + fx- ^^ 



where fo, fs, and/^ are R independent functions and written as 
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(10.14) 

(10.15) 

(10.16) 
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In the large R case, which corresponds to the small color-flux-tube number density 
in the original multi-flux-tube system, the third term fxR"^ is dominant. Hence, the 
free energy F is minimized as x ~ "^j that is, the monopole tends to condense, and 
then the color electric fleld is localized only around ^ = (north pole) and 6 = ir 
(south pole) as shown in Fig. |10.3| . On the other hand, in the small R case, the 
second term Je/R"^ is dominant. There is a constraint on the total flux penetrating 
on the upper sphere. 



$ = / E {6)271 R' sine de 
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= C with t = COS0. 
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Hence, one flnds the equation. 



fE oc / E{e)hmede 

Jo 



E{tfdt 



fdt= I {{E{t)~Cf}dt + C''>C\ (10.20) 
10 Jo 

This condition leads to the uniform color electric fleld E = C, which provides 
the minimum of /g. Thus, the color electric fleld tends to spread over the space 
uniformly. 

Here, we consider the critical radius Re of the phase transition to the normal 
phase, where the monopole condensate disappears. There are three useful inequali- 
ties on /e, fx, and F, 



Ie > 27r( 



(10.21) 



< /y < ^ttXv^ 



(10.22) 
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Figure 10.3: The color-electric field E{6) (solid curve), the monopole condensate 
x{0) (dashed curve) and the dual gauge field B{9) (dotted curve) are plotted as 
functions of the polar angle 6 for the low density case with R =4.0fni. 
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The equality is satisfied in Eq. (|10.23| ) 
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Using the inequalities equations( |10.21| )-( |10.23| ), R^ is larger than the critical R^ 
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For R > Re, there exists a nontrivial inhomogeneous solution, which differs from 
the normal phase. For R < R^, homogeneous normal phase provides the minimum 
of F. Thus, Re is the critical radius of the phase transition from the fiux-tube phase 
to the normal one. In this case, the critical radius and the electric field are given by 
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respectively. 
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10.2 Numerical Results on Multi-flux-tube Sys- 
tem 



We start with showing the low density case of two flux tubes system in Fig. l0.3| . The 



parameters of the DGL theory are flxed as A = 25, v = 0.126GeV, which lead the 
masses tjib = O.SGeV and m^ = 1.26GeV |Q. This parameter set provides the flux- 



tube radius rpT ~0.4fm and the suitable interquark potential with the string tension 
as a =lGeV/fm. The monopole condensate x(^)) the dual gauge fleld B{9) and the 
color electric fleld E{6) are plotted as functions of the polar angle 6. The electric 
fleld E{9) is localized around the two poles(6'=0 and vr) and drops suddenly as 9 
deviates from the two poles. The monopole condensate vanishes at the two poles and 
becomes constant in the region away from these poles. This behavior corresponds 
to the case of independent two vortices in superconductivity. Different from these 
physical quantities, the dual gauge fleld B{9) is discontinuous at ^ = 7r/2. There 
should be Dirac-string like source to provide, because the electric flux leaving out 
from the two poles. It should be noted that the system has the reflection symmetry 
on 6^ = I plane. 

We show now in Fig. l0.4| the number density dependence of the flux-tubes. For 



the large radius of the sphere, {R >2fm), the color electric flux is localized at 6' = 
and 71 and there the monopole condensate vanishes, while the x becomes constant 
X — V around 6* = |. As the radius R decreases, the electric flux, locahzed at 6^ = 0, 
9 = IT, starts to overlap, and the value of the monopole condensate x becomes small. 
The electric fleld E{9) becomes constant and x vanishes below a critical radius Re. 
We show in Fig. |10.5| the monopole condensate at 6* = | (the maximum value of 
the monopole condensate) as a function of the sphere radius R (flux-tube number 
density p = l/{2nR^)). The (flrst order) phase transition occurs and the system 
becomes homogeneous normal phase above the critical value of the flux-tube number 
density. 

Here, we compare the free energy of two flux tube system with that of inhomo- 



geneous system in Fig. p.0.6| . At large i?, the former is smaller and the system favors 
the existence of two flux tubes. As R becomes smaller, the energy difference of the 
system becomes smaller. Two flux-tubes melt and the electric fleld are changed to 
be homogeneous below the critical radius Re = 0.35fm, which corresponds to the 
critical density pc = l/(27ri?^) = 1.3fm~^. This critical density agrees with the 
analytic estimation in Eq.( |10.26| ), Pc = y^^ = 1.3fm~^. 

We discuss now the system of flux-tube and anti- flux-tube with opposite direction 



placed at 9 =0 and 9 = tc respectively as shown in Fig. p.0.7| . At low flux-tube number 
density ( R > 2fm), the flux-tube is localized at 6' = 0, while the anti-flux-tube at 
9 = TT. The monopole condensate x(^) vanishes at the two poles {9 = 0, tt) and 
becomes constant away from these poles. As R decreases, the electric flux starts to 
cancel each other and the monopole condensate becomes small. We also compare 
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Figure 10.4: The flux-flux system on S^. The color-electric field E{6) and the 
monopole condensate x(^) are depicted as functions of the polar angle 6 for the 
three radii; R =2.0fm, 0.5fm and 0.347fm. Below the critical radius -Rc=0.347fm, 
the color-electric field E becomes constant and the monopole condensate x vanishes 
entirely. 
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Figure 10.5: The R dependence of the monopole condensate. Here, the flux tube 
density is given by l/(27ri?^). The monopole condensate at ^ = | decreases, as the 
radius R becomes smaller, and vanishes below R = Re. 
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Figure 10.6: The free energy in flux-flux system (solid curve) and uniform system 
(dashed curve). Because the lower free energy system is realized, the inhomogeneous 
system is changed into a homogeneous system below the critical radius. Re = 0.35fm. 
Here, the critical flux-tube density reads pc = 1.3fm~^. 
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the free-energy of the flux-tube and anti-flux-tube system with the homogeneous 
system, in which both monopole condensate and electric field are vanished. The 
critical radius, Re = 0.31fm, is similar to the value of the two flux-tubes system. 

Thus, we have found in both cases that the solution behaves as two independent 
flux-tubes at a large R, i.e., a small number density p. As the radius R decreases, 
the monopole condensate decreases and eventually vanishes at a critical density, 
where the color-electric field E becomes uniform. The critical density pc is found 
as pc =1.3fm~^ = (1.14fm)~^ for the two color-flux tube system; pc =1.7fm~^ = 
(1.3fm)~^ for the flux-tube and anti-flux-tube system. Such similar values in both 
cases suggest that an actual flux-tube system would become uniform around similar 
density to pc ~1.5fm~^, because realistic flux-tube system includes the flux-tubes 
and anti-flux tubes randomly, which would correspond to an intermediate system 
between the above two ideal cases. Thus, the configuration of the color-electric field 
and the monopole field depend largely on the flux-tube density. 

As discussed above, many flux-tubes are melted around pc ~ 1.5 fm~^. Let us 
discuss here in case of central collision between heavy-ions with mass number A. The 
nuclear radius R is given hj R = RqA^^^, where RQ=1.2iia corresponds to the nuclear 
radius, and the normal baryon- number density is x^ — T^- ■'■^ ^^^ central region. 



_A 

one nucleon in the projectile makes hard collisions with ^A^/^{= 27rRlA^/^ x pg) 
nucleons in the target, because the reaction volume is nRl x 2R = 2ttRqA^^^. Hence, 
nucleon-nucleon collision number is expected as {^A^^^y = |A^/^ per the single- 
nucleon area ttR^ between projectile and target nuclei. Assuming one flux-tube 
formed in one nucleon-nucleon hard collision, the flux-tube density is estimated as 

9^2/3 ,2/3 

p = '^^^i = Q 4fj^-)2 • For instance, p would be 4.5, 5.8 and 17.5 fm~^ for A = 27, 40 
and 208. This would indicate that the scenario of creating large sizes of flux tube 
becomes much relevant for A-A collision with larger A. 

It would be important to reconsider the process of the QGP formation in terms 
of the flux-tube number density. When the flux-tube density is low enough, the 
flux-tubes are localized. Each flux-tube evolution would be regarded as the multi- 
creation of hadrons in the high energy p-p collision via the flux-tube breaking. In 
this process, q-q pair creation plays an essential role on the QGP formation, which 
is the usual scenario. 

On the other hand, for the dense flux-tube system, neighboring flux-tubes are 
melted into a large cylindrical tube, where monopole condensate disappears. Such a 
system, where the color electric field is made between heavy-ions, becomes approx- 
imately homogeneous and is regarded as the 'color condenser'. In this case, large 
homogeneous QGP may be created in the central region. 

In the actual case, however, the variations and directions of the color-flux-tubes 
are random. For instance, in the peripheral region, flux-tubes would be localized and 
are broken by quark-antiquark pair creations. In the central region, dense flux-tubes 
are melted by annihilation or unification [RO] of flux-tubes, which will be discussed 
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Figure 10.7: The flux the and anti-flux-tube system on S"^. The color-electric fleld 
E{6) and the monopole condensate xiO) are depicted as functions of the polar angle 
9 for the three radii; R =2.0fm, 0.5fm, and 0.315fm. Below R =0.315fm, both the 
color electric flux and the monopole condensate vanish entirely. 
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in the next section. In this region, a huge system of dynamical gluons appears, 
because many dynamical gluons are created during this process. Thermalization of 
such quarks and gluons leads to QGP. Thus, the process of QGP formation should 
depend largely on the density of created flux-tubes, which is closely related to the 
incident energy, the impact parameter and the size of the projectile and the target 
nuclei. 



10.3 Flux-tube Interaction and QGP Formation 
Process 

Finally, we investigate the interaction among various color-flux-tubes, which is con- 
sidered to be important for the QGP formation process. Each flux tube is character- 
ized by the color charge Q [^, |91| at its one end. To classify sorts of the flux tube, 
we call the flux tube with a red quark (i?) at its one end as ''R-R flux tube" , and so 
on. In this case, the "direction" of the color-electric flux in the flux tube should be 
distinguished. For instance, R-R flux tube is different from R-R flux tube in terms 
of the flux direction. We study the interaction between two color-electric flux tubes 
using the DGL theory. The color-electric charges at one end of the flux tubes are 
denoted by Qi cind Q2- We idealize the system as two sufficiently long flux tubes, 
and neglect the effect of their ends. We denote by d as the distance between the 
two flux tubes. For d ^ friZ^, the interaction energy per unit length in the two flux 
tube system is estimated as 

i^int ^ ^-^^^^mlK,{msdl (10.28) 

where Kq{x) is the modifled Bessel function. Here, we have used the similar calcu- 
lation on the Abrikosov vortex in the type-II superconductor 



As shown in Fig.p.0.8|, there are two interesting cases on the interaction between 



two color-electric flux tubes. 

1. (a) For the same flux tubes with opposite flux direction (e.g. R-R and R-R), 
one finds Qi = —Q2 i-e. Qi-Q2 = — e^/3, so that these flux tubes are attracted 
each other. It should be noted that they would be annihilated into dynamical 
gluons in this case. 

2. (b) For the different flux tubes satisfying Qi-Q2 < (e.g. R-R and B-B), one 
finds Qi ■ Q2 = — e^/6, so that these flux tubes are attractive. In this case, 
they would be unifled into a single flux tube (similar to G-G flux tube). 

Based on the above calculation, we propose a new scenario on the QGP for- 
mation via the annihilation of the color-electric flux tubes. When the flux tubes 



Figure 10.8: The annihilation process of the color-electric flux tubes during the 
QGP formation in ultrarelativistic heavy-ion collisions, (a) The same flux tubes with 
opposite flux direction are attracted each other, and are annihilated into dynamical 
gluons. (b) The different flux tubes (e.g. R-R and B-B) are attractive, and are 
unified into a single flux tube. 

are sufficiently dense in the central region just after ultrarelativistic heavy-ion col- 
lisions, many flux tubes are annihilated or unified. During the annihilation process 
of the flux tubes, lots of dynamical gluons (and quarks) would be created. Thus, 
the energy of the flux tubes turns into that of the stochastic kinetic motion of 
gluons (and quarks). The thermalization is achieved through the stochastic gluon 
self-interaction, and finally the hot QGP would be created. Here, the gluon self- 
interaction in QCD plays an essential role to the thermalization process, which is 
quite different from the photon system in QED. 

In more realistic case, both the quark-pair creation and the flux-tube annihilation 



would take place at the same time. For instance, the flux tube breaking |^, 0, |9^ 
would occur before the flux tube annihilation for the dilute flux tube system. On 
the contrary, in case of the extremely high energy collisions, these would be lots 
of flux tubes overlapping in the central region between heavy ions, and therefore 
the flux tube annihilation should play the dominant role in the QGP formation. In 
any case, the DGL theory would provide a calculable method for dynamics of the 
color-electric flux tubes in the QGP formation. 



Chapter 11 

Summary and Concluding 
Remarks 



In the basis of the dual superconductor picture, we have systematically studied 
the confinement phenomena using the lattice QCD Monte-Carlo simulation, the 
monopole-current dynamics and the dual Ginzburg-Landau (DGL) theory, an in- 
frared effective theory of QCD. In the dual Higgs theory, color confinement is un- 
derstood by one-dimensional squeezing of the color-electric fiux in the QCD vacuum 
through the dual Meissner effect caused by monopole condensation. For the con- 
struction of the dual superconducting theory from QCD, there are two large gaps 
on speciality of the abelian sector and the appearance of monopoles, however, these 
gaps are expected to be fulfilled by taking the 't Hooft abelian gauge fixing, which is 
defined by the diagonalization of a suitable gauge-dependent variable, <l>[A^(x)]. In 
the abelian gauge, the SU(iVc) gauge theory is reduced into the \J{1)^'=~^ gauge the- 
ory including color-magnetic monopoles, which topologically appears corresponding 
to the nontrivial homotopy group, n2(SU(A^c)/U(l)^'=~^) = Z'^'=~^. In this gauge, 
the diagonal and the off-diagonal gluon components behave as the abelian gauge field 
and the charged matter field, respectively, in terms of the residual gauge symmetry. 
As a remarkable fact, "abelian dominance", irrelevance of off- diagonal gluons, is nu- 
merically observed for the nonperturbative QCD phenomena like confinement and 
dynamical chiral-symmetry breaking in the lattice QCD simulation in the MA gauge, 
which is a sort of the abelian gauge. Monopole condensation has been also suggested 
by the lattice QCD as the appearance of the global network of the monopole current 
in the MA gauge. 

First, we have studied the origin of abelian dominance for the confinement force 
in MA gauge in terms of the gluon-field properties using the lattice QCD. In the MA 
gauge, the gluon-field fluctuation is maximally concentrated in the abelian sector. 
As the remarkable feature in the MA gauge, we have found that the amplitude of 
the off-diagonal gluon is strongly suppressed, and therefore the phase variable of 
the off-diagonal (charged) gluon tends to be random, according to the weakness of 
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the constraint from the QCD action. Using the random-variable approximation for 
the charged-gluon phase variable, we have found the perimeter law of the charged- 
gluon contribution to the Wilson loop and have proved abelian dominance for the 
string tension in the semi-analytical manner. These theoretical results have been 
also numerically confirmed using the lattice QCD simulation. 

Second, we have studied th QCD-monopole appearing in the abelian sector in the 
abelian gauge. The appearance of monopoles have been transparently formulated 
in terms of the gauge connection, and is originated from the singular nonabelian 
gauge transformation to realize the abelian gauge. We have investigated the gluon 
field around the monopole in the lattice QCD. The QCD-monopole carries a large 
fluctuation of the gluon field and provides a large abelian action of QCD. Never- 
theless, QCD-monopoles can appear in QCD without large cost of the QCD action, 
due the large cancellation between the abelian and off-diagonal parts of the QCD 
action density around the monopole. We have derived a simple relation between the 
confinement force and the monopole density by idealizing the monopole contribution 
to the Wilson loop. 

Third, we have studied the monopole- current dynamics using the infrared monopole- 
current action defined on a lattice. We have adopted the local current action, consid- 
ering the infrared screening of the inter-monopole interaction due to the dual Higgs 
mechanism. When the monopole self-energy a is smaller then ac = ln{2D — 1), 
monopole condensation can be analytically shown, and we have found this system 
being the confinement phase from the Wilson loop analysis. By comparing the lat- 
tice QCD with the monopole-current system, the QCD vacuum has been found to 
corresponds to the monopole-condensed phase in the infrared scale. We have con- 
sidered the derivation of the DGL theory from the monopole ensemble, which would 
be essence of the QCD vacuum in the MA gauge because of abelian dominance and 
monopole dominance. 

In the second half part of this paper, we have studied the deconfinement phase 
transition and the hadron flux-tube system using the DGL theory. Deconfinement 
phase transition is one of the most interesting subject in the QCD phase transition 
both for ultra-relativistic heavy-ion collisions and for early universe. We have for- 
mulated the effective potential in the DGL theory, introducing the quadratic source 
term to study the QCD vacuum at finite temperatures. We have found the re- 
duction of the monopole condensate at finite temperatures, and found a first-order 
deconfinement phase transition at the critical temperature Tq — 0.49GeV using the 
temperature-independent parameters. The monopole condensate vanishes and the 
broken dual gauge symmetry is restored above Tc- We have considered the tem- 
perature dependence of the monopole self-interaction noting Tc = 0.28GeV as the 
lattice QCD simulation indicates. We have found a large reduction of the monopole 
self-interaction near the critical temperature. We have investigated the temperature 
dependence of the masses for the relevant mode like niB and m^, and have found 
their large reduction near the critical temperature Tq: mB,my ~ Tc- We have 
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calculated also the string tension at finite temperatures. The results agree with 
the lattice QCD data both in the variable A(T) and in the variable v{T) cases. In 
particular, the glueball mass reduction at high temperatures would be an important 
ingredient in the QCD phase transition. In the pure gauge, there are only glue- 
ball excitations with the large mass (~ IGeV) at low temperatures, and therefore 
it seems unnatural that the QCD phase transition takes place at a small critical 
temperature, Tq — 0.28GeV. This problem would be explained by the large reduc- 
tion of the glueball mass near the critical temperature as is demonstrated. In other 
words, this glueball-mass reduction may determine the magnitude of the critical 
temperature Tc in the QCD phase transition. 

Based on this effective potential at finite temperature, we have investigated the 
properties of hadron bubbles created in the early Universe and discussed the hadron 
bubble formation process. From numerical results, we can imagine how the QCD 
phase transition happens in the big bang scenario, (a) Slightly below T^, only large 
hadron bubbles appear, but the creation rate is quite small, (b) As temperature is 
lowered by the expansion of the Universe, smaller bubbles are created with much 
formation rate. During this process, the created hadron bubbles expand by radiating 
shock wave, which reheats the QGP phase around them.(c) Near Ti^^, many small 
hadron bubbles are violently created in the unaffected region free from the shock 
wave, (d) The QGP phase pressed by the hadron phase is isolated as high-density 
QGP bubbles,which provide the baryon density fluctuation. Thus, the numerical sim- 
ulation using the DGL theory would tell how the hadron bubbles appear and evolve 
quantitatively in the early Universe. 

We have studied the multi-flux-tube system in terms of the DGL theory. We 
have considered two ideal cases, where the directions of all the color-fiux-tubes 
are the same in one case and alternative in the other case for neighboring flux- 
tubes. We have formulated the system of multi color-fiux-tubes by regarding it as 
the system of two color-fiux-tubes penetrating through a two dimensional sphere 
surface. We have found the multi flux-tube configuration becomes uniform above 
some critical flux-tube number density pc = 1.3 ~ 1.7fm~^. On the other hand, the 
inhomogeneity on the color electric distribution appears when the flux-tube density 
is smaller than pc. We have discussed the relation between the inhomogeneity in the 
color-electric distribution and the flux-tube number density in the multi-flux-tube 
system created during the QGP formation process in the ultra-relativistic heavy- 
ion collision. When the flux-tube number density is low enough, the system can 
be approximated as the incoherent sum of the individual flux tube. Its evolution 
would lead to the creation of q-q pairs via flux-tube breaking. In this process, this 
q-q creation plays a essential role on the QGP formation. On the other hand, for 
a dense flux-tube system, neighboring flux tubes are melted into a large cylindrical 
tube, which is approximately homogeneous in the whole space and is regarded as 
the "color condenser". In this case, large homogeneous QGP may be created in 
the central region. Thus, the process of QGP formation is expected depending 
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largely on the density of created flux tubes. We have finally investigated also the 
interaction of two color-flux-tubes in the DGL theory and showed that there is a 
possibility of annihilation or unification of two color-flux-tubes when the flux tubes 
are sufficiently dense. This process would provide an enough energy to create many 
dynamical gluons and to form QGP. 
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Appendix A 

Monte Carlo Method for Lattice 
QCD 

A.l Gauge Field on the Lattice 

Quantum Chroinodynamics (QCD), a nonabelian gauge theory, is difficult to solve 
analytically due to the large gauge coupling in the low-energy region. The Monte 
Carlo simulation based on the lattice gauge theory is one of the promising methods 
for the direct calculation of the QCD partition functional. We review here the 
fundamentals of the lattice QCD and the numerical simulation [oU . 

The Minskowski space-time is transformed into Euclidean space-time by replac- 
ing x° = —ix"^, where x° and x'^ are the Minskowski and the Euclidean time, respec- 
tively. A similar transformation is made in the time component of all four-vectors. 
In the Euclidean space, there is no distinction between upper and lower indices of 
four-vector, i.e, x^ = (x^, x^, x^,x^) = {xi,X2,X3,X4). 

In principle, to latticize the theory is performed by the replacement of derivatives 
by finite differences with the lattice spacing a. This method, however, does not 
preserve gauge invariance, which is an essential attribute of the theory. Wilson 
introduced a gauge- invariant lattice action using the path representation of the gauge 
group G. The most elementary paths on the lattice are links out of which any path 
can be constructed. The link is identified by the site s and the direction fi. With 
each link on the lattice, we associate a group element, called the "link variable" as 

U^{s) = exp{iaeAf,{s)} = exp{i9^{s)} e G, (A.l) 

with the gauge coupling constant e. The lattice angle variables 6'^ defined as 6^ = 
mrpa ^ aeAfj^ is dimensionless, and link variable Uf^{s) is transformed as 

U,{s)^V{s)U,is)V^is + fi), (A.2) 

where V{s) and V{s + ft) are the gauge functions located at the starting and the 
end points of the link U^{s). Thus, the gauge transformation is described by simple 
multiplication of the group elements in the lattice formalism. 
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In the continuum, the field strength tensor 0^,^ = G" T" E g is defined in terms 
of an infinitesimal closed path. By analogy, we defined it on the lattice in terms of 
a "plaquette" n^j,, a square bounded by four links. The field tensor Gfj,^, which is 
associated with the oriented plaquette specified by the links {/x, z/} attached to site 
s, is defined through the relation 

n^,(s) = U^{s)U,{s + fi)Ul{s + i))Ul{s) = exp{tea^G^,{s)) G G. (A.3) 

For a ^ 0, we recover the continuum field tensor, 

G'^.(s) = -[(Ms + /i) - A,{s)) - {A^{s + z>) - A^{s))] + ie[A^{s), A,{s)] + 0{a). 

The lattice action is constructed by using the relation, 

^^trn^,(s) = EEtr[l - ^ea'G,, - \e\\G,,f + ■■■], (A.4) 

S fl,U S fl,U 

where n and u are summed from 1 to 4. The first term on the right hand side is a 
constant, and the second term vanishes because Gf^u is traceless. The sum over s, 
H, V on the left hand side is twice the sum over all plaquettes, so that one finds 

^ ^ trD^, = Y.Y. (trn^. + trniJ = 2Re ^ J] trD^,. (A.5) 

Because interchanging fi and u changes the orientation of the designated plaquette. 
Thus, we have 

^ ^ Re trD^, ^ -le' j d^xY,ii{G^,f + const., (A.6) 

S fl>U ^,u 

as a -^ 0. Here, the sum on the left hand side extends over all plaquettes on the 
lattice. 

Finally, we get the lattice QCD action as 



S = /rf'^E^tr(G,,r = ^EERetr[l-n,,(s)] 



= /5EERei^tr[l-n,,(s)], (A.7) 

where (3 = ^^ is the control parameter relating to the lattice spacing a. The 
continuum limit a ^ corresponds to e — > or /5 — ^ oo in the theory with the 
asymptotic freedom like QCD. Thus, the QCD partition functional in the lattice 
formalism is given as 

rff/^e-5[t/.W] = /[n,,^rff/^(s)]e-^[^-(^)]. (A.8) 
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For SU(2), the hnk variable U^{s) is parameterized as U^ = f/° + iT"'U°^, and the 
measure takes a form 

dU,{s) ^ dUldUldUldUpiUl + Ul + Ul + Ul - 1). (A.9) 

The expectation value of arbitrary operator O is calculable as 

using the numerical simulation with the Monte Carlo method. 



A.2 Monte Carlo Method 

In this section, we briefly summarized the Monte Carlo method of the lattice QCD. 



In the lattice QCD with the action (|A.7| ), the ensemble of gauge configuration is 



given by the canonical ensemble, characterized by Boltzmann distribution, 

P{U) = |e-^^ (A.ll) 

where we put S\U] = l3S[U] and regard S[U] as a "Hamiltonian" . In principle, the 
lattice QCD partition function Z can be calculated by the infinite number of Monte 
Carlo simulations using any random sampling of U^. In practical, however, almost 
all gauge configurations with large S do not contribute to the lattice QCD partition 
functional Z. Therefore, it is necessary to sample the "important gauge configu- 
ration" with small 5* effectively in terms of the numerical calculation with finite 
number of operations. Such a "important sampling" can be realized by generating 
the random number with the weight 6"^*^^^', and the Monte Carlo method enables 
us to calculate the lattice QCD partition function, numerically. 

In the Monte Carlo method, our goal is to generate a time sequence of configura- 
tions such that the configuration U occurs with probability P{U) after a sufficiently 
long time. Thus, the time average of any quantity would be the same as its average 
over the canonical ensemble. 

The time sequence is generated by a stochastic process. A configuration U is 
updated to U' with a transition probabihty T{U', U), which has the following general 
properties: 

1. T{U',U) >0 

2. JdU'T{U',U) = JdUT{U,U') = 1, 

3. T{U', U)P{U)dU = T{U, U')P{U')dU'. 
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The first two are just properties that any probabihty should have. The last is 
known as the condition of detailed balance. Integrating over U' on both sides of the 
last condition, we obtain 

P{U)= j dU'T{U,U')P{U'), (A.12) 

which states that the probability P{U) is an eigenfunction of T(f/', U). This means 
that transition probability preserves the equilibrium ensemble. Thus, we get the 
gauge configuration ensemble obeying P{U) as a "time" sequence of "thermaliza- 
tion". Using the obtained gauge configuration, the ensemble average of arbitrary 
operator 0{U) is numerically estimated by 

{O) = j dUO{U)P{U). (A.13) 

On the actual numerical simulation, there are two cautions to be carefully checked. 
One is the achievement of the thermo-equilibrium, the other is the vanishing of the 
cancellation among the sampling gauge configurations. 



Appendix B 

Procedure of Maximally Abelian 
Gauge Fixing 



The maximally abelian (MA) gauge is the special abelian gauge exhibiting infrared 
abelian dominance in the lattice QCD, and provides the theoretical basis of dual 
Higgs picture from QCD. In the SU(2) lattice formalism, the MA gauge is defined 
so as to maximize 

Rma[u,] = Wnis) 

^ s 

^ J2tT{U,{s)rsUl{s)rs} 

S,IJ, 
S,fJ, 

= 2Y^[l-2{Ul{sf + Ul{sr}] (B.l) 

S,IJ, 

by the gauge transformation. Here, R{s) is a local scalar variable defined as 

R{s) = J2tT{U^,{s)rsUUs)rs} with U^^{s) = U^s - fi). (B.2) 

Here, R{s) is manifestly invariant under the lattice rotation and the reflection. In 
the MA gauge, the operator 

^s)^J2U±,is)TsUUs) (B.3) 

is diagonalized. In this appendix, we show the procedure of the MA gauge fixing on 
the lattice. 

To begin with, we introduce a "local" gauge transformation, whose gauge func- 
tion Ko(s) is not unity at the site Sq only, 

f Ko(s) = V{so) fors = so /g^N 

1 Ko(s) = 1 for s^ So- 
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In order to maximize the value Rma[U^], one may consider to maximize the local 
variables R{s) at each site s by the local gauge transformation Vg. However, since 
i?(so) at the site Sq is changed not only by the gauge transformation K(,(<s) but also 
by the gauge transformation VsQ-fi{s) with neighboring sites so—fi, one cannot obtain 
the MA gauge configuration Uf^i^so)^^ only by simple local gauge transformation 
Vsy{s). After the local gauge transformation at all sites on the whole lattice, one 
has to repeat this procedure until Rma is maximized. 

Now, let us derive the gauge transformation K()(s) to maximize -R(so). After the 
gauge transformation by Vs^Xs), -R(so) is changed as 

R^'iso) = Y.^T{V{so)U,{so)rsUl{so)V^{so)r, 

+U^{so - il)V\so)TsViso)U\so - il)T,} 
= tT[J2{U^iso)TsU\so) + Uliso - il)rsU^iso - fi)} ■ V\so)r3V (so)] 

= tTMso)S{so)]. (B.5) 

Here, we define 

S{s) = 5'^(s)r'^ = S-T = V\s)tsV{s) e su(2), (B.6) 

$(s) ^ ^-^s)t^ ^ $ . ^ ^ ^ U^,{s)rsUUs) G su(2), (B.7) 

which are both elements of Lie algebra and satisfy relations tr($) = tr(S') = and 
S"^ = 1. To maximize R{sq) by this gauge transformation, S{so) is taken to be 
the same direction as $(so) in the SU(2) ~ 0(3) parameter space, 5*//$. After 
this gauge transformation, $(so) is diagonalized as <l>^(so) = V {so)^{sq)V'^ (sq) = 
$^(so)t3, and S{sq) becomes T3. Here, $ plays a similar role as the Higgs field in 
the 't Hooft-Polyakov monopole. 

In the abelian gauge, the gauge function V{s) is an element of the coset space 
SU(2)/U(1)3 using the residual U(l)3 degrees of freedom. We take the representative 
element of V{s) so as to satisfy V^{s) = 0, or 

V{s) = V\s) + i{V\sy + V^sy}. (B.8) 

Because of (V^)^ + (V^)^ + {V^Y = 1, we can parameterize V{s) as 

= cos 6'(s) 

= sin^(s)cos0(s) (B.9) 

= sin^(s) sin0(s), 

and then S{s) is expressed as 

S = S^t" = VW3V = sin2^cos0ri +sin2^sin0r2 + cos2^r3. (B.IO) 




Figure B.l: The gauge function V^ver used in the over-relaxation method. The 
vector Vtarget correspouds to the gauge function which maximizes R{so)- 

Since $(so) = $"(so)t'^ is obtained from the original gauge configuration U^{s), we 
get Ko(so) as 

tan Ati[So) - (^3^2 - (^3)2 fB.ll) 



tan0(so) 



52 ^ ^ 
51 $1- 



Thus, the gauge function V{s) which maximizes R{so) is obtained so as to obey 
S'(so)//$(so). This procedure makes -R(so) defined in (|B.5D maximum by Vsq{s). 
This gauge transformation, however, influences R{s) of the neighboring sites, s = 
So ± ft, and in fact does not make them maximum. Therefore, we have to perform 
this procedure to the neighboring sites. By doing this, however, the original R{so) 
gets some change and hence R{so) is no more in its maximum. This fact forces us 
to repeat the local-gauge transformation many times. 

To optimize the convergence, in the practical simulation, we take an over-relaxation 
method. We show the vector (cos 6, sin 6 cos 0, sin 6 sin 0) in Fig. |B.l| , corresponding 



to the gauge function V{so) in Eq.( p.9|) . In the over-relaxation method, we take the 
angle value Q9 instead of 9 obtained in Eq.( |B.9D , 



[ K°er(50) = 


= cos(fie) 


Kier(So) = 


= sin(r26') cos 


I VLriSo) - 


= sin(r2^) sin 



(B.12) 
This overrelaxation parameter Q is taken as 1 ~ 2. 
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